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Introduction 

The NAVSTAR Global Positioning System (GPS) satellites are being used in a wide range of 
civilian applications besides their intended use for navigation by the Department of Defense’s 
military applications. Applications can be divided into two broad categories, namely, navigation and 
static positioning (absolute or relative). Knowledge of the position of the satellites is fundamental to 
all the applications of GPS. The required knowledge of the GPS satellite position accuracy can vary 
depending on a particular application. Application to relative positioning of receiver locations on 
the ground to infer Earth’s tectonic plate motion requires the most accurate knowledge of the GPS 
satellite orbits. Research conducted at the University of Texas Center for Space Research (CSR) 
under the NASA grant NAG5-940 was directed towards improving and evaluating the accuracy of 
GPS satellite orbits. This reports summarizes the activities of CSR in improving GPS satellite 
orbits, during the grant period from Jan 1988 to April 1991. Major focus was on the understanding 
and modeling of the forces acting on the satellites besides paying close attention to reference frame, 
time system, measurement modeling, parameterization and other aspects of orbit determination. 
Gravitational forces were modeled by truncated versions of extant gravity fields such as, GEM-L2, 
GEM-T1 or TEG-2 and third body perturbations due to Sun and Moon. Nongravitational forces 
considered were the solar radiation pressure, and perturbations due to thermal venting and thermal 
imbalance. At the GPS satellite orbit accuracy level required for crustal dynamic applications, 
models for the nongravitational perturbations play a critical role, since the gravitational forces are 
well understood and are modeled adequately for GPS satellite orbits. With these considerations, 



CSR’s activities during the grant period were in the following areas: 

1. Computation and distribution of weekly GPS satellite ephemerides 

2. Global versus regional orbit comparison 

3. Signal multipath effects 

4. CASA UNO campaign data analysis 

5. Southwest Pacific campaign data analysis 

6. Reference frame comparison 

7. Modeling of thermal imbalance force perturbation 

8. Eclipsing satellite orbit prediction error 

9. Effects of mixed receiver types 

10. Effects of selective availability (SA) 

In the following, each of the above topics will be described briefly, and relevant publication/meeting 
presentation will be quoted. 

1. Weekly Ephemerides 

Pseudo-range measurements from all the GPS satellites collected at the CIGNET sites (three in 
late 1986 expanding to nine by 1988) were processed routinely to produce precise ephemerides of 
the satellites. The ephemerides were written in the ’NGS Format* and were distributed to the 
participants of data collection and the sponsors of the network, namely, NGS at Rockville, 
Maryland, USGS at Menlo Park, California, Texas State Highway Department at Austin, Texas, 
CGS at Ottawa, Canada, (through NGS to) Institut fuer Angewandte Geodaesie (IFAG), Germany, 
and the Norwegian Geodetic Survey at Tromso, Norway. Apart from regular weekly distribution of 
the GPS ephemerides to the above agencies, from time to time requests for precise ephemerides for 
specific GPS weeks by other institutes in the U.S. (CCAR, UNAVCO etc.), England and Australia 
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were fulfilled. Precise ephemerides computed using GPS carrier phase measurements were 
distributed for specific weeks covering GPS campaigns. For examples of characteristics of orbits 
distributed see Appendix I. 

2. Comparison of Orbits Obtained Using Data From Global and Regional Network 

The Cooperative International GPS Network (CIGNET) started with three ground stations in the 
continental U.S. and with time expanded in the U.S. and to Europe, Japan and Australia. Since the 
accuracy of GPS satellite ephemerides is fundamental to baseline computation applications, effects 
of changing the ground station network distribution on the orbit accuracy is an important aspect to be 
considered. At the time of this analysis the distribution of stations contributing data for routine GPS 
orbit computation was not really global, but were restricted to U.S. and northern Europe. Five 
stations in the north American continent (Yellowknife, Mojave, Westford, Richmond, Austin) and 
three in Europe (Wettzell, Onsala, Tromso) were considered in orbit determination experiments 
using only the pseudo-range observations. Solutions for several weeks were computed using data 
from the five stations in the North American continent only or from all the eight stations mentioned 
above. Trajectories were computed for individual GPS satellites using solutions of these two cases 
and compared in order to find any common characteristics that be present. Depending on the relative 
amounts of data and the geometry of coverage, the differences between the ephemeris determined 
using data from North American stations only and using data from all the eight stations, varied from 
about 4 meters to about 40 meters. Although these numbers are not the accuracies of the orbits, they 
are an indication of the effects of network distribution, geometry and the amount of data. Orbit 
accuracy itself could be inferred through baseline estimates or through comparison to orbits with 
known accuracy computed by other organizations. However, during the period of this analysis such 
evaluations were not attempted. Details of this study were reported in a Master’s Thesis by Vince 
Pollmeier, and the summary and conclusions of this thesis are presented in Appendix n. 
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3. Signal Multipath and Choke Ring Experiment 

Degradation of accuracy of the GPS measurements, namely, pseudo-range and carrier phase, due 
to signal multipath is an aspect that has been treated with considerable importance in the hardware 
(receiver and antenna) design and in the data analysis. Measurement accuracy directly impacts the 
orbit accuracy and hence every attempt is made to improve accuracy of GPS measurement. Signal 
multipath is a phenomenon in which the radio waves coming directly from the GPS satellite are 
corrupted by their reflections off of any nearby surfaces which can reflect radio waves. Since the P- 
code pseudo-ranges are computed by time correlation of the P-code signal with that generated at the 
receiver, its accuracy is affected by multipath much more than that of carrier phase which involves 
only the measuring of the phase of the carrier. The difference could be as high as two orders of 
magnitude. Although in geodetic applications only carrier phase is used as the measurement of 
choice due to its precision, the pseudo-range is essential to one or more of the following necessities 
of data processing, namely, phase ambiguity resolution, cycle-slip detection and fixing, measurement 
time-tag correction. Jn order for the pseudo-range to be effective in its above mentioned uses, its 
accuracy should be as high as possible; higher the accuracy, lesser would be be the averaging time 
needed. Pseudo-range accuracy level that is less than the wavelength of its carrier is desirable. For 
this reason several antenna designs have been studied and experimented by JPL. One of the 
experiments conducted by JPL was to use a choke ring antenna in place of the regular TI-4100 
antenna to collect data at the Mojave CIGNET site for a few days. This data set was analyzed at 
CSR, and some of the results were presented at the GPS Workshop held at JPL during the Spring of 
1988. A part of the presentation is included in Appendix EL CSR also cooperated with personnel at 
USGS Menlo Park, California, in their experiments and analysis of data to determine the effect of 
multipath error on accuracies of baseline determinations. The abstract of the paper presented by 
USGS personnel at the spring 1989 meeting of the American Geophysical Union at Baltimore, 
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Maryland, is included in Appendix in. 


4. CASA UNO Data Analysis 

During January and February of 1988, the Jet Propulsion Laboratory, in cooperation with other 
agencies (national and international), conducted an epoch measurement experiment in South 
America and in South Pacific. Data were collected from stations in North and South America, 
Hawaii, South Pacific and Australia in conjunction with the CIGNET stations (in terms of tracking 
scenario and collection interval). The campaign period was about three weeks during the ’GPS 
Weeks’ 419, 420 and 421. Precise GPS ephemerides using pseudo-range and phase measurements 
and baselines between various sites were determined using this data set. Extensive analysis has been 
performed at CSR using this data set over the past two years concentrating on various aspects of 
determining precise GPS satellite orbits and accurate baselines between tracking sites. Force model 
fidelity, numerical integration accuracy, proper reference frame for geodetic applications of GPS and 
estimation strategies were some of the topics that were scrutinized closely, using this and other data 
sets. Results of analysis of which this data set formed a major part were published in various forms 
(meeting presentations, journal article, thesis, dissertation) part or whole of which are included in 
Appendix IV. 

5. Southwest Pacific Campaign Data Analysis 

Tonga Trench in the southwest Pacific region is considered to be one of the most active fault 
zones in the world, and GPS campaigns (sponsored jointly by NSF, NASA and UNAVCO) were 
conducted during three consecutive summers of 1988, 1989 and 1990, to study plate motion around 
this zone. CSR participated in the campaign planning and preparation and in data processing. 
Precise GPS ephemerides were computed for selected weeks during each of these campaign periods. 
These include weeks 444,445 in 1988, weeks 498 to 501 in 1989 and weeks 549 to 551 in 1990. 


rr.r* 'sr , Kr_'-r.??3*r'T<r 
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Extensive analysis of some of the baselines in this region and evaluation of GPS satellite orbit 
accuracy have been performed using these data sets and are being continued. Although neither SLR 
nor VLBI measurements of baselines in this region is available for comparison with GPS 
determinations, orbit accuracies can be inferred by considering the baseline repeatabilities and other 
criteria. Results of experiments performed at CSR using these data sets were presented in various 
meetings and have published in symposium proceedings. Copies of some of the publications are 
included in Appendix V. 

6. Reference Frame Comparisons 

Most of the GPS applications including some researchers involved in geodetic applications use 
the GPS ephemerides broadcast from the satellites as their source. Consequently, the solutions 
(absolute or relative positions) are defined in the coordinate system of the GPS ephemerides used 
which is the WGS-84. Often such results (especially geodetic) are compared with solutions obtained 
in other coordinate systems such as defined by SLR or VLBI. These coordinate systems should be 
coincident for comparisons to be meaningful or their relative orientation must be known in order to 
account for differences that may be present in the results. In order to facilitate such comparisons or 
to transform solution in one system to other a set of seven parameter transformation was estimated. 
This was performed using two sets of tracking data. 

1. The first experiment used tracking data from the five Operational Control Segment (OCS) 
stations and data from CIGNET sites collected during the CASA UNO campaign (weeks 419 
to 421). The data type used in the analysis was pseudo-range since that was the only type 
available from the OCS stations. The coordinates of the OCS stations were available in the 
WGS-84 reference frame whereas those of the CIGNET sites were available in a reference 
frame defined by SLR and VLBI systems. Holding few of the CIGNET sites fixed, 
coordinates of the OCS sites (along with other parameters) were estimated which were in the 
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SLR/VLBI system. Comparing these to the coordinates of the OCS sites given in WGS-84 
system, the seven parameter transformation set was determined using least-squares technique. 

2. This same procedure was repeated in a second experiment, but using carrier phase GPS data 
from the five Defense Mapping Agency (DMA) sites whose coordinates were available in the 
WGS-84 system, and from the CIGNET sites augmented by a global set of stations. The data 
used in this experiment were part of those collected during 1989 south west pacific campaign, 
namely, the weeks of 499 and 500. 

Results of these experiments were presented in various meetings and published in symposium 
proceedings. Abstracts/summaries of these publications are provided in Appendix VI. 

7. Thermal Imbalance Perturbation 

When the surfaces of the GPS satellite directly exposed to the Sun are heated up, most of the 
energy is radiated back into the space. The surface away from the Sun also gains temperature due to 
conduction and then radiates the energy away into space. Depending on the conductive properties of 
material between the surfaces and the radiation characteristics of both the surfaces the satellite as a 
whole may or may not be in thermal equilibrium and as a result may experience a small force if there 
is a thermal imbalance. During the transient phase when the satellite enters or exits the Earth’s 
shadow, a thermal nonequilibrium and hence a net force is certainly possible. Depending on the 
characteristics of the transient phase (rate of heating/cooling) and the duration of the shadow period, 
the magnitude of this perturbation could be significant. Simulation experiments were performed at 
CSR to determine effects of such thermal imbalance perturbations on satellite orbits assuming such 
simple models as flat plates or rectangular boxes for the satellites. Such simple ’satellites’ were 
assumed to be in orbits of full Sun or maximum shadow configurations and the heat equation 
(second order partial differential equation) was considered in one dimension only. The major 
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difficulty in including this perturbation in oibit computation is due to the fact that the second order 
ordinary and the partial differential equation (finite difference or finite element techniques) have to 
be integrated simultaneously. Some of these aspects were studied closely but complete 
implementation of this perturbation for a typical satellite has not been achieved (mainly due to the 
difficulty mentioned above). Results of this study have been reported in Masters Thesis and in 
presentations at conferences. Appropriate summaries are presented in Appendix VII. 

8. Eclipsing Orbits 

Solar radiation pressure is one the most important nongravitational perturbations acting on the 
GPS satellites. Because of the altitude of the satellite, the magnitude of this perturbation is high 
(being of the order of 10~ 7 m/sec/sec) and hence play a major role in the accurate computation of the 
GPS satellite orbits. Its modeling is critical to achieving orbit accuracies needed for geodetic 
applications. Acceleration due to this perturbation is computed using the "ROCK4" models when 
the satellite is exposed to the Sun and is assumed to be zero when the satellite is in the shadow of the 
Earth (or the Moon). Thus when the orbit geometry is in such eclipsing configuration, computed 
acceleration due to SRP has a large jump in its time history, at the shadow crossing points (entering 
and leaving the shadow). Such large discontinuities in the acceleration values cause considerable 
error in the values of satellite state computed using constant step numerical integrators such as the 
KSG method. Manifestation of this problem has been clearly identified in situations where GPS 
satellite orbit predictions were performed in order to evaluate ephemeris accuracy. Variable step 
integrators such as the Runge-Kutta method are computationally expensive and need restarting of the 
numerical integration algorithm at such discontinuities. The approach taken to remedy this 
deficiency is to perform special operations at times when the satellite enters or leaves the shadow. 
This point usually occurs within the interval of a numerical integration step and hence the 
acceleration function values used by the integrator at the ends of this step have a large jump. To 
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avoid this discontinuity, the function values are assumed to be valid beyond the shadow crossing, on 
either side of it. In other words, when the satellite enters the shadow, the SRP acceleration value 
computed before shadow is assumed to be valid up to the first nodal point (of the integrator) after 
shadow boundary, and SRP acceleration value of zero is assumed backward up to the first nodal 
point before the shadow boundary. With this assumption the "back-difference table" of the KSG 
numerical integration algorithm is modified and appropriate corrections to the satellite state at these 
nodal points are made. This procedure is called the "modified back-difference table" (MBD) 
approach. Information on more details about this approach and on the improvements obtained in the 
prediction accuracy of satellite orbit can be found in Appendix VIE and the references thereof. 

9. Effect of Mixed Receiver Types 

The earlier tracking sites of the Cooperative International GPS Network (CIGNET), namely, 
Mojave, Westford, Austin, Wettzell and Richmond, used the TI-4100 (dual frequency P-code) 
receivers. Later in 1989, additional CIGNET sites collected data using Mini-Mac 28 16- AT (dual 
frequency codeless) receivers. Further, NGS started replacing the TI-4100 receivers with the Mini- 
Macs. As a result, for a period of time data from two types of receivers had to be processed 
simultaneously for orbit determination. The differences in the data types and especially in way the 
time tag of observations were handled in both of these receivers had a considerable impact on 
preprocessing procedures and on the accuracy of the estimated orbits. In order to evaluate the 
impact of such mixed receiver data type on orbit accuracy a limited experiment was conducted. This 
was considered important since some of the field campaigns (such as the South West Pacific 
Campaign) used more than two types of receivers to collect data and if orbit accuracy is affected due 
to mixed receiver data, it could impact the baseline estimates. 

During specific weeks at Mojave and at Wettzell data were collected, as part of CIGNET 
tracking, simultaneously from both the TI-4100 and the Mini-Mac receivers. Several orbit 
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determination experiments were performed using these data sets. In one experiment, only TI-4100 
receiver data collected at Mojave would be used along with data from other sites (all TI) Mojave 
would be replaced by Mini-Mac data collected at Mojave during the same period, thus providing a 
mixed data set. Orbits determined in these two cases were compared and in addition were used in 
determining a known baseline. The rms of orbit differences was as much as 4 meters for some 
satellites; however, the impact on accuracy of the particular baseline tested, was not significant. 
Results of this study were presented in meetings and were published in meeting proceedings, copies 
of which are included in Appendix IX. 

10. Effects of Selective Availability (SA) 

The influence of Selective Availability (SA) on geodetic positioning was presented at the Fall 
1990 meeting of the American Geophysical Union ("Observations of SA and the Effect on Precision 
Geodesy"). An example of SA is given in Appendix X, which illustrates phase residuals obtained by 
removing the effect of the orbit, the ionosphere and the troposphere. The examples shown were 
taken simultaneously from different locations and receivers. The similarity of The residuals indicates 
that the source of the variation was common and not produced by local effects, such as clock 
variations. In fact, the effect is SA and is characterized by a spectrum of frequencies with periods of 
several minutes. It was also demonstrated that SA is common to LI and L2; thus the effect cancels 
in LI minus L2 differencing. Double difference phase residuals on a zero baseline, i.e., Trimble 
4000SST and TI4100 receivers were operated with the same antenna. During the period of the 
experiment. Week 549, SA was active on Block II satellites. As a consequence, double differences 
formed with the zero baseline using only Block I satellites can be characterized by white noise 
residuals. However, because of the data recording offset between the Trimble and TI receivers (the 
TI recorded 0.92 sec prior to the Trimble), the effect of SA does not cancel in the differencing 
process, and residuals that exhibit SA-induced periodicities remain. Other results, using baselines of 
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1000 km, showed that if the receivers record at nearly the same time (to within a few milliseconds), 
the effect of SA can essentially be removed. 



APPENDIX I 


1. Report on GPS orbit determination for Week 378 using pseudo-range data 

2. Report on GPS orbit determination for Week 378 using double differenced phase data 


REPORT ON GPS ORBIT DETERMINATION 


GPS WEEK 378 

APRIL 5 TO APRIL 11, 1987 
(DAYS 95 - 101 ) 

RELEASE: 

EW378VA . EPH 


PREPARED 

FOR 

U. S. GEOLOGICAL SURVEY 
MENLO PARK, CA 

BY 

CENTER FOR SPACE RESEARCH 
UNIVERSITY OF TEXAS AT AUSTIN 

WITH THE COOPERATION 
OF 

TEXAS STATE DEPARTMENT OF HIGHWAYS AND PUBLIC TRANSPORTATION 

AUSTIN, TX 

NATIONAL GEODETIC SURVEY 
ROCKVILLE, MD 

CANADIAN GEODETIC SURVEY 
OTTAWA, CANADA 


For further information: 

Dr. B. E. Schutz or Dr. P. Abusali 
Center for Space Research 
The University of Texas at Austin 
Austin, Texas 78712 


i-iRPANET : 

SCHUT2V.UTCSR . DECNET8CHPC . BRC . UTEXAS . EDU 


Telephone: 

(512) A71-A267 
FTS 770-52A0 
Telex : 

70A265 CSRUTX UD 
BITNET: 

A0VS751 S> UTCHPC 
SPAN: 

UTSPAN : : UTCSR : : SCHUTZ 


GPS Week 


378 


Contributing Stations : 


Locat ion 


ID 


Receiver 


Reference Agency 

Qsci 1 lator 


Austin, TX 
Kauai , HI 
Mojave, CA 
Westford, MA 
Yellowknife, NWT 


A TI-4100 
K TI— 4100 
M TI-4100 
W TI-4100 
Y T I -4 100 


Cesium 
H-Maser 
H-Maser 
H-Maser 
Ces i um 


SDHPT 

UNAVCO 

CSR 

NGS 

CGS 


NOTE: SDHPT - State Department of Highways and Public Transportation (Texas) 

CSR - Center for Space Research 

NGS — National Geodetic Survey 

CGS - Canadian Geodetic Survey 

USGS — U. S. Geological Survey 

UNAVCO- University Navstar Consortium 
MODELS : ^ 

Gravi ty 

MERIT Standard ( GEM-L2 , modified) 

Luni-solar from DE-200 


Nongravi tat ional 

Constant area, estimate reflectivity, estimate y-bias 
Polar Motion, UT1 

Lageos derived, reported on GE Mark III 


Station Coordinates 


Location 

Aust i n : 

Kauai : 
Mojave: 
Westford : 

Yel lowknif e : 


x ( m ) 

-743773.949 
-5543S1 7 . 634 
-2356214.421 
1492233.223 
-1224064.162 


y <m) 

-5460643.710 

-2054588.618 

-4646733.994 

-4458091.535 

-2689834.572 


z ( m ) 

3200347.492 

2387854.437 

3668460.630 

4296045.911 

5633432.720 


Reference 

Mader Cl 987b 1 
Mader C 1987a] 
Mader C 1987b] 
Mader C 1987a] 
McArthur C 1937 


The Conventional Terrestrial System is defined by the NGS VLBl-defioed 
frame. 

References 

Mader, G. L., "GPS Antenna Connections to VLBI Reference Points at 
Kokee Park and Westford," NGS Memo, July 17, 1987a. 

Mader, G. L., "GPS Antenna Connection to VLBI Reference Point at Mojav 
and the Position of Austin GPS," NGS Memo, October 13, 1987b. 

McArthur, D., "Yellowknife GPS Monitor Site," CGS Memo, February 4, 19 


Data Distribution : 


DAY 3 

A 

K 

Sunday M X 

Y X 

A 

K 

Monday M X 

Y X 

A 

K 

Tuesday M X 

Y X 

A 

K 

Wednesday M X 

Y X 

A 

K 

Thursday M X 

Y X 

A 

K 

Friday M 

Y X 

A 

K 

Saturday M 

Y X 


Satel 1 i te PRN 
6 9 11 12 13 

X X X X X 

X X X X X 

X X X X X 

X X X X X 

X X X X X 

X X XXX 

x X X X X 

X X X X 

X X X X 

X X x.x X 

X X X X X 

X X X X X 

X X X ' X X 

X X X X X 

X X X X X 

X X X X X 

X X X X X 

X X X X X 

X X X X X 

X X X X X 

X X X X X 

X X X X X 

X X X X X 

X X X X X 


X X X X x 


NOTES: 

o The locations , e.g.« M» are defined on the preceding page 
in the column labelled ID. 

o Data from the Westford receiver suggests that the receiver Mas 
operated on the internal oscillator instead of the hydrogen 
maser; consequently, Westford data were not used in this ephemeris 
release. 



SOLUTION CHARACTERISTICS (NORMAL POINT DATA): 


Sa tell 
PRN 

i te 

A 

RMS ( m ) 
K M 

Y 

Over - 
all 

Sig 

A 

NP Meas. Processed 

K M Y TOTAL 

3 

- 

- 

0.94 

1 .26 

1 . 18 

11 

- 

- 

55 

142 

197 

6 

1.85 

1.21 

1.05 

1 .36 

1.42 

3 

i 17 

84 

86 

195 

<T62 

9 

2.29 

0.90 

1 . 18 

1.38 

1.58 

2 

172 

150 

96 

198 

616 

1 1 

1 .51 

1 . 17 

1.01 

1.27 

1.30 

4 

182 

71 

99 

220 

572 

IE 

2.29 

1.24 

1 .54 

1.35 

1.75 

3 

169 

70 

78 

178 

A95 

13 

2.06 

1 .36 

1 .58 

1 .30 

1.63 

3 

188 

96 

121 

207 

612 


NOTES: The column labelled M Sig" is the formal standard error based on 
the trace of the satellite position estimate variances? using 
a measurement standard deviation of 1 m. It reflects the 
tracking geometry and the amount of data available. 

The normal point <NP) data were obtained from the raw data by 
compressing the 30 second interval data into 15 minute bins. 


STATE COMPARISON: 


Th.e following table provides a comparison in the epoch positition vector. 
The differences in position are obtained from the position vector reported 
in the EW378V1 solution distributed earlier and the current solution, 
EW378V4. The values in the table correspond to differences at the initial 
epoch only. The primary differences in the two solutions are as follows: 

c Different station coordinates: analysis has shown difference 

in the two coordinate sets of less than two meters in 
individual components 

o Use of Kauai data instead of Haleakala: although the amount 

of data from the two Hawaiian sites is very similar, the Kaua 
coordinates are more internally consistent with the other 
coord i nates . 

Epoch: 5 April 1987 00:00:00 GPS 

Week 378 


Epoch Difference (m) . One Week Difference (m) 

Satellite (EW378V4 minus EW37BV1) (EW378V4 minus EW378V1 ) 


PRN 

X 

Y 

Z 

Mean 

RMS 

3 

o 

• 

10.9 

1S.0 

16.5 

17.0 

6 

-1.5 

-5.7 

8.4 

5.7 

6.0 

9 

-5. 1 

-4.7 

1.6 

6.S 

6.3 

11 

0 

a 

I s - 

H 

1 

15.4 

7.0 

18.7 

19.4 

IE 

-5.1 

-IE. 6 

-5.S 

15.1 

16.0 

13 

-5.4 

-0.6 

4.5 

4.3 

4 . E 
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GPS WEEK 378 

APRIL 5 TO APRIL 11, 1987 
(DAYS 95 - 101) 


DIFFERENCED PHASE 
RELEASE A 
(EW378D1.EPH) 


PREPARED 

BY 

CENTER FOR SPACE RESEARCH 
UNIVERSITY OF TEXAS AT AUSTIN 


Sponsored by: '■ 

U. S. Geological Survey 

National Aeronautics and Space Administration 
National Science Foundation 


For further information: 

Dr. B. E. Schutz or Dr. P. Abusali 
Center for Space Research 
The University of Texas at Austin 
Austin, Texas 78712 


Telephone : 

(512) 471-4267 
FTS 770-5240 
T© lex : 

704265 CSRUTX UD 
BITNET: 

SCHUTZ @ UTADNX 
SPAN: 

UTSPAN : : UTADNX : : SCHUTZ 


GPS Week : 378 

Contributing Stations 


Location 

ID 

Receiver 

Reference 

Oscillator 

Agency 

Austin, TX 

A 

TI-4100 

Cesium 

SDHPT 

Mojave, CA 

M 

TI-4100 

H-Maser 

CSR 

Kokee Park, HI 

K 

TI-4100 

H-Maser 

UNAVCO 

Westford, MA 

W 

TI-4100 

H-Maser 

NGS 

Yellowknife, NWT 

Y 

TI-4100 

Cesium 

CGS 


MOTE: SDHPT - State Department of Highways and Public Transportation (Texas) 

CSR - Center for Space Research 

NGS — National Geodetic Survey 

CGS - Canadian Geodetic Survey 

USGS — U. S. Geological Survey 

UNAVCO— University Navstar Consortium 


MODELS: 

Gravity , 

GM - 398600.440 km**3/s**2 
Earth radius - 6378137.0 m 
MERIT Standard ( GEM-L2 , modified) 

Luni-solar from DE-200 

Nongravitational 

Solar radiation reflectivity estimated 
y-bias estimated 

Polar Motion, UTl 

Lageos derived, reported on GE Mark III 
Center of Mass: 

Phase center in spacecraft axes with origin at center of mass 
x=0.21 m; y=0.0 m; z=0.854 m 

Measurement Model: 

Daily values of zenith delay are estimated for each site. 
Relativity correction included. 

Station Coordinates 

System defined by Lageos Laser Ranging [SSC(CSR)88 L 01]. 

VLBI coordinates [SSC(GSFC)88 R 01] rotated into SLR system 
through ties at Wettzell, Westford, Ft. "Davis and Mojave. 

Both sets of coordinates are given in the BIH Annual Report 
for 1987. The coordinates of the Ll phase center used in 
the solution are: 


Site 


x (m) 


y (m) 


z (m) 


Kokee Park: 
Mojave : 
Westford: 
Yellowknife : 


-5543817.866 

-2356214.840 

1492232.844 

-1224064.507 


-2054588.143 

-4646733.811 

-4458091.718 

-2689832.995 


2387854.349 

3668460.544 

4296045.914 

5633432.550 


The location of the Kokee Park VLBX reference point (1311) in this 
system is: 

Kokee VLB I : -5543846.014 -2054563.735 2387813.666 


The tie between the GPS antenna location and the VLBI was determined 
by Mader, as described in an NGS memorandum dated July 17, 1987. 


SOLUTION SUMMARY: 

!otes : 

All full-rate data were collected at a 30 second interval unless 
otherwise noted. 

Solution for PRN 3 was not computed due to sparcity of data. 

Austin data were not included until final coordinate solution is 
performed. 

The report for the pseudo-range orbit determination ( EW378V1 . EPH ) 
contains information on the data distribution. 

Full-rate phase data were edited and corrected for cycle slips. 

Solution was obtained with double differences only for entire week. 

File EW378D1.EPH is the ephemeris of the respective spacecraft 
center of mass. 


Number of Double Differences: 9521 


RMS of Double Differences 


0.031 m 


STATE ADJUSTMENT: 


The following table provides the adjustment in the epoch positition vector. 
The differences in position are obtained from the position predicted from 
the preceding week solution to the epoch of this report. The values in the 
table correspond to differences between the position vectors at the last 
time in the ephemerides reported for the preceding week and the first 
vectors of the current week, that is, the point of commonality between the 
current week and the preceding week. 

Epoch: 5 April 1987 00:00:00 GPS 

Week 378 


Position Difference (m) 


Satellite 

(Week 378 epoch minus 
last point of Week 377) 

PRN 

R 

T 

N 

6 

2.04 

6.43 

3.08 

9 

1.57 

5.82 

-0.29 

11 

3.97 

10.89 

-1.02 

12 

0.40 

-2.97 

-1.13 

13 

1.23 

0.75 

0.74 


R - Radial 

T - Along-track, transverse 
N - Cross-track, normal 



APPENDIX II 


Pollmeier, V. M., "A Comparison of GPS Orbits Estimated Using Observations From North 
America and Europe," Masters thesis, The University of Texas at Austin, December 1988 
(Summary and Conclusions) 



A COMPARISON OF GPS ORBITS ESTIMATED USING 
OBSERVATIONS FROM NORTH AMERICA 
AND EUROPE 

by 


Vincent M. Pollmeier, B.S. 


THESIS 

Presented to the Faculty of the Graduate School of 
The University of Texas at Austin 
in Partial Fulfillment 
of the Requirements 
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CHAPTER 8 

SUMMARY/CONCLUSIONS 


The results of the comparison of the N. American /European 
estimates with those of the N. American estimates show that there is 
some consistent and significant differences between the two sets of 
estimates in terms of the initial condition estimates, the estimated 
parameters (y-bias and solar radiation), the week long trajectories, or 
the coincident point offsets. The points of note are that the transverse 
components of most of the differences are the largest components. 
This suggest the possible existence of some type of unmodeled 
phenomena that is causing an effect in the along-track direction. This 
unmodeled force is probably non-gravitational, in nature. 

The N. American formal standard errors are 54% larger on the 
average than the corresponding errors for the N. American /European 
estimate. It is possible that this increase in the formal standard error is 
attributable to the fewer number of normal point observations used in 
the N. American estimate as opposed to the N. American/European 
estimate. However, it seems unlikely that this would result in this 
large a change and that some portion of the difference in the formal 
standard error must be due to the larger range of coverage afforded by 
the inclusion of the European observations. A future examination to 
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determine just what part of this reduction in the formal standard error 
is attributable to the better coverage could be performed by reprocessing 
the N. American/European solutions with an observation set 
somehow reduced to an equal size of N. American solution may be in 
order, but is impossible to include here. 

The large differences in orbit solutions that were possibly 
indicated by the earlier experiment (Stolz, et al., 1987), were not 
apparent. The largest differences between the two sets of estimates 
were those for Satellite 8. This is not surprising due to the nature of 
Satellite 8's clock and the behavior of other values for this satellite. 
The solutions did not consistently show any significant differences that 
would indicate a major departure in the solution. Inclusion of GPS 
observations from other locations which are even more globally 
distributed than those discussed herein, is likely to improve the orbits 
even further. 

It does appear that there is some unmodeled force acting on the 
GPS satellites. Its exact cause is not apparent, however. It is possible 
that it is a clock related problem. If this is the case, a thorough 
examination of the estimated clock parameters (not discussed in any 
depth here, due to time and length considerations) is in order. 


APPENDIX HI 


"GPS antenna performance comparison,” presented at the Third Annual Workshop on 
GPS Geodesy, Jet Propulsion Laboratory, Pasadena, Calif., March 1988 

"How much does multipath affect the TI-4100 antenna, and how much difference 
does it make?" presented at the Spring Meeting of the American Geophysical Union, 
Baltimore, Maryland, May 1989 (Abstract) 
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COMMENTARY 


The antenna experiment at Mojave was done February 6-12, 1988 (GPS week 422). The 
data for February 9 and 10 (days 40 and 41) and also for two days in the previous week, 
February 2 and 3 (days 33 and 34) of week 421, were analyzed. Table-1 shows, for all the 
satellites, the comparison between the pseudo-range noise rms of day 33 and day 40. 
Although it is not a one-to-one comparison, the rms of noise is less on day 40 than that on 
day 33, for all satellites. Table-2 shows similar comparison for days 34 and 41. 

Figures 1 and 2 show the pseudo-range noise for PRN-3 for two consecutive days in week 
421. The time origins of the plots were adjusted for the 4-minute delay in day-to-day 
repeatability of GPS, so that the figures can be laid one on top of another. When that is 
done, there is a clear indication of multipath effect. This is not the case in week 422 as 
indicated by comparison of Figures 3 and 4. Further, comparison of Figures 1 and 3 and 
Figures 2 and 4 shows the reduction in rms and the change in structure of the noise in 
pseudo-range measurements. 


T AB L- 9 - i 


dumpo 1/U4M of 1/JWulo M oiUL RMS (™) 


SV-PRN 


WKKBBM 

VJeel 455 Cd+0) 

^ObS 

RMS 

nobs 

RMS 

A 04 

R-1172 

70 5 

0-8267 

3>£8 

(-7731 

.?08 

0- 6876 

*57 5- 

^•8333 

377 

0-8735 

“14 

5- 0175 

370 

0-7582 

588 

3863 

537 

1-127/ 

35) 

5-035-5 

306 

0-7475 

541 

5-77)5 

565 

0-8764- 








TABi- £ ' 2 - 


Co-mfxxnhoA of TjmxAo Kcm^e. Moae Kk\s(W 


SV-PRN 

W&t/ic 4-2>i (T^ 34) 

W-jeeA 4-2-2 CD 4*0 

Mobs 

•RMS 

N03S 

RMS 

3 

115 

a. 5511 

534 - 

0>n2l2 

b 

54*0 

2- 5051 

4-58 

0-1037 

8 

501 

5-3517 

319 

49709 

°l 

430 

f. ma 

314 

0- 5993 

11 

635 

5-3654 

330 

0-9515 

15 

448 

1- 144A 

305 

0- 644/ 

15 

510 

5-1104 

43b, 

1 - 024? 









RNGDIF (M) 



RNGDIF (M) 




RNGDIF (M) 


(PHS R - PSD R) - AVERAGE (SV- 3 :GS-81 
WEEK NO. 422; WEEK SEC: 189209. 1 



TIME (MIN) 




RNGDIF (M) 


(PHS R - PSD R) - AVERAGE (SV- 3 :GS-81 
WEEK NO. 422; \NEEK SEC: 275699.1 



350.0 


370.0 


390.0 410.0 

TIME (MIN) 


430 . 0 


450 . 0 





ORIGINAL PAGE IS 
OF POOR QUALITY 


o ns 




V the 
VLBI 
• have 
i r ihe 
ornate 
nodel. 
m and 
ormal 
acific 
notion 
The 
i at a 
e near 
of the 
e also 


aBi 


>phy«- 


with 
mea- 
that 
>m, if 
ues, 
y the 
ed to 
\(t be 
data 
•y de- 
/LBI 
s the 
s. A 
eden 
nenta 
iting, 
-ram- 
VLBI 
it re- 
abin- 
data 
0 km 
Dr 18 
of 18 
ights 
ation 
stem 
itjr is 


skan 

ned, 

site 

(feta 

and 

arch 
The 
9 ps 


ured 
with 
LAI . 
diak 

Vyr 

che 

the 

such 

mge 

ie*r 

■?rxe 


Oopt. 
s, CA 


racy, 

have 



made GPS measurement at several triangulation 
sites with observations as early as 1857. The 
map below shows the more important 
triangulation sites. Filled triangles are 
sites where we have made high-accuracy GPS 
measurements, while open circles show sites to 
be measured in March 1989. Results will reveal 
the extent of convergence and right lateral 
shear across the Los Angeles basin in the last 
century. 
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How Much Does Multipath Affect the TI-4100 Antenna, and 
How Much Difference Does it Make? 

JL Svarc (U.S. Geological Survey, 345 Middlefield Road, MS/977, 
Menlo Park., CA 94025; Phone 415-329*4832; Biinet address 
WPRESCOTT@USGSRESV) 

PAM Abusali (University of Texas, Center for Space 
Research/WRW 402, Austin, TX 78712; 512-471-5573) 

L Young (JPL/Mail Stop 238-600, 4800 Oak Grove Dr., Pasadena, 
CA 91109; 818-354-5018) 

WH Prescott, JL Davis (also at USGS, Menlo Park) 


gun in Septrmber-Oc sober 1988 to directly monitor preaenr-day defor- 
mation throughout the eastern Mediterranean region, The September 
observations included establishing a network of 18 GPS in Tur- 
key, 4 of which were observed previously by Satellite l «i-r Ranging 
(SLR) under the W LG ENEH/M LOLAS Project. Coumuous GPS ob- 
servations were maintained at 3 SLR sites in Greece (Dionysos. As- 
kites. Rhodes), i in southern Italy (Mmera). and 3 in northern Eu- 
rope (Wcuzell.Onsala, Tromso) throughout the Turkey Cmnpugn. 

The 18 sites established in Turkey were observed during a 23 day pe- 
riod using 4 Tl-4100 dual frequency receivers provided by the Univer- 
sity Navstar Consortium (UNAVCO). The GPS site at Ankara 
(Hacettepe University) was monitored continuously throughout the 
Turkey Campaign, white most other sites in Turkey were ob served a 
minimum of 3 days each. Sites in Turkey were located to "wi iipr 
broad scale plate tectonic movements and intraplate deformations, de- 
formation* along the western and central pari* of the North Anato- 
lian fault (including strain accumulation, fault creep n ere Karabuk, 
and the distribution of strain where the fault bifurcates in northwest 
Turkey), and lithospheric extension in the Aegean Trough region of 
Southwest Anatolia. The GPS network will be extended into eastern 
Turkey during the Fall of 1989 to begin monitoring deformation asso- 
ciated with collision of the Arabian Plate with Eurasia. Selected 
baselines will be reobserved yearly, with a complete reoccupation ev- 
ery 2 years. 
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Continental extension on sets of parallel faults 

Rob_Wcstawav (Department of Geological Sciences, 
University of Durham. South Road, 

DURHAM DH1 3LE, England) 

Regions of distributed continental extension, such as Italy, 
central Greece and western Turkey, typically comprise an 
upper-crustal brittle layer containing sets erf parallel active 
faults that bound blocks that are rectangular in plan view 
and like tilted dominoes in cross-section. Evolution of such 
regions can be quantitatively predicted as a function of 

extensionai strain rate Ej and vertical vorticity x z beneath 


We have studied the effect of several Global Positioning System 
(GPS) antenna/backplane configurations on P-code multipath, and 
the effect of multipath on relative positions obtained from earner 
phase data. Two experiments were conducted at the National 
Geodetic Survey tracking station located at Mojave, California. A 
third experiment was conducted near the town of Parkfield, 
California. Data were collected using a standard TI-4100 antenna, a 
Fixed Radiation Pattern Antenna (FRPA-2), a Rogue amen- 
na/“choke-ring” assembly, and a TI-4100 antenna surrounded by mi- 
crowave-absorbing foam. In the table below, we present estimates of 
pseudorange noise determined by differencing the phase and pseudo- 
range measurements. The primary contributors to this noise are ran- 
dom pseudorange measurement error and P-code multipath. We have 
established a test network on the roof of a USGS building in Menlo 
Park, California. Even though this network appears to be in a high 
multipath environment, the relative positions of the roof stations are 
correctly determined to less than 1 mm. The difference between 
lengths determined by GPS and Geodolite (a high precision laser dis- 
tance measuring device) is 0.6 ± 0.5 mm. Although there are applica- 
tions where very accurate pseudorange measurements are needed, we 
conclude that the effect of multipath on relative positions determined 
by carrier phase is at the level of 1 mm. 


Antenna Configuration 

HaQtennfr ("cood” location - brushv area at MoiaveT 

RMS <m) 
0.52 

ILam&ina. (“bad" location - roof at USGS1 

0 87 

hEPA-2 antenna (roof at Moiave) 

0.58 

XI antenna with absorbine.foam (Parkfield) 

0.59 

HjLQlSnrra (roof at Moiave) 

0.77 

Rogue antenna/choke fines (roof at Moiave) 

QJA 


the brittle layer. The ratio x z /(2Ej), or II, is a dimension- 
less parameter, Holmes number, which describes relative 
importance of extensionai strain and rotation around vertical 

axes during evolution of such regions. H*eoc(y), where y is 
the angle between strike and common slip vector azimuth 
for faults in a parallel set. 

In central Italy, H— 0.15, and clockwise rotation during 
extension over the past -2 Myr has been negligible. H— 0.7 
in central Greece and -+0.7 in western Turkey. This inplies 
that the present active fault generation in both regions has 
taken up -10° rotation during the past -2 Myr whilst the 
region has extended by -1.2. A single previous fault 
generation took up -30° rotation during the previous -5 
Myr, whilst the region extended by -1.7. Thus, observed 
paleomagnetic rotation of -40° clockwise in Greece and -40* 
anticlockwise in Turkey during the past 7 Myr can be 
explained as a consequence of vertical vorticity causing 
oblique slip on two generations of parallel fault sets. 
Financial support is being provided by UK Natural 
Environment Research Council grant GR3/6966 


Tectonics of the Mediterranean 
and Circum-Mediterranean-I 
(G41A) 

Chesapeake III THURS AM 
R Reilinger, MIT 
K Hastens, L-DGO 
Presiding 
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1988 Global Positioning System (GPS) Crustal Deformation Measure- 
ments m Turkey 

E Reilinger O) M N Toksoz (l) A Barka (1) E Kasapoglu (2) P Wil- 
son (3) H Secger (3) J Stowcll (4) B Stephens (4) 

(l)Earth Resources Laboratory, Cambridge, MA 02139; (2) 

Dept, of Geological Engineering, Hacettepe Univ„ Ankara, Turkey; 
(3) IF AG. Frankfurt, Fed. Rep. Germany; (4) UNAVCO. Univ. of 
Colorado. Boulder. CO 

A multi-ins titutionai experiment involving scientists from the 
United States, Germany, England, Turkey. Greece, and Italy was be- 
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Kinematics of Mediterranean tec tome deformation 

Bab Westawav (Department of Geological Sciences, 
University of Durham. South Road, 

DURHAM DH1 3LE, England) 

Mediterranean active tectonics are dominated by northward 
convergence of the African plate relative to Europe. It is 
well known that this convergence is taken up in the eastern 
Mediterranean by subduction of Africa beneath Greece and 
Turkey, and in the western Mediterranean by shortening of 
the continental lithosphere in NW Africa. Recent work in 
Italy and adjacent regions of Africa has enabled the plate 
boundary position in the central Mediterranean and the sense 
of deformation along it to be deduced for the first time. 
Deformation within the African plate in Tunisia and northern 
Libya deflects its local direction of motion from northward 
towards the NE. consistent with the observed sense of 
extension in southern Italy. Relative motion between the 


This jKijjjt" may !>c freely copied. 


Vol. 70, No. 15, April 11, 1989 Eos 


Adriatic Sea, central and northern Italy and Yugoslavia can 
be regarded as anticlockwise rotation of the rigid Nonhem 
Adriatic microplate relative to Europe.This model is 
summarised in the map below, in which dots and crosses 
denote regions that are extending or shortening. 

Financial support is being provided by UK Natural 
Environment Research Council grant GR3/6966. 
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Mediterranean Geodynamics from Satellite Laser Ranging 

D. E. Smith . R. Kolcnkiewicz (Laboratory for Terrestrial Physics, 

NASA/GSFC, Greenbelt, MD 20771) 

M. IL Torrence, P. J. Dunn, R. G. Williamson, S. M. Klosko, J. W. 
Robbins (ST Systems Corp, Lanham, MD 20706) 

E. C. Pavlis (Univ. of Maryland-Astronomy Proeram, NASA/GSFC. 

Greenbelt, MD 20771) 

S. K. Fricke (RMS Technologies, Lanham MD 20706) 

The Mediterranean campaign, conducted by three transportable laser 
systems tracking the Laser Geodynamic Satellite (LAGEOS) since 
1985, has provided accurate epoch positions for sites in Italy, Greece 
and Turkey. These three-dimensional locations can be set in a 
European reference frame defined by fixed tracking instruments. The 
fixed sites in Italy, France, Austria, Germany and England have been 
tracking LAGEOS for as long as five years, enabling the 
determination of their inter-siation baseline values to an accuracy of 5 
mm/year. The most significant deviation from rigidity assumptions for 
the European plate is a compressive motion at Matera, Italy towards its 
neighbors to the north. The accuracy of baselines between the 
Mediterranean sites and Matera reaches approximately 3 cm in each of 
three dimensions for an annual solution. Similar levels of accuracy are 
obtained for relative distances between any pair of locations occupied 
by the transportable systems. 
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► Evolution of the Mediterranean Region 

M i Helman and j F Dewey (Deoartment of Earth Sciences, 
University of Oxford. Oxford, 0x1 3PR. England i 
(Sponsor: Kim hastens! 

Convergent motion of Africa, Ad'ia, and Eurasia has long 
been recognized as the cause of the Alpine orogeny 
i sensu lato), but the amount and direction of the motion 
could only be suggested by geology of the orogemc belt 
until the significance of oceanic magnetic anomalies and 
fracture cone^ was recognized. The vector (finite; 
difference between the independent seafloor spreading 
systems of the Central Atlantic t Af nca-Ngrth America 
separation! and the North Atlantic (North America- 
Eurasia separation) is the relative motion between 
Africa and Eurasia. 

The relative motion of Africa with respect to Eurasia 
has been the object of several studies since the early 
1970’s. Sharp, sudden changes in the direction of rela- 
tive motion characterize the earliest of these. Such 
changes have since been recognized as the result of the 
misinterpretation or rack of data. More recent studies 
have tended to show a smoother path, with more gradual 
changes in direction and general agreement on parts of 
the relative motion path; e.g., Late Jurassic-Earl / 
Cretaceous ESE directed, simstral strike-slip and 
Tertiary northward directed compression. Other aspects 
of Africa's motion are less widely agreed upon. These 
include the change to NE compression in the Cretaceous, 
the possibility of E-W dextral strike-slip in the mid- 
Tertiary, and NW directed motion from the Late Miocene. 

Parts of Africa’s motion path can be directly correlated 
with geologic data from the Mediterranean area, whereas 
other parts are at odds with the geologic evidence, 
part ’ :u lar ly kinematic indicators from the Alps. This 
suggests the importance of other factors Imicroplates, 
gravity spreading) in the evolution of this region. 

G41AHJ7 lUibH Mllti) 

Reactivation of Convergence Structures and Scismogenesis in the 
Trans-tensional Regime of the Apennines of Central Italy 

Leonardo Scebcr (Lamom -Doherty Geol. Obs. Palisades. NY 10964), 
and Peter A. Geiser (Dept, of Geol.and Geophy.. Univ of 
Connecticut, Sum, CT 06268) 

The Apennines of central Italy are dominated by thrusts and folds that 
accommodated a large amount of shortening in the Neogene. Although 
aspects of the depth extrapolation of some of these structures is still 
controversial, subsurface data from oil exploration and structural 
mapping strongly suggest that shallow-angle bedding -parallel thrusts 
with large displacements are the controlling structures in the upper 
crust. The current tectonic regime in the central Apennines is charac- 
terized by extension, as expressed by several large and often destruc- 
tive earthquakes with large normal faulting component and by range- 
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CASA UNO GPS ORBIT AND BASELINE EXPERIMENTS 
B. E. Schutz, C. S. Ho, P.A.M. Abusali, and B. D. Tapley 
Center for Space Research, The University of Texas at Austin 


Abstract . CASA UNO data from sites distributed in 
longitude from Australia to Europe have been used to 
determine orbits of the GPS satellites. The characteristics of 
the orbits determined from double difference phase have 
been evaluated through comparisons of two-week solutions 
with one-week solutions and by comparisons of predicted 
and estimated orbits. Evidence of unmodeled effects is 
demonstrated, particularly associated with the orbit planes 
that experience solar eclipse. The orbit accuracy has been 
assessed through the repeatability of unconstrained 
estimated baseline vectors ranging from 245 km to 5400 km. 
Both the baseline repeatability and the comparison with 
independent space geodetic methods give results at the level 
of 1-2 parts in 10®. In addition, the Mojave/Owens Valley 
(245 km) and kokee Park/Ft. Davis (5409 km) estimates 
agree with VLBI and SLR to better than 1 part in 10 8 . 

Introduction 

The purposes of the CASA UNO experiment have been 
described in detail by Kellogg et aL [1989]. Although the 
primary goal was to make epoch measurements in areas 
subjected to significant crustal motions from which strain 
information could be ultimately inferred, the experiment 
also offered an exceptional opportunity to study other 
phenomena. The establishment of a global network of 
stations to track the Global Positioning System (GPS) 
satellites in support of the experiment provided a data set 
from civilian sites that could be used for assessment of GPS 
orbit accuracy and force model fidelity. Such studies are 
important for the application of GPS to the accurate 
measurement of intersite vectors of hundreds or thousands of 
kilometers in length. 

Even though the GPS satellites orbit at high altitude and, 
as a consequence, are less influenced by errors in the 
modeling of both gravitational and nongravitational forces 
than lower altitude satellites, the orbits are still influenced by 
forces that are not completely understood. The well-known 
H y-bias" force [Fliegel et al., 1985], which is defined along 
the solar panel axis to account for observed orbit 
characteristics is possibly caused by thermal radiation or 
misaligned solar panels. Such effects are small, 
approximately 10” 9 m/s 2 in magnitude, but the complete 
time history and precise directions of the force (or forces) 
have not been established. These characteristics, among 
others, can limit the accuracy of the determined satellite 
ephemerides which, in turn, can limit their application to 
determination of intersite vectors used for geophysical 
studies. 

The purposes of this paper were to examine the fidelity of 
the models used for the analysis of CASA UNO data and to 
assess the resulting orbit accuracy. The evaluation process 
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was based on various criteria, including the comparison of 
estimated vector baselines with those determined by other 
space geodetic techniques and the repeatability of GPS- 
determined baselines. Furthermore, ephemerides resulting 
from different arc lengths have been compared and 
information about the model fidelity has been inferred. 

Data and Site Descriptions 

The CASA UNO experiment was conducted in the period 
from January 18 through February 5, 1988, corresponding to 
GPS Weeks 419 through 421. A significant portion of the 
global tracking data was provided by the Cooperative 
International GPS Network (CIGNET). Although CIGNET 
has been augmented since the CASA UNO experiment, the 
network consisted of the sites shown in Table 1 during the 
experiment. Several additional sites in Table 2 were used to 
extend the global coverage offered by CIGNET as well as to 
provide coverage to support evaluation and assessment of 
results. 


TABLE 1. CIGNET Sites: January 1988 


Site 

Instrument/Software 

Oscillator 

Austin, TX 
Mojave, CA 
Onsala, Sweden 
Richmond, FL 
Tromso, Norway 
Westford, MA 
Wettzell, FRG 
Yellowknife, NWT 

TI-4100/CORE 

TI-4100/CORE 

TI-4100/CORE 

TI-4100/CORE 

TI-4100/CORE 

TI-4100/CORE 

TI-4100/CORE 

TI-4100/CORE 

Cesium 

Hydrogen Maser 
Hydrogen Maser 
Hydrogen Maser 
Cesium 

Hydrogen Maser 
Hydrogen Maser 
Cesium 

TABLE 2. Extended Network 

Site 

Instrument/Software 

Oscillator 

American Samoa TI-4100/GESAR 

Black Birch, NZ TI-4100/GESAR 

Ft. Davis, TX TI-4100/GESAR 

Kokee Park, HI TI-4100/GESAR 

Owens Valley, CA TI-4100/GESAR 

Tidbinbilla, Australia TI-4100/GESAR 

Rubidium 

Cesium 

Hydrogen Maser 
Hydrogen Maser 
Crystal 

Hydrogen Maser 


The fiducial concept has been adopted for the study in 
which coordinates from the space geodetic techniques of 
satellite laser ranging (SLR) and very long baseline 
interferometry (VLBI) were used to establish the coordinates 
of the GPS sites. Both SLR and VLBI have been 
demonstrated to have internal consistencies of a few 
centimeters and the comparison of mutual sites have 
suggested comparable accuracy. Most CIGNET sites have 
been established in the vicinity of VLBI sites and local 
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surveys have been performed to enable determination of the 
GPS receiver coordinates in a VLBI system. The relative 
vectors used for this study between the respective VLBI 
reference points and the GPS ionospherically corrected 
(L1/L2) phase centers are given in Table 3. The specific 
markers used for GPS were RM4 at Ft. Davis, BP ARIES 3 
at Owens Valley and NM/C/194 at Tidbinbilla. Because of 
problems encountered with processing of data collected at 
Richmond and Austin, these sites were not used in the 
analysis reported in this paper. 


TABLE 3. Relative Position Vectors 


VLB I/SLR Marker 


Vector from VLB I/SLR to L1/L2 
Phase Center (m) 


* y z 


Ft. Davis 

7216 

18.751 

-38.199 

-74.106 

Kokee Park 

1311 

6.868 

-24.034 

-3.930 

Mojave 

7222 

-43.667 

21.844 

-10.209 

Onsala 

7213 

53.521 

-40.269 

-42.432 

Owens Valley 

7207 

-43.025 

79.840 

36.259 

Tidbinbilla 

1642 

-80.923 

-45.438 

81.976 

Westford 

7209 

26.296 

38.824 

30.409 

Wettzell 

7224 

6.213 

90.527 

-29.829 

Yellowknife 

7285 

59.992 

-302.568 

-122.822 


A Kokee Park became a regular CIGNET site in 

Apnl 1988, the site used for CASA UNO was a different 
location (JPL Marker 3028-S). To establish the Kokee Park 
position with respect to the VLBI shown in Table 3, an 
additional experiment was conducted during GPS Week 441. 
The regular CIGNET site coordinates were determined 
?™g Week 431 by the National Geodetic Survey 
(NGS). The Week 441 experiment was conducted by JPL 
and the USGS to determine the CASA UNO coordinates 
with respect to the CIGNET coordinates, from which the 
relative coordinates of the CASA UNO site were inferred as 
given in Table 3. 

The site at American Samoa was located near a previously 
occupied SLR geodetic marker (7096). However, the 
marker 7096 could not be occupied by the GPS instrument 
and the local survey from the GPS location to the marker has 
wen Performed. The site at Tromso was occupied in 
July 1989, with a mobile VLBI system but the results are not 
yet available. The New Zealand site was previously 
occupied with a doppler system, but uncertainty exists in 
relating the reference frames. These considerations resulted 
m the treatment of these three sites as nonfiducial sites that 
were estimated simultaneously with other parameters 
Recent results of Prescott et al. [1989] that demonstrate 

m CndCnCy ° f the 12 phase center on m °del of the 
• ^ tCn t na were used to establish the relative vectors 
n Tabl ® 3. In particular, the height of the L1/L2 phase 
center above the LI phase center was 3.9 cm for TI-4100 
benes 100 antennas and 2.3 cm for Series 4000. 

Further information concerning the local site surveys is 
given in the CSTG GPS Bulletin. Details concerning local 
vectors between the VLBI or SLR geodetic markers and 

fou " <1 NASA T ““ cal 

The phase data collected by the receivers were used in a 
double difference mode to eliminate most of the influence of 
satellite and receiver clock variations. Pseudo-range data 
were used to ensure an accuracy of one microsecond with 


respect to GPS time in the phase time tag. Preliminary 
orbits were generated from the pseudo-range data to support 
editing and correction of cycle slips. All solutions for one or 
two-weeks began at Sunday 00:00:00 of the respective week 
and the data were collected at 30-second intervals. 

Reference System and Models 

A reference system with origin at the Earth center of mass 
was used to describe the GPS satellite ephemerides. Since 
the VLBI system is insensitive to the center of mass, a 
system defined by SLR was used to establish the appropriate 
origin. Using sites that are common to both VLBI and SLR 
solutions, a seven parameter transformation was determined 
(tnmslation, rotation and scale) between selected VLBI and 
SLR solutions. The approach was similar to a technique 
used by Kolenkiewicz et al. [1985] and Murray and King 
[1988], The adopted VLBI system was produced at the 
Goddard Space Flight Center, labeled GLB484 [C. Ma and 
J. Ryan, private communication, 1989], and the SLR 
solution was obtained at The University of Texas at Austin, 
labeled CSR 8901 and based on analysis of the LAGEOS 
satellite [R. Eanes and M. Watkins, private communication, 
1989]. These solutions were generated for an epoch of 
January 1, 1988, and the coordinates were adopted for the 
CASA UNO period with no further adjustment. After 
application of the transformation parameters, the common 
VLBI and SLR sites agree below 3 cm RMS in all three 
components [J. Ray et al., in preparation, 1990], The 
resulting GPS site coordinates for the L1/L2 phase center 
were obtained by adding the relative vectors in Table 3 to 
the VLBI/SLR coordinates. The result for the GPS 
coordinates is in Table 4, identified as FSC-4. 

TABLE 4. FSC-4 GPS Site Coordinates (m) 


Site x 


Ft. Davis -1324192.413 -5332061.177 3232044.361 
Kokee Park -5543839.212 -2054587.748 2387810.010 
Mojave -2356214.813 -4646733.808 3668460.510 
Onsala 3370659.923 711877.022 5349788.065 

Owens Valley -2409643.917 -4478269.461 3838639.584 

Tidbinbilla -4461015.908 2682720.628 -3674299 714 
Westford 1492232.872 -4458091.687 4296045.902 
Wettzell 4075546.502 931825.518 4801599.280 

Yellowknife -1224064.474 -2689832.978 5633432.538 


ATAVSOb lUVAIbld 


r ^ uccil USCC1 

extensively with SLR analysis of the LAGEOS satellite 
[Tapley et al., 1985], The gravitational force model 
included the MERIT Standards [Melbourne et al 19831 
with GM - 398600.440 km 3 /s 2 and the GEM-L2 gravity 
field truncated at degree and order 8. The nongravitational 
force model used the ROCK4 solar radiation pressure model 
and included a y-bias force [Fliegel et al., 1985], Polar 
motion and UT1 were based on LAGEOS analysis [Schutz 
et al., 1988]. The Chao [1974] troposphere model was used 
and the ionosphere corrections were applied using the dual 
frequency measurements. The satellite phase center with 
respect to the center of mass was x = 0.21 m, y = 0.0 m and 
z -0.854 m. The results have been obtained with a multi- 
satellite derivative of the University of Texas Orbit 
Processor (UTOPIA) [Schutz and Tapley, 1980], Various 
parts of the data analysis have been accomplished using 
VAXstation 2000, VAX 8200, CDC Cyber 170/750 and 
Cray X-MP/24 computers. 
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Results 

As previously noted, several sites did not have antenna 
coordinates available in the adopted VLBI/SLR coordinate 
system, thereby precluding their use as fiducial sites in the 
strict sense. These sites included Black Birch, American 
Samoa and Tromso. Although the Owens Valley and the Ft. 
Davis coordinates were known in the adopted system, they 
were treated as unknown to enable their use for evaluation 
purposes. With these considerations, all sites were fixed to 
the coordinates in Table 4, except Black Birch, American 
Samoa, Tromso, Owens Valley and Ft. Davis which were 
adjusted in the estimation process without constraints. 

Regardless of the arc length used, the epoch position and 
velocity vectors were estimated at Sunday 00:00:00 GPS 
time of the appropriate week. In addition, a scale parameter 
for the ROCK4 solar radiation pressure model [Fliegel et al., 
1985] and one or more y-bias parameters were estimated for 
each satellite. A zenith troposphere delay parameter was 
estimated for each three hour interval at each station during 
the observing periods. In all cases, the orbit was determined 
from double difference phase measurements, with pseudo- 
range measurements used to correct phase time tags to GPS 
time. A bias parameter was estimated also for each 
respective combination of two satellites and two ground 
stations in the double difference mode. No a priori 
constraints were applied to any estimated parameter, i.e., all 
parameters were completely free to simultaneously adjust 

One-week solutions, or one-week arcs, were computed for 
each respective week of the experiment. In addition, two- 
week arcs were computed using Weeks 419 and 420 as well 
as 420 and 421. The double difference (DD) residual root- 
mean-square (RMS) and the number of measurements in 
each week are given in Table 5. The statistics include cases 
resulting from the modeling of multiple y-bias parameters 
for the two-week arcs. The one week arcs used a single y- 
bias for each satellite. The higher RMS in the two-week 
arcs compared to the one-week arcs is indicative of 
remaining unmodeled effects. The experiments with 
increased number of y-bias parameters demonstrate that the 
unmodeled force cannot be satisfactorily explained by a 
force along the solar panel axis, the spacecraft-fixed y- 
direction, since only modest reduction in the RMS was 
achieved. Although the DD RMS level for the two-week 


TABLE 5. Solution Statistics 


Number of Double Differenced DD Residual 
Week No. Phase Measurements RMS (m) 

419 

82,039 

0.0297 

420 

91,115 

0.0327 

421 

78,632 

0.0277 

Two- week cases, one y-bias for the arc*: 


419/420 

141,812 

0.09702 

420/421 

134,785 

0.05645 

Two-week cases, two y-bias for the arc*: 


419/420 

141,812 

0.05567 

420/421 

134,785 

0.05252 

Two- week cases, four y-bias for the arc*: 


419/420 

141,812 

0.05208 

420/421 

134,785 

0.05047 


♦The number of y-bias parameters applies to each 
satellite. The two-week arcs did not include Owens 
Valley and Tromso. 


arcs is at the 5 cm level, experience has shown that the error 
distribution represented by this magnitude can adversely 
affect baseline vector estimates. 

The characteristics of the solutions described in Table 5 
were further investigated by comparing the two-week 
solutions with the individual one-week orbit solutions. The 
two-week solutions in which two y-bias parameters per 
satellite were estimated were chosen for the comparison. As 
expected, the differences exhibit a behavior characterized by 
a twelve-hour period (once per revolution). The two-week 
orbit determined for Weeks 419 and 420 was compared to 
the one-week solution for Week 419. The differences for all 
satellites were confined to the range of ± 10 m, however, 
satellites with the PRN identity of 3, 8, 11 and 13 also 
exhibited either a long period (>14 days) or quadratic trend 
in the along-track direction in addition to the once per 
revolution periodicity. It is presumed that the one-week 
orbit is a more accurate representation of the actual orbit, 
hence, the differences are one illustration of problems in 
force model fidelity on the two-week arc. 

Another approach to the assessment of the orbit accuracy 
and model fidelity is to use the estimated parameters from 
one-week for the prediction of the orbit into the next week. 
The predicted orbit can be compared with the orbit estimate 
for the second week. Errors in either the state estimates or 
the force model will introduce errors in the predicted states. 
The orbits determined for Week 419 were predicted into 
Week 420 and compared with the determined orbit for Week 
420. While PRN 6, 9 and 12 exhibit differences bounded by 
± 5 meters with a periodicity of once per revolution, the 
other satellites exhibit a quadratic-like growth in the along- 
track component reaching 50-100 m after one week. The 
satellites share the common characteristic that PRN 3, 8, 11 
and 13 experience eclipsing of the Sun by the Earth and the 
other satellites do not. As a consequence, it is tentatively 
concluded that an effect associated with the eclipse season 
remains unmodeled. 

The one-week orbits described in the preceding 
paragraphs were evaluated further with applications to the 
estimation of baseline vectors, the essential test of geodetic 
accuracy. In some cases, the assessment was based on the 
repeatability of the baseline estimate from different 
experiments. While this is a necessary test for accuracy, it is 
not sufficient Still other cases allow a comparison with 
values determined by independent techniques such as SLR 
or VLBI which are subjected to different error sources. 

The site solutions obtained for each individual week are 
given in Table 6 as the vector relative to the fixed site used 
for forming the double differences. The solution is given as 
the average of the three weeks and the repeatability from 
week to week is indicated by the RMS. The results indicate 
length repeatability at the level of one or two parts in 10 8 , 
except for Tidbinbilla to American Samoa which is about 
four parts in 10 8 . It should be noted that if the Week 419 
estimate for American Samoa is eliminated, the remaining 
two solutions agree at the level of one or two parts in 10 8 . 

As noted previously, only the Kokee Park/Ft. Davis and 
Mojave/Owens Valley have been independently determined 
by another technique, namely, VLBI and SLR. Table 7 
shows the comparison of the estimated baseline (Table 6) 
with VLBI/SLR vectors determined from Table 4, resolved 
into longitude, latitude and height differences. Both cases 
show agreement with the independent determination of one 
part in 10 8 or better in the baseline magnitude. These results 
have been determined without a priori constraints in 
solutions that were allowed to freely adjust to a converged 
solution. The quality of the 245 km baseline results is 
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TABLE 6. Estimates Averaged from One Week Solutions (m) 


Baseline 

Ar 

Ay 

A 2 

Length 

Tito AS 

RMS 

-1638933.040 

0.044 

-3680013.690 

0.144 

2105673.934 

0.101 

4545598.419 

0.177 

TItoBB 

RMS 

-278055.991 

0.005 

-2168360.092 

0.072 

-551140.618 

0.020 

2254519.151 

0.066 

ON toTR 

RMS 

-1267718.477 

0.023 

9692.532 

0.029 

608405.800 

0.026 

1406186.936 

0.011 

KO to FD 

RMS 

4219646.772 

0.024 

-3277473.446 

0.047 

844234.346 

0.022 

5409249.735 

0.048 

MO to OV 
RMS 

-53429.130 

0.012 

168464.348 

0.033 

170179.065 

0.026 

245348.370 

0.007 

AS-American Samoa 
KO-Kokee Park 
OV-Owcns Valley 

BB-Black Birch 

MO-Mojave 

TI-Tidbinbilla 

FD-Ft. Davis 

ON-Onsala 

TR-Tromso 




TABLE 7. Difference (m) of Table 6 with Table 4 


Baseline 

Longitude 

Latitude 

Height 

Length 

MO to OV 
KOtoFD 

-0.014 

-0.013 

-0.007 

-0.008 

0.004 

0.017 

0.000 

0.012 


consistent with results of other investigators [Blewitt, 1989; 
Dong and Bock, 1989]. 

Conclusions 

Evidence has been presented for the existence of 
unmodeled GPS force characteristics that may, in part, be 
associated with eclipse season. While understanding these 
forces is important, results have been obtained for vector 
baselines that demonstrate 1-2 parts in 10 8 repeatability and 
comparison with SLR and VLBI determinations. 

Compensation of the error sources by stochastic techniques 
or improved understanding of the nature of the forces could 
contribute to further improvements in the baseline results. 
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ABSTRACT 

Experiments have been conducted with GPS data collected during the CASA UNO campaign 
(January, 1988). This campaign collected data from an enhanced CIGNET by including data from 
Australia, New Zealand and American Samoa. The experiments have focused on the model fidelity 
associated with one and two week arcs computed during the campaign period. The orbit accuracy 
was evaluated, in part, through comparison of the satellite positions at the time point common to 
adjacent one week arcs and by examining the week to week consistency of station 
coordinate/baseline estimates. Based on the consistency tests and the comparison of baseline 
estimates with those determined by other techniques, the accuracy has been estimated at a few parts 
in 10 8 . 

INTRODUCTION 

Since the last Geodetic Symposium held in Austin, TX, in 1986, significant progress has been made 
x * n the organization of a civilian network of GPS receivers, primarily dedicated to regular tracking of 
the GPS Block I constellation. In September, 1986, with the sponsorship of the U.S. Geological 
Survey (USGS), The University of Texas at Austin Center for Space Research (CSR) organized a 
cooperative network of three TI-4100 receivers. This network consisted of the following locations 
and cooperating agencies: 

• Mojave, CA: receiver owned by CSR, operated by National Geodetic Survey (NGS) using a 
hydrogen maser frequency reference; 

• Austin, TX: receiver owned and operated by the Texas State Department of Highways and 
Public Transportation using a cesium reference; 

• Westford, MA: receiver owned and operated by the NGS using a hydrogen maser reference. 

This network was used to support regional geodetic applications and orbit determination. 

Through the efforts of the NGS, an expanded and more global network has been organized. This 
network of cooperating stations, the Cooperative International GPS Network (CIGNET) now 
consists of the sites given in Table 1 (CSTG Bulletin, 1989). 

Early orbit determination studies were limited by the regional nature of the Continental U.S. 
network. In the presence of no model errors and precise measurements, accurate orbits can be 
determined using such a regional network. However, in the presence of complex GPS clock 
variation and uncertainties in force models, a regional network is inadequate for orbit determination. 
Some difficulties can be diminished by use of a global tracking network which allows observation of 
the satellite over a greater portion of the orbit, thereby providing enhanced information about the 
model error characteristics and enabling improved averaging of model errors. Nevertheless, model 
errors continue to exist and are also dependent on the temporal data span, or arc length. It is 
generally acknowledged that as the arc length is lengthened, the effects of model errors will become 
more pronounced. On the other hand, the pronounced effect can be used to study the phenomena 
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and, hence, improve the modeling and understanding of it. 

The purpose of this paper is to examine the fidelity of the models through the use of GPS tracking 
data in different arc lengths. The evaluation process uses various criteria, including the arc to arc 
consistency for selected baseline estimations and the differences in orbits for two adjacent arcs. This 
paper is not intended to be an exhaustive treatment, but instead to document interim results obtained 
from a selected set of cases in an ongoing study. 

DATA 

For the study described in this paper, the CASA UNO Campaign (Kellogg, et al., 1989) was chosen 
to augment the regular TI— 4100 data set of the CIGNET. This campaign was conducted for a period 
of almost three weeks and used sites which enhanced the global distribution. The additional sites, 
which used TI-4100 receivers, are given in Table 2. While the purpose of this campaign was to 
obtain a set of epoch measurements in Central and South America, additional global tracking 
contributions were organized to provide improved orbit determination support. It should be noted 
that the CIGNET site at Kauai (Kokee Park) was not in operation during CASA UNO, however, a 
TI-4100 receiver was operated in the vicinity of the permanent CIGNET site which was installed 
about three months after the campaign. 

The results described in this paper used the first two weeks of the CASA UNO, namely, GPS Week 
419 (January 17-23, 1988) and Week 420 (January 24-30, 1988). All stations used GESAR or a 
closely related product (CORE) developed by the Applied Research Laboratory of the University of 
Texas at Austin. The CIGNET sites used either hydrogen maser or cesium frequency standards and 
collected data at 30 sec intervals. 

The pseudo-range data were used from all sites to determine preliminary orbits for phase editing and 
determination of receiver time tag corrections. Phase processing required several steps, including 
the tasks of locating and fixing cycle slips. The phase data were used in a double difference mode 
for orbit and parameter estimation. 

The arcs considered in this paper were one or two weeks in duration. All arcs began at Sunday 
00:00:00 GPS Time of the respective week. 

REFERENCE SYSTEM 

Most of the sites listed in Table 1 are near Very Long Baseline Interferometry (VLBI) or Satellite 
Laser Ranging (SLR) sites whose coordinates have been determined to an accuracy of a few 
centimeters. Local surveys, performed with classical and GPS techniques, provide the tie between 
the GPS antenna and the VLBI or SLR reference point. Since VLBI is insensitive to the Earth center 
of mass required for the description of satellite orbits, a combination of the VLBI/SLR reference 
frames was used for the GPS reference frame. 

An analysis similar to a technique described by Kolenkiewicz etal., (1985) and Murray and King 
(1988) was applied to obtain the GPS reference frame. In the approach used in this paper, the SLR 
coordinates SSC (CSR) 88L01 (Schutz et al., 1988) were mapped to an adopted epoch of 
January 1, 1988, using the tectonic motion model of Minster and Jordan (1978). Similarly, the 
Earth-fixed VLBI coordinates of Goddard Space Flight Center, GLB206 (Ryan et al., 1988) were 
mapped to the same epoch. Although the GLB206 was replaced by the solution GLB223, the 
differences between the two solutions at the sites of interest is small. 
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Several mutual or nearby collocations of SLR and VLBI exist. The sites used for this study and the 
corresponding geodetic markers are given in Table 3. The relative locations of appropriate geodetic 
markers were obtained from NASA Tech. Ref. 1198 (1988). By adding the marker relative position 
vector to the VLBI coordinates, a set of coordinates for the SLR sites expressed in the VLBI 
reference frame was formed. 

The SLR coordinates expressed in the VLBI system were compared with those derived from SLR 
analyses and a seven parameter transformation set was obtained. This transformation was applied to 
other VLBI coordinates in order to map them into the SLR reference frame. The resulting 
differences between the coordinates of the VLBI sites transformed into the SLR system compared to 
the SLR coordinates are given in Table 4. 

Analysis documented by Chin (private communication, 1989) for the local ties was used to obtain 
the GPS antenna coordinates after the use of the transformation parameters. The resulting GPS 
coordinates for the fixed (unadjusted) CASA UNO LI -phase centers, identified as CSR Fiducial 
Station Coordinates, FSC-2, are given in Table 5. The coordinates for the special CASA UNO fixed 
sites used operator measurements of the antenna height above the respective geodetic marker. 
Additional information on local surveys at the Tidbinbilla site was provided by the Australian Lands 
and Information Group (Steed, personal communication, 1988). The Tidbinbilla site was the marker 
identified as NM/C/194. The Ft Davis and Owens Valley sites refer to markers RM-4 and BP 
Aries-3, respectively, and the coordinates of these markers with respect to VLBI sites can be found 
in NASA Tech. Ref. 1198 (1988). 

MODELS 

The models used for the GPS analysis were based primarily on experience derived from Lageos laser 
ranging investigations (Tapley et al., 1985) and the adoption of the SLR geocentric coordinate 
system described in the preceding section. The adopted force model for GPS analysis is given in 
Table 6. To a large extent, the models have been used and validated on several satellites, including 
SLR targets and radiometric satellites. 

The orbit analysis software, MSODP, has been developed at the University of Texas at Austin, partly 
derived from UTOPIA (Schutz and Tapley, 1980). Various parts of the computations have been 
done on a VAX 8200, CDC Cyber 170/750 and a Cray X-MP/24. 

RESULTS 

For the two weeks examined, most data collected with the enhanced CIGNET was of good quality, 
although some sites experienced excessive cycle slips at some times. Only one site (Richmond) 
experienced significant problems that resulted in little usable data. 

During the CASA UNO, several sites did not have antenna coordinates available in the adopted 
VLBI/SLR coordinate system. These sites were Black Birch, American Samoa, Kauai and Tromso. 
Although Owens Valley coordinates were known in the adopted system, the site coordinates were 
treated as unknown to enable use of the site for evaluation purposes. During Week 446, a local GPS 
survey was conducted that enabled determination of the CASA UNO Kauai site in the adopted 
system. A VLBI occupation of the Tromso site is planned for summer, 1989, thus enabling 
assessment of the coordinates given later in this paper. 

Regardless of the arc length, the epoch position and velocity vectors of each satellite were estimated. 
In addition, a scale parameter on the ROCK4 model and one or more y-bias parameters were 
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estimated also. A zenith delay parameter was estimated for each three hour interval during 
observing periods. In all cases, double differenced phase measurements were used to estimate the 
parameters of interest. These measurements require estimation of a bias parameter for each 
respective combination of two satellites and two ground stations. 

The solution statistics obtained for Weeks 419 and 420 are given in Table 7. As readily seen from 
this table, over 65,000 double differenced measurements were used in the one week arcs. The 
respective one week residual RMS is a little more than 3 cm. In the case where the arc length is 
extended to two weeks, the residual RMS increases by almost a factor of two, although the signal 
level expected to cause such an increase is a little more than the RMS, namely 4 cm. The cause of 
the increase is under investigation, but the nongravitational forces are suspected of being major 
factors. 

The orbit position differences at the common time point between the arcs of Weeks 419 and 420 are 
given in Table 8. While the differences at common time are not a definitive assessment of the orbit 
accuracy, they are an indicator of the level at which unmodeled effects still exist Even though the 
two week arc will remove this common time difference, the comparison at the initial and final 
epochs demonstrates several meter differences. 

Another means of assessing the orbit quality is by examining the week to week consistency of the 
baseline estimates and the comparison of baselines with VLBI/SLR determined coordinates, where 
available. Table 9 illustrates die week to week comparisons for the coordinates of the sites estimated 
in the one week arc experiments. In this table, the second site was estimated and the first held fixed 
to the coordinates given in Table 5. With the exception of American Samoa, the week to week 
consistency is at the level of a few centimeters on baselines as long as 5409 km, or a few parts in 
10 . The behavior of the American Samoa estimate requires further investigation. 

In some cases the ability to estimate accurate three-dimensional coordinates using a long baseline is 
of interest. With the local survey performed at Kauai during Week 441 by JPL and USGS, the 
coordinates in the FSC-2 system have been determined and can be compared with the estimates 
obtained with the analysis of Weeks 419 and 420. The results are given in Table 10. While 
differences between the estimates and the survey are about 10 cm, the differences are consistent with 
a few parts in 10 8 . Nevertheless, the discrepancy is under further investigation. 

With the VLBI determination of Tromso coordinates planned for summer, 1989, the GPS determined 
coordinates of Tromso can be assessed. The results from Weeks 419 and 420 are given in Table 1 1 
relative to Onsala. These coordinates refer to the GPS receiver antenna phase center which will have 
to be corrected to an appropriate reference mark for comparison with VLBI. 

CONCLUSIONS 

Comparison of the Week 419 and 420 orbit solutions at a common time point shows differences at a 
few meter level, even for PRN 8 which is difficult to handle because of the clock modeling required 
with pseudo-range. It has also been shown that the baseline week to week consistency is a few parts 
in 10 8 for baselines up to 5400 km. The two- week solution, as well as the common point 
differences, suggest the existence of remaining unmodeled effects, possibly caused by 
nongravitational forces. Comparison of the three-dimensional coordinates of Kauai with a later GPS 
local survey shows agreement at a few parts in 10 8 from a 5000 km baseline. Tromso coordinate 
estimates are given for future comparison with a VLBI determination. 
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Future analyses will focus on the third week of CASA UNO as well as the fidelity of the force model 
required on two and three week arcs. It is believed that the results of the long arc analyses will 
contribute to the short arc results through improved understanding of the forces and development of 
improved models. 
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TABLE 1. C1GNETT1-4100 SITES (JANUARY, 1989) 


Kokee Park, Kauai, HI Tiomso, Norway 
Mojave, CA Westford, MA 

Onsala, Sweden Wett 2 ell, FRG 

Richmond, FL Yellowknife, NWT 


TABLE 2. ADDITIONAL CASA UNO SITES (TI-4100) 


Black Birch, New Zealand 
Canberra (Tidbinbilla), Australia 
Ft. Davis, TX 
Owens Valley, CA 
Pago Pago, American Samoa 


TABLE 3. VLBI/SLR TRANSFORMATION SITES 

Location 

VLBI 

Mark 

SLR 

Mark 


Ft. Davis, TX 
Mojave, CA 
Westford, MA 
Wettzell, FRG 


TABLE 4. TRANSFORMATION SITE DIFFERENCES AFTER APPLICATION 

OF 7-PARAMETER TRANSFORMATION 

SLR Differences (m) 

Location Mark Ax I Ay I Az 


Ft. Davis, TX 
Mojave, CA 
Westford, MA 
Wettzell, FRG 


7834 

0.043 

RMS 

0.051 
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TABLE 5. FSC-2 COORDINATES (m) 

Site 

X 

y 

z 

Mojave, CA 

-2356214.8400 

-4646733.8110 

3668460.5440 

Westford, MA 

1492232.8440 

-4458091.7180 

4296045.9140 

Yellowknife, NWT 

-1224064.5070 

-2689832.9950 

-5633432.5500 

Wettzell, FRG 

4075546.3550 

931825.4590 

4801599.1700 

Onsala, Sweden 

3370659.8770 

711876.9860 

^349788.1050 

Ft. Davis, TX 

-1324192.4360 

-5332061.1350 

3232044.3480 

Owens Valley, CA 

-2409643.9260 

-4478269.4230 

3838639.5860 

Canberra, Australia 

-4461015.8080 

2682720.5220 

-3674299.5410 

| Epoch: January 1, 1988 Coordinates refer to phase center 


TABLE 6. GPS FORCE KINEMATIC AND MEASUREMENT MODELS 

Gravitational Force: 

Earth 

GM- 398600.440 hn^s 2 
R- 6378137.0 m 

GEM-L2, MERIT Standards (Melbourne, et al.* 1983) truncated at degree and order 8, with 
the value of J2 corrected for the permanent tide 
Other Bodies 

Sun and Moon from DE-200 

Nongravitational Force: 

Radiation Pressure 

Direct solar, ROCK 4 model, scale factor estimated 
Other 

y-bias, estimated parameter in spacecraft y-axis 

Kinematic Models: 

Precession, Nutation 

IAU (see MERIT Standards) 

Earth Rotation 

UT1 tide variations (see MERIT Standards) 

Polar motion and UT1, determined by LAGEOS SLR (reported weekly on GE Mark 3) 

Measurement Models: ' - 

Relativity 

Chao troposphere model, zenith delay estimated for each day 
Ionosphere correction applied from dual frequency 
Phase center in spacecraft axes with respect to center of mass: 
x m 0.21 m; y = 0.0 m; z - 0.854 m 
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TABLE 7. SOLUTION STATISTICS 

Week 

Number of Double 
Differenced Measurements 

Residual 

RMS(m) 

419 

h “ 66,735 

0.03127 

420 

75,077 

0.03447 

419/420 

141,812 

0.05223 


TABLE 8, POSITION DIFFERENCES AT COMMON TIME 

Satellite PRN 

Component 1 f~6 8 f~9 TTi TT5 T 


x 2.58 

y 0.36 

z 0.08 


9 

11 

12 

13 

2.48 

■ 

.25 

2.01 

1.96 

-1.22 

-1 

.89 

0.72 

-4.18 

-2.36 

0 

.56 

0.87 

0.87 


TABLE 9. WEEKS 419 AND 420 COMPARISONS 



Week 420-4 19(m) 


Baseline 

Ax 

Ay 

Az 

Ab 

Canberra/ A. Samoa 
b- 4545598m 

-0.147 

-0.212 

40.087 

0.265 

Canberra/Black Birch 
b = 2254519m 

-0.012 

-0.070 

40.045 

0.048 

Ft. Davis/Kauai 
b = 5409249m 

-0.006 

+0.085 

-0.038 

0.062 

Mojave/OVRO 
b = 245348m 

-0.040 

-0.030 

0.034 

, 0.011 

Onsala/Tromso 
b = 1406186m 

-0.021 

-0.027 

-0.023 

-0.009 


TABLE 10. WEEKS 419 AND 420 KAUAI ESTIMATION 


Week 419 -SM38.3y.250 

Week 420 -5543839.256 

Week 441 Survey | -5543839.248 
Ft. Davis/Kauai Baseline: 


Week 419 
Week 420 
Week 441 Survey 


-5543839.250 -2054587.639 2387809.937 
-5543839.256 -2054587.554 2387809.899 
-5543839.248 -2054587.750 2387810.051 


' 5409249.809m 
5409249.871m 
5409249.722m 
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Precise geodetic measurement can provide an excellent study of plate tectonic motion. 
The potential of Global Positioning System (GPS) has been shown in determining the 
precise baseline. The precise GPS satellite orbit is one key to high-accuracy long baseline 
(>200 km) solutions. This study examines GPS orbit determination using two L-band 
signals transmitted by GPS satellites. 

The research concentrates on the analysis of the pseudo-range and phase data 
collected during the CASA UNO campaign. The fidelity of the applied force and derived 
observation models has been tested by processing the tracking data from the campaign with 
the unconstrained (i.e., no a priori state covariance) batch estimation approach. The orbit 
accuracy has been assessed through the baseline solution and the orbit comparison. 

Based on the results from the tracking data experiments, the precise GPS orbit 
resulted in baseline accuracy of a few parts in 10 8 can be obtained by using a weekly data 
set from a globally distributed tracking network. In addition, the potential that the precise 
baseline can be independently determined has been demonstrated if a global orbit is 
available. Evidence has been shown for the existence of unmodeled GPS force 
characteristics that is, in part, associated with Earth shadowing. The approach with several 
subarcs to model the acceleration in satellite body-fixed Y-direction improves the multi- 
week orbit quality, especially of eclipsing orbits. 
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ABSTRACT j 

In July 1988, four TI-4100 receivers occupied sites in the Southwest Pacific to collect epoch 
measurements for tectonic studies. Three UNAVCO receivers operated for about ten days at Upolu 
(Western Samoa), Rarotonga (Cook Islands) and Vava’u (Tonga). A TI ROM receiver on loan to 
UNAVCO was operated at Tongatapu (Tonga). Other sites were occupied for a period of two days. 
Fiducial station support was provided with the operation of a receiver at OrToral, Australia, and a 
Mini-Mac receiver operated at Tsukuba, Japan. The CIGNET sites provided data from locations in 
the United States, Canada, Federal Republic of Germany, Sweden and Norway. The Southwest 
Pacific baselines vary in length; the longest baseline is about 1600 km. This paper describes the 
overall quality of the first six days of the 1988 Campaign data and presents the results of geodetic 
analyses of these data. Various cases have been considered and the internal consistency of the 
baselines was found to be a pan in Hr on 1600 km baselines. 

INTRODUCTION 

In July 1988, a GPS campaign was conducted in the Southwest Pacific to collect epoch 
measurements to be used for crustal motion studies in the plate boundary complex of the Tonga 
region. In this region where the Indo-Australian and Pacific plates create a complex crustal 
deformation, the relative plate velocities are expected to be in the vicinity of 18 cm/yr. 

The epoch measurements of July 1988 were collected at sites which are a subset of sites planned for 
occupation in the Main Campaigns of 1989, 1990 and 1992. The Main Campaigns will extend from 
Rarotonga on the Pacific Plate to New Caledonia on the Indo-Australian Plate. The purposes of the 
1988 Campaign were: 

• Collect epoch measurements at a few selected sites 

• Gain experience in Southwest Pacific logistics, GPS data collection and data analysis 

• Use experience to plan the Main Campaign. 

The 1988 Campaign was conducted during July 11-21 (Days 193-203) and the participating 
institutions are given in Table 1. The sites visited in 1988 are illustrated in Fig. 1. Further 
information on details of the July 1988 Campaign are given by Bevis (1989) who serves as the 
Project Coordinator. 


The purpose of this paper is to summarize analysis and results of data collected during the first six 
days of the 1988 Campaign corresponding to GPS Week 444. The analysis includes an examination 
of the sensitivity of the results to various parameters. 

DATA 

Data collected by the CIGNET (CSTG Bulletin, 1989) played an important role in the 1988 
Campaign for the determination of GPS orbits. These data were augmented by equally important 
contributions from Orroral, Australia, and Tsukuba, Japan. Further informadon concerning the 
fiducial sites and the field sites are given in Table 2 and die global distribution of sites is illustrated 
in Fig. 2. In the studies performed, Tromso was not treated as a known, fixed site but was allowed to 
adjust in the same manner as the field sites. 

In general, most receivers performed excellently during die the campaign. One receiver (Richmond, 
FL) experienced problems and contributed little data. Occasional problems were experienced at 
some CIGNET sites. Because of a software problem, the TI ROM receiver located at Tongatapu 

recorded data for some periods that was out of time synchronization with the other receivers by nine 
seconds. 

The pseudo-range data from all sites were used to provide corrections to the phase time tag. Phase 
data were used in a double difference mode from the current operating Block-I satellites, namely, 
PRN 3, 6, 8, 9, 1 1, 12 and 13. 

MODELS AND REFERENCE FRAME 

The force and kinematic models used in the analysis of the July 1988 Campaign are described by 
Schutz, et al. (1989). The technique used to obtain the reference frame and the resulting fiducial 
coordinates can be found in the same source. The adopted coordinates for the Southwest Pacific 
campaign, identified as FSC-2, are given in Table 3. Information provided by AUSLIG (Steed, 
private communication, 1988) and GSI (M. Tsuji, private communication, 1988) enabled 
establishment of the Orroral and Tsukuba coordinates in the FSC-2 system. 

The set of estimated parameters associated with the dynamical description of the satellite motion are 
given in Table 4. In addition, coordinates of all experiment sites in the Southwest Pacific were 
estimated in addition to the coordinates of Tromso, Norway. A priori constraints were not used with 
any estimated parameters, i.e., the a priori covariance on estimated parameters was infinite. 

RESULTS 

A single orbital arc spanning the entire Week 444 was used to determine the GPS orbits for analysis 
of the field sites. Experiments were performed to investigate the sensitivity of the results to force 
model characteristics and data contributions. Consideration of data contributions included the 
influence of the Tongatapu receiver synchronization which was investigated by treating the data with 
less than one second synchronization as a separate and independent "site" from the data collected 
with a nine second synchronization offset. The cases considered are summarized in Table 5. 

The three-dimensional coordinates for the field sites were estimated, but the results are given in 
terms of the intersite distance (baseline) from the "base" station, Rarotonga, in Table 6. Rarotonga 
was chosen as the base station because of its distant location on the Pacific plate. Although the orbit 
changes from case to case by a few centimeters to several meters, the 1500 km baselines show 
changes at the one to two centimeter level or better than a part in 10 . 



To demonstrate the influence of the fiducial coordinates, the fixed coordinates of Kauai were 
changed by about 10 cm and Case A was repeated. This change in the coordinates produced a 
change in the absolute length of the baseline of two centimeters. The other cases produced 
comparable shifts in the absolute value and the same level of internal consistency from case to case 
as shown in Table 6. 

The estimated three-dimensional coordinates of Tromso arc shown in Table 7. For comparison, the 
estimate is compared with a result obtained for GPS Week 419, as given by Schutz, et al. (1989). 
The comparison between Weeks 419 and 444 show agreement at the few centimeter level in all three 
components for a 1400 km baseline. As a further comparison, the Case A orbit, which excluded 
Tromso in the orbit determination, was held fixed and the Tromso coordinates were estimated. The 
result is in less agreement with the cases in which Tromso simultaneously contributed to the orbit 
determination and the estimation of the site coordinates, although the baseline length is consistent. 

CONCLUSIONS 

Based on the analysis of available CIGNET and field data collected during GPS Week 444, the 
results presented in this paper show internal consistency on the 1500 km intersite baselines at the 
level of 1-2 cm, corresponding to a part in 10 . No external standard is currently available in this 
region to serve as a comparison. Although in a different region of the world, a determination of 
Tromso, Norway, was performed that produced a result in good agreement with determinations from 
CASA UNO data. A VLBI visit scheduled for summer, 1989, will enable an evaluation of the 
Tromso results. 

Lack of receiver synchronization at the 10 second level produced a small difference in baseline 
length (3 cm) when compared to results obtained with less than one second time synchronization. 
Other factors may have contributed to the 3 cm difference, such as the influence of data loss from 
fiducial stations, especially Kauai, experienced during some intervals. 

Future analysis will focus on the remaining data in Week 445, which comprises about 40% of the 
total data set. Some differences exist in the global tracking in the remaining data sets, such as no 
available data from Onsala. Further refinements will be completed in the processing of Week 444, 
such as the use of fiducial coordinates mapped to July, 1988, rather than the coordinates for 
January 1, 1988, used in the analysis described in this paper. It is expected that the absolute values 
of the baselines will slightly change with the modified fiducial coordinates, but the internal 
consistency is expected to be unchanged. 
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TABLE 1 


INSTITUTIONS AND AGENCIES PARTICIPATING 
IN THE SOUTHWEST PACIFIC GPS PROJECT 


North Carolina State University (USA) — — - 

University of Colorado (USA) 

University of Texas at Austin (USA) 

ORSTOM (France, New Caledonia 
National Geodetic Survey (USA) 

Cornell University (USA) 

Jet Propulsion Laboratory (USA) 

Institut Gtfographique National (France) 

University of New South Wales (Australia) 

Australian Surveying and Land Information Group (Australia) 
Department of Survey and Land Information (New Zealand) 
Department of Lands and Survey (Western Samoa) 

Ministry of Lands, Survey and Natural Resources (Tonga) 
Lands and Survey Branch (Niue) 

Mineral Resources Department (Fiji) 

Topographic Service (Vanuatu) 

Service Topographique (New Caledonia) 

Geographical Survey Institute (Japan) 


TABLE 2. DATA 

Fiducial Sites 

Orroral, Australia 

TI4100/GESAR 

Cesium 

Tsukuba, Japan 

Mini-Mac 

Cesium 

Kauai, HI 

TI4100/CORE 

H-Maser 

Mojave, CA 

TI4100/CORE 

H-Maser 

Yellowknife, NWT 

TI4100/CORE 

Cesium 

Westford, MA 

TI4100/CORE 

H-Maser 

Richmond, FL 

TI4100/CORE 

H-Maser 

WettzeU, FRG 

TI4100/CORE 

H-Maser 

Onsala, Sweden 

TI4100/CORE 

H-Maser 

Tromso, Norway 

TI4100/CORE 

Cesium 

Field Sites 



i 


Rarotonga (Cook Islands) 
Upolu (Western Samoa) 
Vava’u (Tonga) 
Tongatapu (Tonga) 


Rubidium used on occasion 


TI4100/GESAR Crystal 
TI4100/GESAR Crystal 
TI4100/GESAR Crystal 
TI4100/RQM C 


at field sites. 






TABLE 3. FSC-2 COORDINATES (m) 


Mojave, CA 

-2356214.8400 

-4646733.8110 

3668460.5440 

Westford, MA 

1492232.8440 

-4458091.7180 

4296045.9140 

Yellowknife, NWT 

-1224064.5070 

-2689832.9950 

-5633432.5500 

Wettzell, FRG 

4075546.3550 

931825.4590 

4801599.1700 

Onsala, Sweden 

3370659.8770 

711876.9860 

5349788.1050 

Kauai, HI 

-5543818.2778 

-2054582.7942 

2387858.4651 

Orroral, Australia 

-4446476.5730 

2678105.1735 

-3696262.2145 

Tsukuba, Japan 

-3957193.3880 

3310191.6530 

3737733.3390 


TABLE 4. ESTIMATED PARAMETERS 


Arc Length: 7 days 

Epoch position, velocity of each satellite 

Scale parameter for ROCK4 radiation pressure for each satellite 

y -bias, one for each satellite 

Troposphere zenith delay (3 hour) for each station 

Station coordinates (selected) 

Double difference biases 


TABLE 5. STUDY CASES 


Description 


A Excluded data: Tsukuba, Tongatapu, Tromso 

B Same as Case A, except three y —bias parameters estimated instead of one 
C Excluded data: Tsukuba, Tromso 

D Excluded data: Tromso, Tongatapu 

E Excluded data: Tsukuba, Tongatapu 


Estimated parameters are defined in Table 4, except as noted. 










TABLE 6. BASELINE RESULTS* 



Distance (m): Rarotonga to 

Case 

DD 

RMS(m) 

No. of 
DD 

Vava’u 

1509181+ 

W. Samoa 
1527539+ 

Tongatapu- 1 
1605588+ 

Tongatapu-2 

1605588+ 

A 

B 

C 

D 

E 

0.03696 

0.03594 

0.03729 

0.03737 

0.03572 

63,370 

63,370 

70,997 

63,834 

75,797 

0.635 

0.649 

0.630 

0.635 

0.637 

0.105 

0.116 

0.105 

0.105 

0.092 

0.575 

0.541 

* Phase center to pi 

tiase center, not reduced to geodetic mark 

Note: Case to case orbits change by up to 3 meters 
Tongatapu 1 synchronized < 1 sec 
Tongatapu 2 synchronized < 10 sec 
DD =» Double Difference 
RMS = Root Mean Square 


( 


TABLE 7. OTHER COMPARISONS 

Onsala/Tromso Baseline: 

CaseE: 1406186.915m 

Case A orbit, fixed: 1406186.918m 

CASA UNO 419: 1406186.900 

Three-dimensional Onsala/Tromso: 



Case E 

Case A orbit fixed 

Week 419 

Ax 

-1267718.414 

-1267718.350 

-1267718.396 

Ay 

9692.512 

9692.563 

9692.560 

Az 

608405.882 

608406.024 

"608405.885 


\ 











Figure 1. Southwest Pacific Network: 1988 
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Nearly global distribution of GPS data have been collected to support several campaigns in recent years, 
e.g., the CASA UNO and GOTEX experiments. Such campaigns are becoming increasingly characterized 
by the use of a variety of receiver types as the abundance of new receivers increases. It is not uncommon 
for the global data set collected at the CIGNET sites to be augmented with data from global data outside the 
geographical area of the campaign. As the networks have become more global and as the number of avail- 
able GPS satellites has grown, both opportunities and difficulties have been presented to the analyst com- 
munity. The opportunities presented by the global data sets include an improved understanding of the nature 
of nongravitational forces acting on the satellites. Such improved understanding could contribute to con- 
comitant improvements in baseline determinations through improved ephemerides. Results obtained from 
specific campaigns are described based on double difference measurements. The results show evidence of 
unmodeled effects on the orbits and possible sources are discussed. Evaluation of the orbits is conducted 
through orbit comparisons and assessment of baselines derived from the use of those orbits, where possible. 
For example, a recent comparison of the Onsala/Tromso baseline with recently released VLBI results shows 
agreement at a level of better than 2 parts in 10 8 on the 1406 km baseline. 


INTRODUCTION 


The importance of accurate orbits to support applications of GPS is widely recognized. Yet the orbit accu- 
racy ir one of the more difficult characteristics to ascertain in an absolute sense. For the most part, orbit ac- 
curacy is assessed in geodetic applications (e.g., positioning) by the ability of the orbit to support the desired 
accuracy in the application. In fact, for those applications, the orbit could exhibit large errors in an absolute 
sense and still be tolerated if the orbit is able to produce the desired geodetic accuracy. Even though the orbit 
accuracy does not always map in a one-to-one manner into geodetic accuracy, it is a readily demonstrated 
fact that high orbit accuracy is a necessary condition for high geodetic accuracy from satellite techniques. 

The orbital position accuracy depends on a variety of factors, including the fidelity of models describing 
both the forces acting on the satellite as well as the kinematic models. In a complex way, the model fidelity 
interacts with other characteristics such as the data span, or arc length, and the measurement precision (and 
accuracy) to influence the orbit accuracy as noted, for example, by Lichten [1990] and Schutz, et al. [1990]. 
In addition, software accuracy is a factor that cannot be ignored. 

From the standpoint of assessing the orbit accuracy, several criteria can be used, such as: 

• direct comparison of independently determined orbits using either the same or different data sets 

• direct comparison of independent software models 

• direct comparison of geodetic results obtained by independent sources 

This paper attempts to assess the GPS orbit accuracy for a specific period of time using these criteria. A com- 
plete analysis with all criteria has not been conducted, and unequal emphasis is placed on each of the pre- 
ceding criteria in the following sections. 


DATA 


In the summer of 1989, a crustal motion campaign with GPS was conducted in the Southwest Pacific. The 
data collected from the CIGNET stations was augmented with data from several additional sites as well as 
the Defense Mapping Agency sites. The fiducial sites are given by Abusali, et al. [1990]. The campaign took 
place in GPS Weeks 496-500, however, the results given in this paper concentrate primarily on Weeks 499 
and 500. The coordinates of the fiducial sites were adopted in a reference frame defined by Satellite Laser 
Ranging (SLR) and Very Long Baseline Interferometry (VLBI) as described by Schutz, et al. [1990]. 

A mixture of receivers was used at the fiducial sites. For example, the CIGNET sites used MiniMac 2816AT 
receivers, except for the use of TI-4100 at Kauai, Yellowkinife, Onsala and Tromso. The extended network 
used TI-4100 at Orroral, Huahine and Ft. Davis, a MiniMac at Kwajalein and a Rogue at Wellington. The 
DMA sites all used TI-4100 receivers. The field sites in the Southwest Pacific used TI-4100 and Trimble 
receivers, but these data were not used in the analysis described in this paper. 

The phase and pseudo-range data were collected at 30 second intervals and were used in a double difference 
mode to eliminate most of the influence of satellite and receiver clock variations. Pseudo-range data were 
used to ensure time tag accuracy of one microsecond with respect to GPS time. 


MODELS 


The force and kinematic models used were generally consistent with the IERS Standards [McCarthy, 1989], 
with a few exceptions. These standards have been primarily adopted for the Lageos SLR satellite, but can 
be readily adjusted to accommodate GPS. The primary gravity differences between the Standards and the 
GPS model is the truncation of the gravity field at degree and order 8 for GPS. The nongravitational force 
model used ROCK4 with an estimated scale parameter and a y-bias parameter. A zenith delay was estimated 
every 3 hours at each tracking site. Ionospherically corrected phase measurements were used in the analysis. 

The results have been obtained using the Multi-Satellite Orbit Determination Program (MSODP), a multi- 
satellite derivative of the single satellite University of Texas Orbit Processor (UTOPIA). The data were pre- 
processed using IBM PC, VAXstations and a VAX 8200. Final results have been, for the most part, gener- 
ated on a Cray X-MP/24 computer. 


RESULTS 


The results obtained used all available Block-I satellites (PRN 3, 6, 8, 9, 11, 12 and 13) and one Block-II 
satellite (PRN 14). The statistics for one week arcs coincident with Weeks 499 and 500 are given in Table 
1. As shown in this table, almost 100,000 double difference phase measurements were used in the solution. 
The solution for Week 500 shows fewer measurements because not all available data have been merged into 
the solution process at the present time. 

The resulting Block-I Earth-fixed ephemerides for Week 499 were compared with ephemerides produced 
by the Defense Mapping Agency and the Naval Surface Warfare Center. PRN 8 was not available in the 
DMA ephemeris set. The results of a direct comparison, without any adjustments in the respective ephemer- 
ides, are summarized in Table 2. The differences for PRN 3 arc shown in Fig 1 for radial, along-track and 
cross-track components. As shown in Table 2 and illustrated in Fig. 1, the radial component exhibits a con- 
sistent mean value of 1 .3 m which is caused by the different values of GM used in WGS-84 (398600.5 km 3 / 
s ) and the IERS Standards (398600.44 km /s 2 ). The cross-track differences suggest a reference frame 
source and the along-track components are expected to be dominated by the force model. 


Table 1. Solution Summary 



Week 499 

Week 500 

Number of satellites 

8 

8 

Number of stations 

18 

16 

Number of DD* 

97828 

82702 

RMS of DD (cm) 

3.0 

3.0 

Number of sites estim. 

9 

7 


*DD: Double difference 


Table 2. Week 499 Ephemeris Comparison with DMA, Mean (RMS) 


Satellite PRN 

Radial (m) 

Along-Track(m) 

Cross-Track(m) 

3 

-1.3 (1.6) 

-1.9 (3.1) 

0.5 (4.6) 

6 

-1.3 (1.6) 

-0.3 (2.2) 

0.5 (4.2) 

9 

-1.3 (1.7) 

-0.8 (2.4) 

0.5 (3.6) 

11 

-1.3 (1.4) 

1.0 (1.8) 

0.0 (2.5) 

12 

*1.3 (1.3) 

-2.1 (2.8) 

-0.1 (3.8) 

13 

-1.3 (1.4) 

0.5 (1.5) 

-0.2 (2.4) 

Note: 

the value in parenthesis is the RMS 




The ephemeris differences illustrated in Table 1 and 2 as well as Fig. 1 were further investigated by attempt- 
ing to ascertain the level of software compatibility. In this investigation, the DMA/NSWC ephemerides 
were used as “observations” to MSODP and the satellite position, velocity, ROCK4 scale and a y-bias were 
adjusted for each satellite. While the adjustment reduces the level of differences, residuals at the meter level 
remain as illustrated in Fig. 2, the along-track component for PRN 3. The source of the remaining differenc- 
es are expected to be solely caused by software and model differences, including the reference frames and 
the treatment of nongravitational forces. 

A test of model fidelity was conducted in which the solution for Week 499 summarized in Table 1 was pre- 
dicted forward in time and compared with the ephemeris estimated for Week 500. The results obtained were 
generally consistent with those reported for the CASA UNO experiment [Schutz, et al., 1990] in which the 
satellites experiencing eclipsing showed significant degradation in the along-track prediction performance, 
a strong indicator of model fidelity. The effects of eclipsing have been recently examined by Feulner [1990] 
and Vigue [1990]. Feulner investigated the influence of numerical integrator performance across the discon- 
tinuous shadow boundary and modified the MSODP fixed step integrator to handle the effect. While his re- 
sults showed that the effect was significant, it did not remove the observed characteristics with real GPS 
data. In the second study, Vigue solved the partial differential equations for heat conduction and radiation 
for simplified spacecraft models simultaneously with the orbit equations. She showed that distinctly differ- 
ent temperatures can exist on the solar panel surfaces which can produce a several meter orbit effect over 
one week. Additional investigation is in progress. 

An important question from the geodetic standpoint is the influence of the previously described orbit differ- 
ences on the recovery of geodetic parameters. Detailed comparisons of such recovery using the ephemerides 
described in this paper (CSR and DMA/NSWC) are in progress using the field data in the Southwest Pacific. 
The Onsala/Tromso baseline has been repeatedly used as a verification test, especially since it was measured 
in 1989 by a VLBI visit to Tromso. Using the double difference phase ephemerides described in this paper, 
the estimated baseline (without the application of any a priori constraints) agrees with the VLBI result at 
the 1-2 parts in 10 8 l evel on the 1400 km baseline. This result suggests that the ephemeris differences de- 
scribed above are not a limiting factor at this level. The possibility that the differences are a limiting factor 
in achieving a part in 10 9 geodesy remains to be determined. 
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Fig. 1. Week 499 Differences Between DMA/NSWC and CSR Ephemerides (PRN 3) 
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In order to establish a relationship between the SLR/VLBI and WGS-84 reference frames, GPS tracking 
data collected from stations defined in these two frames during a common period have been processed. Pre- 
vious investigations have used common pseudo-range data from the OCS Monitor Stations and augmented 
CIGNET sites to determine a longitude rotation of about 1 m. For the investigation reported in this paper, 
the CIGNET sites have been augmented with TI-4100 data from five DMA sites and were processed in a 
double differenced phase mode. The coordinates resulting from adjustments of the DMA sites into the SLR/ 
VLBI system used for the CIGNET sites with a priori WGS-84 coordinates, were used to determine the sev- 
en parameter transformation set. The characteristics of the transformations have been investigated for se- 
lected time periods in 1989 when DMA data were available to support particular experiments. 



INTRODUCTION 


Several reference frames are in use for geodetic applications, including the WGS-84 [Decker, 1986; White, 
et al., 1989] and those frames based on Satellite Laser Ranging (SLR) and/or Very Long Baseline Interfer- 
ometry (VLBI). The reference frame defined by the International Earth Rotation Service (IERS) is based 
primarily upon SLR and VLBI [Boucher, 1990], known as the IERS Terrestrial Reference Frame (ITRF), 
and is one of the important products of the IERS. 

Since the GPS ephemerides broadcast from the satellites are in the WGS-84 frame, many applications have 
adopted that frame as the most convenient or appropriate. On the other hand, most applications that require 
comparison with external techniques, such as VLBI and SLR, or that will continue time series of geodetic 
measurements that have been obtained by those techniques, will adopt the ITRF or another system based on 
SLR and VLBI. 

Recent results have shown agreement between the SLR and VLBI reference frames at the less than 3 cm 
level [Ray, et al., 1990]. These results were obtained using sites which have been occupied by both SLR and 
VLBI at different times as well as those sites where both techniques operate within a few kilometers on a 
regular basis. This centimeter level demonstration of accuracy supports the notion that a reference frame 
based on the combined SLR and VLBI systems is capable of supporting high accuracy geodetic applica- 
tions. 

The relationship between the WGS-84 and SLR/VLBI reference frames is required for a variety of applica- 
tions, especially those applications in which results obtained in one reference frame are to be compared or 
applied in the other reference frame. Previous studies in which this relationship was investigated include 
Abusali, et al. [1989] in which the smoothed pseudo-range data from the Operational Control Segment 
(OCS) were used in combination with data from the Cooperative International GPS Network (CIGNET) to 
determine the relation between the reference frames. In this study, the data from the OCS were combined 
with pseudo-range data from the CIGNET. By fixing the OCS coordinates to the WGS-84 values, the orbits 
and the coordinates of the CIGNET sites were estimated using WGS-84 force model parameters as well. 
Conversely, the CIGNET sites were fixed to their coordinates in an SLR/VLBI reference frame and the OCS 
sites were estimated simultaneously with the orbit parameters. With the estimated coordinates in the respec- 
tive systems, the seven parameter transformation (translation, rotation and scale) was determined. The ma- 
jor conclusion of this study was that there appeared to be an approximately 60 milliarc second longitude 
rotation between the WGS-84 and the SLR/VLBI reference frames, while little consistency was observed 
in the other six parameters depending on the estimated parameters. In any case, the observed longitude shift 
was within the uncertainty assigned to the longitude origin by the Defense Mapping Agency of 200 milliarc 
seconds (DMA, 1987). 

The purpose of this paper is to further examine the relationship between the WGS-84 and a SLR/VLBI ref- 
erence frame using double differenced phase data and a different data set than used previously. In this case, 
data from the five DMA sites were used to determine position vectors relative to the CIGNET sites whose 
coordinates are described in an SLR/VLBI reference frame. By comparing the resulting coordinates with 
those obtained in the WGS-84 system, the seven parameter transformation parameters can be obtained. The 
resulting parameters are described and investigated in the following sections. 


DATA 


While the establishment of the relation between the reference frames is important for the previously de- 



scribed reasons, the primary motivation for this study was the need for the DMA data to provide tracking 
coverage during the Southwest Pacific ‘89 campaign in which, for various reasons, incomplete geographical 
coverage was available. The primary purpose of this campaign is the study of crustal motions in an area of 
the Southwest Pacific ranging from the Tonga region to New Caledonia. A nearly global set of CIGNET 
sites was further augmented with additional sites in the Pacific as summarized in Table 1, as well as the five 
DMA sites. The campaign took place during GPS Weeks 496-500, however, the results for this paper con- 
centrate on GPS Week 499. 

The coordinates of the CIGNET and augmented network in Table 1 were defined in an SLR/VLB I reference 
frame using results that differed only slightly from those given by Ray, et al. [1990]. The estimation of the 
DMA sites was based on the double differenced phase data formed between the sites given in Table 2. As 
shown in this table, the baselines varied in length from 900 km to 4000 km. 



TABLE 1. Contributing Sites 


CIGNET 

Augmented 

DMA 

Mojave, CA 

Huahine, Fr. Polynesia 

Buenos Aires, Argentina 

Kokee Park, HI 

Onroral, Australia 

Quito, Ecuador 

Westford, MA 

Tanna 

Smithfield, Australia 

Richmond, FL 

Western Samoa 

Hermitage, England 

Yellowknife, NWT 
Wettzell, FRG 
Onsala, Sweden 
Tromso, Norway 
Tsukuba, Japan 

Kwajalein 

Bahrain 


TABLE 2. Baselines 

Reference Site 

DMA 

Baseline (km) 

Orroral 

Smithfield 

941 

Onsala 

Hermitage 

1078 

Richmond 

Quito 

2841 

Wettzell 

Bahrain 

4050 

Quito 

Buenos Aires 

4252 


RESULTS 


With the network adopted for this experiment, information about the orbit is given by Schutz, et al. [1990]. 
To summarize, a one week solution was computed in which the DMA sites were treated as unknowns and 
a simultaneous solution for orbits, spacecraft y-bias, ROCK4 radiation pressure scale parameter, double dif- 
ference biases and troposphere parameters were obtained. For this solution, approximately 100,000 double 
differenced measurements were used and these were fit in the estimation process to about 3 cm. 

The resulting estimates for the DMA site coordinates, defined in the SLR/VLBI reference frame, were then 
compared with coordinates determined in the WGS-84 reference frame. As shown in Table 3, the direct 
comparison of the coordinates is at the one meter level in all three components. Using these differences to 
solve for translation, rotation and scale parameters results in differences (after applying the estimated pa- 
rameters) of less than 0.5 m. As shown in Table 3, a longitudinal rotation (32 mas) is the dominant angular 
contribution. Furthermore, the scale difference is approximately consistent with the different values of GM 
(WGS-84: 398600.5 km 3 /s 2 ; SLR: 398600.4404 km 3 /s 2 ). 

The longitude rotation deserves additional comment. The observed rotation requires an eastward rotation of 
32 mas of the WGS-84 to coincide with the SLR reference meridian. However, an eastward 16 mas rotation 
was required to make the VLBI meridian coincide with the SLR meridian. 

In summary, the VLBI, SLR and WGS-84 reference frames are compatible at the meter level. Based on the 
results given in Table 3, an eastward longitudinal rotation of 16 mas is required to make the WGS-84 coin- 
cident with VLBI and 32 mas is required to be coincident with SLR. Since the choice of longitude origin is 
arbitrary, these rotations may not be valid for all SLR and VLBI solutions. 


ACKNOWLEDGEMENTS 


This study was supported, in part, by NASA, NSF and the University of Texas at Austin. Computational 
assistance was provided by The University of Texas System Center for High Performance Computing. The 
Southwest Pacific GPS Campaign has been organized by M. Bevis and F. Taylor, and includes numerous 
contributors whose assistance is greatly appreciated. The authors are grateful to DMA for providing the data 
from their sites. 


REFERENCES 


Abusali, P. A. M., B. E. Schutz, H. Fliegel and A. Liu, Coordinate Reference Frame Compatibility Between 
VLB I/SLR Systems and WGS-84, Proceedings of Fifth International Geodetic Symposium on satellite Po- 
sitioning, Las Cruces, NM, March 13-17, 1989. 

Boucher, C. and Z. Altamimi, The Initial IERS Terrestrial Reference Frame, IERS Technical Note No. 1, 
1989. 

Decker, B. L., World Geodetic System 1984, Proceedings of Fourth International Geodetic Symposium on 
Satellite Positioning, Austin, TX, April 28-May 2, 1986. 

Defense Mapping Agency, Department of Defense World Geodetic System 1984, DMA Tech. Report 8350- 



2-B, September, 1987. 

Ray, J., J. Ryan, C. Ma, T. Clark, B. Schutz, R. Eanes, M. Watkins and B. Tapley, Comparison of VLB1 and 
SLR Geocentric Site Coordinates, Geophysical Research Letters , to appear, 1990. 

Schutz, B., C. Ho, R Abusali and B. Tapley, Precision Orbit Determination for GPS, Proceedings of the Sec- 
ond International Symposium on Precise Positioning with GPS (GPS ‘90), Ottawa, Ontario, September 3- 
7, 1990. 

White, H. L., B. L. Decker and M. Kumar, World Geodetic System 1984: A Modem and Accurate Global 
Reference Frame, Proceedings of the Fifth International Geodetic Symposium on Satellite Positioning, Las 
Cruces, NM, March 13-17, 1989. 


TABLE 3. WGS-84 to VLBI/SLR 


Direct Comparison: 


Component 

RMS (m) 


X 

1.211 


y ' 

1.033 


z 

Seven Parameter Transformation: 

0.894 


Component 

Std. Error 

Origin x (m) 

-0.094 

0.12 

Origin y (m) 

0.281 

0.12 

Origin z (m) 

- 0.436 

0.12 

Scale (ppm) 

-0.222 

0.02 

Rotation x (mas) 

-6.018 

4.40 

Rotation y (mas) 

2.279 

5.22 

Rotation z (mas) 

31.700 

4.67 

Comparison after Transformation: 


Component 

RMS (m) 


X 

0.446 


y 

0.212 


z 

0.107 
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In order to establish a relationship between the SLR/VLBI and WGS-84 reference frames, GPS tracking 
data collected from stations defined in these two frames during a common period have been processed. Pre- 
vious investigations have used common pseudo-range data from the OCS Monitor Stations and augmented 
CIGNET sites to determine a longitude rotation of about 1 m. For the investigation reported in this paper, 
the CIGNET sites have been augmented with TI-4100 data from five DMA sites and were processed in a 
double differenced phase mode. The coordinates resulting from adjustments of the DMA sites into the SLR/ 
VLBI system used for the CIGNET sites with a priori WGS-84 coordinates, were used to determine the sev- 
en parameter transformation set. The characteristics of the transformations have been investigated for se- 
lected time periods in 1989 when DMA data were available to support particular experiments. 


INTRODUCTION 


Several reference frames are in use for geodetic applications, including the WGS-84 [Decker, 1986; White, 
et al., 1989] and those frames based on Satellite Laser Ranging (SLR) and/or Very Long Baseline Interfer- 
ometry (VLBI). The reference frame defined by the International Earth Rotation Service (IERS) is based 
primarily upon SLR and VLBI [Boucher, 1990], known as the IERS Terrestrial Reference Frame (ITRF), 
and is one of the important products of the IERS. 

Since the GPS ephemerides broadcast from the satellites are in the WGS-84 frame, many applications have 
adopted that frame as the most convenient or appropriate. On the other hand, most applications that require 
comparison with external techniques, such as VLBI and SLR, or that will continue time series of geodetic 
measurements that have been obtained by those techniques, will adopt the ITRF or another system based on 
SLR and VLBI. 

Recent results have shown agreement between the SLR and VLBI reference frames at the less than 3 cm 
level [Ray, et al., 1990]. These results were obtained using sites which have been occupied by both SLR and 
VLBI at different times as well as those sites where both techniques operate within a few kilometers on a 
regular basis. This centimeter level demonstration of accuracy supports the notion that a reference frame 
based on the combined SLR and VLBI systems is capable of supporting high accuracy geodetic applica- 
tions. 

The relationship between the WGS-84 and SLR/VLBI reference frames is required for a variety of applica- 
tions, especially those applications in which results obtained in one reference frame are to be compared or 
applied in the other reference frame. Previous studies in which this relationship was investigated include 
Abusali, et al. [1989] in which the smoothed pseudo-range data from the Operational Control Segment 
(OCS) were used in combination with data from the Cooperative International GPS Network (CIGNET) to 
determine the relation between the reference frames. In this study, the data from the OCS were combined 
with pseudo-range data from the CIGNET. By fixing the OCS coordinates to the WGS-84 values, the orbits 
and the coordinates of the CIGNET sites were estimated using WGS-84 force model parameters as well. 
Conversely, the CIGNET sites were fixed to their coordinates in an SLR/VLBI reference frame and the OCS 
sites were estimated simultaneously with the orbit parameters. With the estimated coordinates in the respec- 
tive systems, the seven parameter transformation (translation, rotation and scale) was determined. The ma- 
jor conclusion of this study was that there appeared to be an approximately 60 milliarc second longitude 
rotation between the WGS-84 and the SLR/VLBI reference frames, while little consistency was observed 
in the other six parameters depending on the estimated parameters. In any case, the observed longitude shift 
was within the uncertainty assigned to the longitude origin by the Defense Mapping Agency of 200 milliarc 
seconds (DMA, 1987). 

The purpose of this paper is to further examine the relationship between the WGS-84 and a SLR/VLBI ref- 
erence frame using double differenced phase data and a different data set than used previously. In this case, 
data from the five DMA sites were used to determine position vectors relative to the CIGNET sites whose 
coordinates are described in an SLR/VLBI reference frame. By comparing the resulting coordinates with 
those obtained in the WGS-84 system, the seven parameter transformation parameters can be obtained. The 
resulting parameters are described and investigated in the following sections. 


DATA 


While the establishment of the relation between the reference frames is important for the previously de- 


scribed reasons, the primary motivation for this study was the need for the DMA data to provide tracking 
coverage during the Southwest Pacific ‘89 campaign in which, for various reasons, incomplete geographical 
coverage was available. The primary purpose of this campaign is the study of crustal motions in an area of 
the Southwest Pacific ranging from the Tonga region to New Caledonia. A nearly global set of CIGNET 
sites was further augmented with additional sites in the Pacific as summarized in Table 1, as well as the five 
DMA sites. The campaign took place during GPS Weeks 496-500, however, the results for this paper con- 
centrate on GPS Week 499. 

The coordinates of the CIGNET and augmented network in Table 1 were defined in an SLR/VLB I reference 
frame using results that differed only slightly from those given by Ray, et al. [1990]. The estimation of the 
DMA sites was based on the double differenced phase data formed between the sites given in Table 2. As 
shown in this table, the baselines varied in length from 900 km to 4000 km. 


TABLE 1. Contributing Sites 


CIGNET 
Mojave, CA 
Kokee Park, HI 
Westford, MA 
Richmond, FL 
Yellowknife, NWT 
Wettzell, FRG 
Onsala, Sweden 
Tromso, Norway 
Tsukuba, Japan 


Augmented 
Huahine, Fr. Polynesia 
Orroral, Australia 
Tanna 

Western Samoa 
Kwajalein 


DMA 

Buenos Aires, Argentina 
Quito, Ecuador 
Smithfield, Australia 
Hermitage, England 
Bahrain 


TABLE 2. Baselines 


Reference Site 
Orroral 
Onsala 
Richmond 
Wettzell 
Quito 


DMA 
Smithfield 
Hermitage 
Quito 
Bahrain 
Buenos Aires 


Baseline (km) 
941 
1078 
2841 
4050 
4252 


RESULTS 


With the network adopted for this experiment, information about the orbit is given by Schutz, et al. [1990], 
To summarize, a one week solution was computed in which the DMA sites were treated as unknowns and 
a simultaneous solution for orbits, spacecraft y-bias, ROCK4 radiation pressure scale parameter, double dif- 
ference biases and troposphere parameters were obtained. For this solution, approximately 1 00,000 double 
differenced measurements were used and these were fit in the estimation process to about 3 cm. 

The resulting estimates for the DMA site coordinates, defined in the SLR/VLB1 reference frame, were then 
compared with coordinates determined in the WGS-84 reference frame. As shown in Table 3, the direct 
comparison of the coordinates is at the one meter level in all three components. Using these differences to 
solve for translation, rotation and scale parameters results in differences (after applying the estimated pa- 
rameters) of less than 0.5 m. As shown in Table 3, a longitudinal rotation (32 mas) is the dominant angular 
contribution. Furthermore, the scale difference is approximately consistent with the different values of GM 
(WGS-84: 398600.5 km 3 /s 2 ; SLR: 398600.4404 km 3 /s 2 ). 

The longitude rotation deserves additional comment. The observed rotation requires an eastward rotation of 
32 mas of the WGS-84 to coincide with the SLR reference meridian. However, an eastward 16 mas rotation 
was required to make the VLBI meridian coincide with the SLR meridian. 

In summary, the VLBI, SLR and WGS-84 reference frames are compatible at the meter level. Based on the 
results given in Table 3, an eastward longitudinal rotation of 16 mas is required to make the WGS-84 coin- 
cident with VLBI and 32 mas is required to be coincident with SLR. Since the choice of longitude origin is 
arbitrary, these rotations may not be valid for all SLR and VLBI solutions. 
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Direct Comparison: 
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RMS (m) 
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z 

Seven Parameter Transformation: 

0.894 


Component 

Std. Error 

Origin x (m) 

-0.094 

0.12 

Origin y (m) 
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- 0.436 
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Scale (ppm) 
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2.279 

5.22 

Rotation z (mas) 

31.700 

4.67 

Comparison after Transformation: 


Component 

RMS (m) 


X 

0.446 


y 

0.212 


z 

0.107 




To appear in the proceedings of the Fifth International 
Symposium on Satellite Positioning 


COORDINATE REFERENCE FRAME COMPATIBILITY BETWEEN 
VLBI/SLR SYSTEM AND WGS 84 

P.A.M. Abusali and Bob E. Schutz 
Center for Space Research 
The University of Texas 
Austin, Texas 78712 

Henry Fliegel and Anthony Liu 
Aerospace Corporation 
2350 East El Segundo Boulevard 
El Segundo, CA 90245 


ABSTRACT 

The tracking sites of the Cooperative International GPS Network (CIGNET) are located near VLBI 
antennas, hence the coordinates of the GPS receiver antennas can be established in a VLBI reference 
frame. Since the GPS dynamics are described in a center of mass system, ties between the VLBI and 
SLR frames are used to define the reference frame. Precise ephemerides of GPS satellites, computed 
using CIGNET data, are then described in the VLBI/SLR frame which results in estimated user 
positions referred to this frame. Some GPS users, as in surveying applications, would prefer to 
compute geodetic quantities in the WGS 84 reference frame. Hence, transformations between these 
two reference frames are needed. 


The Operational Control Segment (OCS) of GPS operates five globally distributed permanent 
tracking sites whose coordinates have been determined in the WGS 84 frame. Processing of tracking 
data from these sites in conjunction with data from the CIGNET stations facilitates the establishment 
of the transformation between the WGS 84 and VLBI/SLR reference frames. Pseudo-range data 
from OCS stations and CIGNET stations were processed separately and the orbits of all satellites 
differ by less than eleven meters in an RMS sense, which is pardy due to the inconsistency between 
the reference frames. The data sets from the two networks were combined and one set of stations 
were estimated with respect to the other set Using the estimated and a priori coordinates of one set 
of stations, a seven parameter transformation between the two systems was computed. The results 
indicate that there is a small longitude shift in the orientation of the two frames. 


INTRODUCTION 

In late 1986, the University of Texas Center for Space Research (CSR) under contract from the U.S. 
Geological Survey (USGS) Menlo Park, California, established a regional GPS tracking network in 
collaboration with the National Geodetic Survey (NGS) and the Texas State Department of 
Highways and Public Transportation (SDHPT). The tracking sites were at Mojave, California, 
Westford, Massachusetts, and Austin, Texas and used TI— 4100 dual-frequency receivers. The 
receiver at Mojave is owned by CSR and has been operated by NGS, the receiver at Westford is 
owned and operated by NGS and the receiver at Austin is owned and operated by SDHPT. The data 
have been routinely processed at CSR, and precise ephemerides of the GPS satellites have been 
computed weekly. These ephemerides have been distributed to the participating agencies and other 
interested institutions. 
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Since early 1987, several other stations have been added to the network by cooperative agreements 
with NGS in an effort to develop a global network. NGS installed its own receivere at additional 
U.S. sites and made collaborative arrangements with foreign agencies to collect GPS data from their 
stations. The network is known as the Cooperative International GPS Network (CIGNET) and the 
currently contributing stations and the operating agencies are listed in Table 1. Data from these 
stations have been processed by CSR to produce precise ephemerides. All the tracking data and 
computed ephemerides have been archived at NGS. 

* “ e CIGNET antennas were installed in the vicinity of already existing or newly established very 
long baseline interferometiy (VLBI) sites. Consequently, the coordinates of the CIGNET antennas 
were determined with respect to the VLBI site coordinates using conventional survey methods and 
GPS techniques [Mader, 1987a,b]. Since the VLBI determined station coordinates are insensitive to 
the Earth’s center-of-mass location, whereas satellite laser ranging (SLR) determined station 
coordinates are influenced by the center-of-mass location, the reference frame adopted for the GPS 
antenna coordinates is a VLBI/SLR-defined frame, determined by the procedure outlined by 
Murray et al., [1988]. As a result, the precise ephemerides computed using these station coordinates 
are described in the same VLBI/SLR reference frame. Further, the gravity field and the Earth’s polar 
motion and UT1 used in orbit computations at CSR are those determined using analysis of laser 
range data from Lageos satellite and reported to the International Earth Rotation Service. 

The GPS satellites broadcast ephemeris messages referred to the WGS 84 system [Decker, 1986], 
and hence, the users of these ephemerides determine their coordinates in the WGS 84 system. If the 
precise ephemeris computed by CSR is used in determining station locations, then the solutions are 
referred to the VLBI/SLR frame. In some geodetic applications, as in surveying for example, users 
would prefer to determine coordinates of sites in WGS 84. As a consequence, there is a need to find 
the relation between the WGS 84 and the VLBI/SLR reference systems. An attempt has been made 
to establish a transformation between these two frames. 

REFERENCE FRAME COMPARISON METHOD 

One method of comparing two different coordinate systems is to compare satellite orbits determined 
using data collected from sets of stations referred to the individual reference frames. These results 
will be dependent on the associated force and kinematic models. The nature and level of orbit 
differences would reflect differences that may exist in the reference frames. However, it is difficult 
to quantify the specific differences using this approach. Another method to compare and quantify 
the differences between two reference frames is to find a set of seven parameters — three describing 
the origin, three describing the directional orientation and one denoting the scale difference of one 
frame with respect to other. Such a set of seven parameters can be determined, given two sets of 
cartesian coordinates for a globally distributed set of stations specified in the individual frames. This 
procedure was used to investigate the differences between the WGS 84 and the VLBI/SLR systems. 

data used in obtaining results 

The set of CIGNET stations listed in Table 1 does not have a completely global distribution, since 
the sites are concentrated in Europe and North America. A globally distributed set of stations is 
preferred for reference frame characterization. For this reason, a data set from one of the special 
global tracking campaigns of GPS satellites was chosen. The Jet Propulsion Laboratory helped 
organize a North Andes— Western Caribbean GPS Experiment in which data were collected from a 
set of globally distributed stations for a duration of about three weeks (1/17/88 to 2/5/88). This 
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campaign was named the CASA UNO [Kellogg et ai, 1989J and stations in Australia, New Zealand 
and America Samoa augments the CIGNET stations. 

In order to be consistent with the number of stations considered in the WGS-84 frame, only five 
stations were chosen from among the total of thirteen fiducial stations which participated in the 
CASA UNO experiment. Coordinates of these five stations were detertnined in the VLBI/SLR 
system as listed in Table 2 and shown in Figure 1. The VLBI/SLR system' cooidinates are described 
by Schutz et al„ [1989], Although Ft Davis and Canberra are not among the regularly contributing 
CIGNET stations, they are referred to as CIGNET stations in this paper solely for the purpose of 
identifying the reference frame in which they are considered. 

The Operational Control Segment (OCS) of the GPS operates five globally distributed stations which 
collect GPS satellite tracking data continuously. The control segment uses these data to compute 
satellite orbits and upload the information to the satellites. The GPS satellites then broadcast the 
orbit information as part of the navigation message. The coordinates of the five OCS stations are 
known in the WGS 84 system, and hence, the ephemerides of the GPS satellites computed by OCS 
and broadcast by the satellites are in the WGS 84 system. The OCS stations and their coordinates 
[Smith, 1986] are listed in Table 3 and shown in Figure 1. GPS satellite tracking data from these 
stations for the duration of three weeks coinciding with the CASA UNO Experiment were provided 
to UT/CSR by the Aerospace Corporation for the purpose of a reference frame comparison study. 
However, contiguous OCS data were available for only four days in GPS week 419 (1/17/88 to 

1/21/88) and five days in week 421 (2/2/88 to 2/7/88). This set of data will be referred to here as the 
OCS data in this paper. 

The OCS stations operate STI communications receiver^ which collect dual-frequency, pseudo-range 
and carrier phase measurements. The pseudo-range measurements are smoothed using the carrier 
phase, and ionospheric-corrected range measurements are provided at fifteen-minute intervals. The 
CIGNET stations operate TI-4100 dual-frequency receivers which record pseudo-range and carrier 
phase at some specified intervals. Data are usually collected at thirty-second intervals. Since data 
are processed at CSR in batches of one week arc lengths, the number of measurements from five 
stations for all the GPS satellites is very large. Hence, the 30-second interval pseudo-range data set 
from the CIGNET was compressed by computing "normal points" in fifteen-minute bins at a 
preprocessing stage where erroneous observations are flagged and the ionospheric correction is 
applied. Thus, the OCS and CIGNET data sets are approximately comparable in that they both are 
ionospheric-corrected pseudo-range measurements at fifteen-minute intervals collected from 
individual sets of globally distributed stations. 

ORBIT COMPARISON 

GPS satellite orbits were computed in batches of one-week arcs. As mentioned earlier, since the 
concentration of OCS data was in Weeks 419 and 421, only these two weeks were analyzed. Orbit 
determinations of six of the seven GPS Block I satellites currently operational were performed 
simultaneously in a multi-satellite data processing mode. The satellites considered were PRN-3, -6, 
~ 9 ’ -12 and -13. PRN-8 was not used because of the variations observed in its on-board clock. 

A multi-satellite orbit determination version of UTOPIA [Schutz, 1980] software system was used in 
the University’s dual Cyber 170/750 computer system. The following set of force models were used 
in all the orbit determination runs: 
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• Earth gravity model 

truncated (8 x 8) GEM-L2 [Lerch et al. t 1985] when VLB 1/SLR frame is referenced 
truncated (8 x 8) WGS 84 used when WGS 84 frame is referenced 

• Third body perturbation due to sun and moon using DE200 JPL ephemerides [■ Standish , 1982] 

• ROCK4 model for solar radiation perturbation [ Porter , 1976] 

• y-bias acceleration [Fliegel et.al , 1985] 

• solid earth tide perturbation 

• relativistic perturbation 

For the kinematic model, polar motion and UT1 determined by Lageos long arc solution was used 
including the short period variations in UT1. Geometric variations on the station location due to 
solid Earth tides were taken into account. 

The parameter estimation state vector consisted of the position and velocity components of each 
satellite, one scale parameter for each satellite for the ROCK4 model, one parameter for the y-bias 
force on each satellite and two parameters of a clock model for each of the passes defined 
individually for all satellites and stations. The number of these passes depended on the data 
distribution and the clock stability. One of the satellite clocks or a ground station clock was held 
fixed as reference and was not estimated. In addition to these parameters, coordinates of appropriate 
ground stations were estimated in some cases. 

Orbit determination experiments were performed using the above mentioned force and kinematic 
models and the tracking data from the two sets of stations. The data were processed in batches of one 
week arcs for the GPS weeks 419 and 421. In each case a solution was obtained using the least 
squares approach which estimated the states of all the satellites and other dynamic parameters for the 
epoch. First, individual orbit solutions were obtained by processing the data from the OCS stations 
and CIGNET stations separately. To be consistent with the station coordinates, the respective gravity 
fields, namely, WGS 84 for OCS data and GEML2 for CIGNET data were used. The amount of data 
from both the sets of stations were about comparable for Week 421, whereas for Week 419 OCS had 
almost twice the amount of data compared to CIGNET. For this reason the orbits of the GPS 
satellites determined for the Week 421 were compared in Figure 2 and the comparison is 
summarized in Table 4. It can be observed from these figures and the table that the orbits of all the 
satellites computed using the two sets of data are within about 22 meters of each other in all 
components and within 1 1 meters in an RMS sense. Further, there is a mean difference of about 1.3 
meters in the radial components of all satellites. This is believed to be due to the difference in the 
value of GM used in the two different fields, namely, 398600.5 km^/sec 2 for WGS 84 and 
398600.4404 km?l sec 2 for GEML2. The negative mean in the along-track components of all 
satellites is also due, in part, to the different GM values. Several other factors, such as the difference 
in station coordinates, atmospheric correction applied to the data, remaining clock error, difference 
in the amount of data and the geometry of coverage may contribute to the overall differences 
between the two sets of orbits compared here. Further, any difference between the reference frames 
would also contribute to the orbit differences since the comparison is done in the J2000 nonrotating 
frame. 

Data from both sets of stations were then combined and orbits of all the satellites were determined 
simultaneously by processing data from 10 stations. The orbit experiments and the results relevant 
to the discussion here are listed in Table 5. Only in Case-4 and Case-5 were the indicated station 
coordinates estimated and in all other cases only the satellite states and other parameters were 
estimated. Comparing the results of Case-3 with those of Case-6 and Case-7 shows that the 
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observation residual RMS is greater for Case-3 than that for either Case-6 or Case-7. This result 
indicates that there is an inconsistency between the two frames in which the two sets of stations are 
referenced. Note that in Case-3 the station coordinates are not in one consistent frame, whereas in 
all other cases they are. Although the improvement in residual RMS from Case-3 to Case-6 or 

Case-7 is small, the fact that there is an improvement is an indication of coordinate reference frame 
difference. 

REFERENCE FRAME COMPARISON 

In an attempt to quantify the difference between the reference frames and to find a transformation 
between these orbit determination experiments were performed in which coordinates of stations 
referenced in one frame were estimated while fixing the coordinates of stations in the other. These 
experiments are Cases-4 and -5 in Table 5. The estimated adjustments to the station coordinates for 
each of these cases are shown in Table 6 and in Table 7. The formal standard error corresponding to 
all these station coordinate estimates were less than 25 centimeters except for the Canberra estimate 
in Week 419 which had about 41 centimeters. The pseudo-range data from Canberra during Week 
419 exhibited systematic behavior compared to other data sets believed to be caused by clock 
variation. The data were subdivided into a large number passes for clock model parameter 
estimation which caused the higher formal standard error. 

Although the results shown in Table 6 and Table 7 do not reveal a definite pattern in all the 
coordinates, the longitude adjustments show some regularity which can be interpreted as a longitude 
orientation difference between the two reference frames. While interpreting these results it should be 
noted that the pseudo-range data used in obtaining these results have a noise RMS of about 1.5 
meters. Further, the station coordinates themselves in either frame may have some uncertainty. 
While the uncertainty of the CIGNET station coordinates are believed to be a few centimeters based 
on estimates of the VLB I and SLR coordinates, the accuracy of OCS coordinates is unknown. 
Nevertheless, the orbit comparison results and the station estimation results indicate that there is a 
small difference between the two reference frames. 

A transformation between the two frames which will minimize the difference can be obtained in 
terms of a seven parameter set consisting of three parameters for translation to make the origins 
coincide, three for rotation to make the direction parallel and one to make the scales equal in both 
the frames. Such a set is obtained as a least squares solution using as observations the differences 
between the estimated and a priori coordinates of the stations referred to one frame. Thus when the 
CIGNET station coordinates are estimated fixing the OCS stations, the estimated CIGNET 
coordinates are with respect to the WGS 84 frame. Since the a priori coordinates of these stations are 
in the VLBI/SLR frame using these two sets of coordinates for the same set of stations the seven 
parameter transformation can be obtained. Table 8 and Table 9 show the transformation parameters 
obtained using the results of Case-4 and Case-5 respectively. In Tables 8 and 9, AX, A F, A 2 
correspond to the origin translation, A L corresponds to scale and Am, A\j/, Ae correspond to rotation 
of one frame to the other. The quantities RMSX, RMSY and RMSZ are the RMS values of the 
differences between the a priori and the estimated coordinates of one set of stations, in the X,Y JZ 
directions respectively. The subscripts B and A refer to the RMS values before and after the 
transformation is applied to the station coordinates. The transformation is applied as follows: 
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X 

y 


'ax 

AY 

+ (1+AL) 

1 Aco -Ay 
-A© 1 Ae 

X 

Y 

z 


AZ 


Ay -Ae 1 

Z 


where X,Y£ are the coordinates of a station in one reference frame and x ,y,z are the transformed 
coordinates of the same station referred to the new frame. It can be seen from the Tables that the 
differences between the station coordinates are reduced by applying the transformation thereby 
indicating that the computed transformation parameters are meaningful and have significance. 
However the values do have some uncertainties and it is difficult to arrive at one set of parameters 
with increased level of confidence using the current results. 

CONCLUSION 

Pseudo-range observations of GPS satellites obtained from two sets of global stations referred to 
VLBI/SLR and WGS 84 systems, were processed in order to check for compatibility between the 
two systems. Comparison of the satellite orbits computed in the individual systems indicates that 
there are small differences in gravity fields and in the orientation of the reference frames. Major 
effect of gravity field differences appears to be due to the differing values of GM. 

The relation between the two reference frames was considered in terms of the seven parameter 
transformation. Estimated values of the seven parameters, except for longitude rotation, appear to be 
somewhat random in nature but the variations are within die noise level of the measurements 
processed. However, die longitude rotation from one frame to the other (Aco) exhibits consistency 
although the magnitude has some uncertainty. From dre results presented here one can determine 
that there is a longitude shift of about 60 milliarc seconds with an uncertainty of about 15 milliarc 
seconds, between the two frames, the sense of rotation being eastwards from WGS 84 to VLBI/SLR 
frame. It should be noted that this value of longitude shift is within the uncertainty (200 milliarc 
seconds) assigned to the WGS 84 X-axis orientation [DMA, 1987J. 
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TABLE 1. CIGNET GPS TRACKING SITES 


Location 

Receiver 

Clock 

Agency 

Austin, Texas 

TI-4100 

Cesium 

TSDHPT 

Westford, Massachusetts 

TI-4100 

H-Maser 

NGS 

Mojave, California 

TI-4100 

H-Maser 

CSR 

Yellowknife, Canada 

TI-4100 

Cesium 

CGS 

Richmond, Florida 

TI-4100 

H-Maser 

NGS 

Wettzell, West Germany 

TI-4100 

H-Maser 

IFAG 

Onsala, Sweden 

TI-4100 

H-Maser 

NMA 

Tromso, Norway 

TI-4100 

Cesium 

NMA 

Hawaii 

TI-4100 

H-Maser 

NGS 


N ote : IFAG — Institut fuer Angewandte Geodaesie 

NMA — Norwegian Mapping Agency 
CGS — Canadian Geodetic Survey 
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TABLE 4. SUMMARY OF ORBIT DIFFERENCES EPH (OCS) - EPH(CIGNET) 

1 I „ i _ _ * ' ' f " ' 1 ~ 1 i i — — , -i ■ i ■ 


Radial 

A-Track 

C-Track 


PRN-3 PRN-6 


-0.071 


PRN-9 


MEAN (m) 


1.378 

-2.926 

0.114 


RMS (m) 


PRN-11 PRN-12 PRN-13 


Radial 

1.677 

1.519 

A-Track 

7.866 

8.341 

C-Track 

8.313 

5.421 


1.282 

-3.316 

0.016 


1.307 

-4.730 

-0.168 


1.693 

1.521 

4.498 

5.356 

2.130 

6.940 


1.352 

-9.998 

0.012 


TABLE 5. ORBIT DETERMINATION EXPERIMENTS 


Comments 


WEEK 419 


OCS data only (WGS 84) 

CIGNET data only (GEML2) 

(OCS + CIGNET) data (GEML2) 

Fix OCS and Estimate CIGNET (WGS 84) 

Fix CIGNET and Estimate OCS (GEML2) 

OCS + Estimated CIGNET Coordinates (WGS 84) 
CIGNET + Estimated OCS Coordinates (GEML2) 


RMS (M) 

No. OBS 


1.111 

1.972 

1.602 

1.557 

1.540 

1.557 

1.540 

3034 

1771 

4805 

4805 

4805 

4805 

4805 


WEEK 421 

OCS data only (WGS 84) 

mwm 

CIGNET data only (GEML2) 


(OCS + CIGNET) data (GEML2) 


Fix OCS and Estimate Cignet (WGS 84) 


Fix CIGNET and Estimate OCS (GEML2) 


OCS + Estimated CIGNET Coordinates (WGS 84) 

1.590 

CIGNET + Estimated OCS Coordinates (GEML2) 

1.556 


TABLE 6. ESTIMATED ADJUSTMENTS (METERS) OF 
CIGNET COORDINATES FIXING OCS STATIONS 


Station 


Yellowknife 
Kauai/CASA 
Ft. Davis 
Wettzell 
Canberra 


Week 419 (Case-4 A) 


Ht 


-0.822 

-1.386 

-2.690 

-2.071 

-0.195 


Lon 


0.459 

1.787 

1.228 

0.821 

2.175 


Week 421 (Case-4B) 


Ht 


-1.127 


Lat 

Lon 

0.189 

-0.454 

m 

-0.509 

1.447 

-0.301 

0.878 

0.939 

1.221 
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TABLE 7. ESTIMATED ADJUSTMENTS (METERS) OF 
PCS COORDINATES FIXING CIGNET STATIONS 


Station 


Kwajalein 
Diego Garcia 
Hawaii/OCS 
Ascension 
Falcon-AFB 


Week 419 (Case-5 A) 


Ht 


-0.086 

1.014 

-1.402 

-0.354 

-1.869 


Lon 


-0.646 

-2.084 

-1.745 

-0.411 

-1.728 


Week 421 (Case-5B) 


Ht Lat Lon 


-1.422 0.374 -0.823 

-0.232 0.851 -3.102 

-1.256 0.828 -1.143 

-1.226 0.770 0.396 

-1.282 0.714 -1.769 


TABLE 8. SEVEN PARAMETER TRANSFORMATION 
FROM VLB I/SLR TO WGS-84 


Parameter Week 4 1 9 Week 42 1 


AZ (m) 

-0.204 

-0.611 

AL (j5p>w>) 

-0.248 

-0.128 

A© (rad) 

-0.075 

-0.063 

Ay (rad) 

-0.014 

-0.015 

Ae (rad) ' 

-0.021 

-0.018 

RMSXr (m) 

1.606 

1.573 

RMSYg (m) 

1.398 

0.799 

RMSZg (m) 

1.235 

1.093 

RMSX A (m) 

0.522 

0.753 

RMSY A (m) 

0.530 

0.519 

RMSZ a (m) 

0.421 

0.288 


TABLE 9. SEVEN PARAMETER TRANSFORMATION 
FROM WGS 84 TO VLBI/SLR 

Parameter 

Week 419 

Week 421 

AX (m) 

0.374 

0.585 

AY (m) 

0.877 

0.495 

AZ (m) 

0.404 

0.664 

AL (ppm) 

-0.081 

-0.171 

A© (arc sec) 

0.041 

0.040 

Ay (arc sec) 

0.008 

0.014 

Ae (arc sec) 

-0.009 

-0.003 

RMSX b (m) 

1.236 

1.718 

RMSYg (m) 

1.466 

1.247 

RMSZ B (m) 

0.569 

0.599 

RMSXj | (m) 

0.628 

0.784 

RMSY a (m) 

0.767 

0.675 

RMSZ a (m) 

0.248 

0.255 
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Datenalned by SI* and VLBI Techniques 
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C. Ms, J.U. Ryan, 6, T.A. Clark (NASA/CSFC, Code 926.9 
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R.J. Eanes , M.M. Watkins, B.E. Schutz, 6 B.D. Tapley 
(Univ. of TX, CSR/WRU402 , Austin, TX 78712; SPAN 
UTS PAN : ; UTCSK : : <las t na«>) 

NASA's Crustal Dynanics Project has devoted a 
significant portion of its mobile SLR and VLBI 
observing programs to collocation activities. These 
permit a direct comparison of the geocentric 
coordinates determined by the two independent 
techniques. A SLR solution (rms residual 4.4 cm) was 
produced at CSR using 16 months of Lageos full-rate 
data spanning 1988. A VLBI solution (rms residual 1 4 
cm) was generated at CSFC using all data from 1979 
through 1989; site coordinates and velocities were 
adjusted as global parameters. The coordinates for 18 
pairs of collocated stations were compared using local 
survey ties to translate the VLBI reference points to 
the matching SLR points. The lengths of the local ties 
range from about 14 m co 387 km. Formal standard 
errors were used for weighting all quantities, 
including the ties. 

After allowing a seven-parameter frame transformation, 
the weighted rms differences are: 17.5 mm, 26.4 mm, and 
23.3 mm for X. Y, and 2 components, respectively. The 
most significant outlier is Medicina, Italy (residual 
of 189 mm), for which there are only 11 mobile SLR 
passes. The overall reduced chi 2 is 11.7 which implies 
that the formal errors for each technique should be 
increased by a factor of 2.4, assuming an equal 
distribution of the errors and no added contribution by 
the local survey ties. However, it is likely that the 
survey ties are a significant source of discrepancy for 
at least a few of the comparison sites. 

The scale of the VLBI frame is fractionally smaller 
than that of the Slit frame by about 10 ppb. Formally, 
the scales of each system are thought to be accurate to 
about 2 ppb. The observed scale difference cannot be 
accounted for by any known or suspected defect in 
either technique. 
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Frame of the Crustal 


The VLB: Terrestrial Reference 
DynarrJ.es Project 

£_ila (Geodynamics Branch, NASA/Goddard Space Flight 
Center, Greenbelt, MD 20771; 301-286-3992) 

The terrestrial reference frame realized by the Crus- a. 
Dynamics Project using VLB: data through the end = £ 
1989 includes 80 points in North America. Europe, 
Africa, Asia, the Pacific and Australia. The origin of 
e frame is within 1 mm of the origin of ITRF89, th<» 
orientation agrees with ITRF89 within 2 mas, but -.he 
scale is determined by the VLBI models. 3-dimensior.a 1 
positions and velocities have been determined. The 
distribution of root-sum-squared 1-sigma formal errors 
in Cartesian position at 1990.0 is: < 1 cm - 14 sites 
1-3 cm - 45 sites, 3- 10 cm - 17 sites. > 10 cm - 3 
sites Comparison with SLR positions suggests that the 
.ormal errors are perhaps too small by a factor of 2-3 
The distribution of 1-sigma formal errors in horizontal 
velocity is: < 1 mm/ yr - 20 sites. 1-3 mm/yr - 22 
sites, 3-10 mm/yr - 9 sites. The remaining sites have 
i sufficient data to compute meaningful velocities. 
The reference frame positions and associated earth 
orientation time series are available electronically 

Svstem ■ wT* 1 Dyn *'* ic5 Project Data Information 

System and will be updated periodically. 
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Reference Frame for Monitoring Crustal 
Deformation wuh the Global Positioning System 

^£ric H , a ^riP i RAY ’ R ° BERT w - *INC (Dept, of Earth, Atmos- 
fis m 7 % C So Sc, t n ^L MrT 54 * 52 1 ’ Abridge, MA 
Pptcd i mhrn@benioff.mit.edu), and 

P (Sif^ C R^ AN (Schoo ‘ of Applied Science, University of 
Canberra, Bclconnen. A.C.T. 2616, Australia) 
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COMPARISON OF SLR /VLBI AND WGS-84 REFERENCE 
FRAMES ' 

P.A.M. Abusali (Center for Space Research, The University of Texas 
at Austin, Austin, Texas 787 1 2) 

B. E. Schutz and B. D. Tapley 

The CIGNET tracking data, combined with data collected during 
the Southwest Pacific Campaign (summer 1989) from the field sites 
and a few additional sites collocated near SLR/BLVI sites, provide a 
global set of GPS data. Coordinates of these stations are known 
(some are estimated) in the SLR/VLBI reference frame. GPS data 
collected by the Defense Mapping Agency (DMA) at the five DMA 
sites, whose coordinates arc known in the WGS-84 frame, are 
combined and processed with the data from augmented CIGNET 
data. The two reference frames mentioned above exhibit 
inconsistencies between them at a level of about 1 m as determined 
by earlier analyses. A 7-parameter transformation betwen the two 
frames is determined using the above data set and compared with 
earlier results. Several issues arc addressed in regard to unifying the 
various reference frames. 
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Combining and Comparing GPS Global Results with 
Other Space Geodetic Techniques 

T P Yunck (both at Jet Propulsion Laboratory 
3540614) InSt,tUtC ° f Tcchnol °gy. Pasadena. CA 91109; 8 18- 

S C Pogorelc (University of Colorado, Dept, of Aerospace and 
Engineering Sciences, Boulder, CO, 80307) P 

The GJobal Positioning System (GPS) will ultimately include 24 

° r n" 1 P 3neS Wi ‘ h ,2 * hr Presently more than 

a dozer GPS satellites are in orbit and available for tracking. Even 

3 i Pan,a ^ onst f , l ation and limited ground coverage, present- 

fn,ercontfai n H 0 f g ‘° bal V T m SUCh as carth dentition and 
ntercontinema 1 baselines are showing consistency with VLBI at the 

level of - 5 cm. Our studies indicate that by early 1991 GPS-based 
estimates for these global parameters arc expected to improve to the 
rw-cm level. With .the : full CPS cons.ell.don and S,tS 
network, accuracy should eventually reach the 1cm level. 

GPS data can be used to refine estimates of the eeocenter earth 
onentatton. global and regional baselines, as well a®prov,de a S 

TOP^^E?DON°rK ,r M k i? e sys,em for ™ si °" 5 « 

Heiminldnn of ,? Nt GPS 0ffcr " nMn unic l uc "Polities in 

(~ h™ HowoUrTh^m,''”' incl "‘J m8 ,CT >P<»I resolution 
• P l he mr>st PO^rful tool for studying elobal 

VLm y Lh rp? T P k r0C ,' S “ S Wl " ** "’ e combmod system of SLR 
VLBI, and GPS. The International Earth Rotation Service (IERS) 

SST lifaen « COm,>,mn8 d,ffcrcm dau and dcfmi "8 

consistent w.th the covariance analyses for present-day eipeued 
performance. However expected accuracy enhancements m the near 
HnhTi «!- ret ! u,re ® better understanding of the rclaiionship between 
GPsJeSS °b ta >ned from different space geodetic techniques. 

YLBI an^or SI R W," 3 prion modcl constraints from 
VLbt and/or SLR. We will discuss issues to be resolved when some 
constraints are relaxed so that the full strength of these different 
geodetic techniques can be realized through their combination. 
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GPS-Derived Earth Orientation Paraaeter* Based 
on Both Block I and Block II Sataliitas 

£ B Ssi (Spaca and Gaodaay Branch (K12). 
Haval surface Warfare Center, Dahlqren. VA 
22448-5000 ? 703-663-7200) 9 ' ' 

From 1987 through 1989 Earth orientation 
estimates were produced as a byproduct of the 
GPS precise epheaeris computations for Block I 
satellites. In order to provide more accurate 
estimates for 1990, the Naval Surface Warfare 
center has started solving for Earth 
orientation on a weekly basis using up to 
twelve satellites - 6 Block I satellites and 
the first 6 Block II satellites. Startinq 
Mmi Ues ^ or t * ie two P°l® coordinates and UT1- 
UTC are obtained using the Defense Mapping 
Agency's predictions generated the previous 
and rate corrections for the 
t he two P° le coordinates and for the 
rate of change of UT1-UTC are being estimated 
along with all satellite orbits and clocks and 
all station clocks simultaneously. This 
estimation procedure uses smoothed pseudoranae 
a«« collected et a global network of ten 
tracking stations. The daily GPS pole 
coordinates and the weekly UT1-UTC rates are 
compared against the final IERS values for the 
data available so far this year. Earth 
orientation estimates based on four-day and 
fit spans fer a five-week period in 

a l EO co ”'P ared against final 
IERS and NEOS values. 
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Gravity Probe B u ax Opportunity to Compart GPS and La*er 
Reference Frames 

M J ? T a pley (Hansen labs GP-B, S unfort University, Stanford CA 94305- 
4080; (415)-725-5697; brtnei TAPLEY%RELGYRO<® STANFORD) 

The Gravity Probe B (GP-B) mission, an experiment in banc physics, will fly in the 
mid 1990's. It will carry a GPS receiver and a set of laser retroreflecion aboard a 
drag-free platform in a 600 km. aliunde polar Eanh obit We show that GP-B will 
»1k)wa good comparisonand strengthening of ihe laser and GPS tracking systems and 
coordinate frames. 

Because of the drag-free system, the dynamics of GP-B's motion can be modeled 
to a high degree of accuracy. By moving .long a well-defined trajectory from one set 
of station, to another, GP-B allows precise <fe ■coronation of the station's relative 
positions; and since it carries both tracking systems, it allow, ihi, desenmnrtcn in 
both coordinate systems. This is why represent, an excellent opportunity to compare 
the User and GPS tracking systems and reference francs. 

As a preliminary step in the evaluation of this opportunity, we conducted a 

evaluation of ihe accuracy wuh which Icast-aquares solutions can determine 
ground stations' positions. We sunuUted a one-month are of data, with User ranges 
from pound station, to the satellite and with GPS position, for both the satelhte and 
the ground station. We assumed colocated laser and gps receivers at the ground 
stations. We generated least squares solutions for the orbit eiemenu of the satellite, 
the positions of the ground stations, and a few geodetic parameters. In some cares -re 
used only laser data, in some cases we used only GPS data, and in some cares wt 
used both. In all cases where we used CPS data, we also adjusted the orbits of the 
GPS satclltes. Generally the combined solution is signifigamly better than either Jjta 
set alone: part of this is due to the increased number of measurements, but partis due 
to the fact that the measurements are some what complementary. We present results 
and suggestions for further research. 
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The Doris System : Status and Prospects 

j . P Chass^n g(CLS, 18 av. E. Belin. Toulouse France), A 
n h o® av ' E B e |ir >. Toulouse. France). C 

Boucher (IGN, 2 av. Pasteur, St Mancte, France) M 
Dorrer (ONES), M Lefebvre (ONES). F Nouel (ONES). 

a u ne ^R re 5 ise orbit determination syste m 
by CNEs to meet the high accuracy -1 


wv.v.v^w M; vrvww \W IMCCl U 19 I ll^i 1 ^CCUTaCY ~1J 

centimeters or better- required for oceanographic altimethc 
system is based upon a worldwide network 

“one trsntnnitlinn nn O t ' . 


missions. The system is oasea upon a worldwide networi 
of ground beacons transmitting on 2 frequencies onboart 
Doppler Receiver and a ground segment to control thr 


TOPEX/POSEIDON M992) DORIS is flown on SPOT2 
(launched Jan 22/90) and is operational with an actual 
network of 36 stations. The paper describes the 
system/principle, the actual status and prospects 
Emphasis will be put on the inflight assesment and will 
adress the system evaluation : instrument performances, 

SSPSf r . pi 5 c,s . l0n -- ! iftk bud 9 et - u,tra stable oscillator 
stability, technological parameters, disponibility data quality 
. A'? , i ab, l ,ty -Propagation models improvments 
(ionosphere, troposphere) with DORIS overall system 
performances . control center, data bank, neWork 
management. Explicit relationships to IERS are requiral 
?, nd "lade through volocations. Lifetime, accuracy and 
flexibikty will allow to contribute to density IERS and to 
case the access from remote sites. 
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ABSTRACT 


THE EFFECTS OF THERMAL IMBALANCE FORCES 
ON SIMPLE SPACECRAFT 
by 

Richard Andrew Cook, B.S. 

Supervising Professor: Dr. Bob.E. Schutz 

The increased orbital accuracy requirements of modem space missions have 
prompted the need to study more subtle forces that perturb the motion of a satellite. 
One such perturbation is due to the reradiation of energy from an asymmetrically 
heated spacecraft. The purpose of this investigation is to obtain some fundamental 
insights into this thermal imbalance effect . Unfortunately, modeling this 
phenomena for a real satellite can be quite difficult because the complete 
temperature distribution within the spacecraft must be determined. Consequently, 
this discussion concentrates on a number of simple spacecraft models that are 
somewhat easier to model. These models are similar to real satellites, but do not 
require an enormous amount of thermal analysis. Each case is investigated to 
determine the maximum impact of the thermal force on the satellite orbit. Some 
attention is also spent on understanding what material properties of a spacecraft are 
significant in this phenomena, and how these parameters can be determined through 
an estimation process. 


v 


Chapter 6 


Conclusions and Recommendations 

The objective of this investigation was to study the influence of the thermal 
imbalance force on geometrical figures representing simple spacecraft. In a broader 
sense, however, the work has entailed a study of two distinctly different subjects. 
The thermal imbalance effect is unique because it combines heat transfer theory 
with satellite orbit dynamics and determination through the thermal force, which is 
the result of heat exchange between a satellite and the surrounding space. It has 
been shown that the rate of this energy transfer is dependent upon the temperature 
distribution of the body. Chapter 3 uses a flat plate model to demonstrate this 
distribution in a simple geometric object. The advanced models in Chapter 4 show 
some of the problems that arise in the study of three dimensional bodies. Thermal 
analysis of actual spacecraft is much more complex than the examples considered in 
this investigation. In fact, determining the temperature behavior of different 
geometrical bodies is a significant topic in itself. There are numerous books 
dedicated to this subject with no mention any of the possible space applications. It 
was very tempting to spend all of this report discussing the temperature behavior of 
increasingly complex satellites. 

A concerted effort was made to not acquiesce to this temptation. After all, 
the prime emphasis of this study was to determine the orbit behavior of a satellite 
subject to this perturbation. In order to do this, the simple thermal models have 
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each been implemented into an orbit computation program. These programs are 
specifically set up to determine the worst case influence of the thermal force. The 
assertion made here is that the simplified thermal models and orbit assumptions do 
not obscure the fundamental physics of the effect. Fortunately, this contention is 
supported by the results. In addition to the specific examples treated, an entire 
chapter is specifically devoted to discussing how to use this effect to improve the 
force model for a real satellite. Obviously, the thermal model for such a satellite 
will be extremely complex, but there are some key implementation aspects that are 
common for every situation. The result of this work is a basic introduction to all 
aspects of the thermal imbalance effect The study was split almost evenly between 
orbital mechanics and thermal analysis. As a result, neither topic is covered 
completely , but significant progress in the interdisciplinary aspects has still been 
made. 

Section 6. 1 Conclusions 

There are a few significant conclusions that can be made from this 
investigation. First, it is apparent that the thermal imbalance force can have a 
measurable effect on earth orbiting satellites. This effect is observed regardless of 
the satellite model or thermal analysis techniques that are used. The significance of 
this perturbation depends on the application. It is definitely important for projects 
that require precise orbit determination. Worst case results indicate that meter level 
radial errors can be generated in fourteen days. The estimation of other force model 
parameters can probably absorb some of this, but differences greater than 10 cm 
can still occur. This is obviously a problem for the TOPEX mission, with desired 
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radial accuracies of about 10 centimeters. The transverse orbit effect due to this 
perturbation can range into the tens of meters, which may be significant for some 
applications. 

The thermal imbalance effect can also vary significantly depending on what 
type of satellite is being modelled. The results in Chapter 3 indicate that the thermal 
force generated by a solar panel can be quite important. The temperature imbalance 
in such an object is sizable because one side receives all of the incident radiation. 
In addition, the heat capacity is small because the mass is assumed to be small. The 
thermal imbalance effect is most significant when the heat capacity is small. The 
discussion in Section 3.5 indicated that the temperature in a body with high heat 
capacity is nearly constant. This implies that the thermal forces generated by heavy 
objects is probably negligible. If the main body of a satellite is represented as a 
single thermally connected object, then it does not have to be included in the 
thermal model. On the other hand, if there are thermally isolated objects like 
antennas or retroreflectors, then these objects must all be modelled independently. 

The particular thermal analysis technique used to determine the temperature 
distribution of a satellite is also important. One of the most significant problems 
encountered in this study was how to simultaneously solve the heat equation and 
the satellite equations of motion. Coupled partial and ordinary differential 
equations are especially difficult to treat numerically, and every solution method has 
to use some type of approximation. The strength of the coupling determines the 
usefulness of a particular technique. The method introduced in Chapter 3 is useful 
for cases where the equations are closely related. Unfortunately, this method is 
slow because the ordinary differential equations describing the temperature 
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variations have to be numerically integrated. In addition, the method is not easily 
extendable to more complicated systems. The finite difference method is generally 
more useful for weakly coupled systems because of its speed. However, it is not 
as useful for strongly coupled equations because of stability limitations. This 
method can be used for multi-dimensional bodies, although the computational 
requirements are excessive. 

The discussions in Section 3.5 and Chapter 5 underscore the role of the 
material properties of the body. Both the magnitude and direction of the thermal 
force are dependent upon the values of the model parameters. The emissitivity is 
important because the size of the force depends directly on the emissitivity 
imbalance. The thermal conductivity partly determines the size of the temperature 
imbalance, so it is also significant. The values of these parameters can be obtained 
by estimation, but the amount of computations is significant. It is possible to use 
the scale factor method introduced in Section 5.1, but only as a check on the 
accuracy of the entire model. 

Section 6.2 Recommendations 

This report has been a general introduction to the thermal imbalance effect. 
As a result, it is more a work of breadth rather than of depth. There are still many 
aspects that need to be studied in more detail. It is important, though, to not spend 
too much time investigating subjects that are insignificant. The following 
recommendations are observations about those particular areas that deserve future 
attention. The first order of business in future research should be to use the simple 
models from this report to learn more about the realistic influence of this 
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perturbation. All of the work in this report has been directed towards determining 
the largest possible influence of the thermal force. In order to accomplish this, 
several rather unrealistic assumptions were made. If the sun, earth, and satellite are 
all allowed to move in a more natural way, then significant information would be 
gained concerning the long-term behavior of a real satellite perturbed by thermal 
forces. 

Future research into the thermal behavior of more complex objects is 
probably not warranted. Most real satellites can be separated into !a main body and 
extended solar panels. The results in this report indicate that the thermal force 
generated by the main body is negligible. The thermal imbalance force can be 
completely determined by the force on the solar panels. Of course, the assumption 
here is that the satellite body is one thermally' continuous body. If a particular 
satellite has isolated components that have significant temperature variations, then 
they probably should also be modelled. Parts like this will vary from satellite to 
satellite, so developing thermal models for them should wait until a specific 
spacecraft is discussed. 

This is not to say however, that there should not be time spent developing 
more sophisticated thermal analysis techniques. Both solution methods used in this 
discussion give reasonable results, but they do have some limitations. Emphasis in 
this area should be placed on developing fast and accurate approaches that can be 
used on a regular basis. The finite difference method has the most potential, but the 
stability of the results must be improved. There also exist more advanced ways to 
solve partial differential equations that are not discussed here, and they may prove 
to be more useful. Special consideration should probably be given to modem finite 
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element approaches. If the study of a more complex body is warranted, then these 
methods should definitely be studied. The finite difference technique is applicable 
to multi-dimensional bodies, but the computational requirements are excessive. 

Finally, there should also be some effort made to determine accurate values 
for the parameters used in the thermal imbalance model. One way that this can be 
done is by extensive research into the material properties of a specific satellite. In 
addition, instrumentation on the satellite may be able to provide direct 
measurements of some quantities. Unfortunately, however, some values may be 
difficult to obtain and will only be valid for that one satellite. Consequently, more 
time should be devoted to developing the tools required to estimate these parameters 
from real observational data. This again relates back to the need for faster and more 
accurate solution techniques. The results in Chapter 5 should be closely studied 
before deciding which parameters to estimate. 
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ABSTRACT 


THERMAL IMBALANCE EFFECTS ON 
A GPS SATELLITE 

by 

Yvonne Vigue, B.S.A.E. 


In recent years, the Global Positioning Satellite System orbit 
estimation and prediction requirements have prompted studies in non 
gravitational orbit perturbations which were previously considered 
insignificant. One of these perturbations, the thermal imbalance force, is 
believed to be caused by an unequal temperature distribution throughout the 
satellite. This discussion concentrates on obtaining an accurate thermal 
force model for a GPS satellite in orbit using the Finite Element Method to 
solve the transient heat transfer problem. Two Finite Element Method 
programs are used to determine satellite temperatures and thermal forces 
which are important in the force modelling used in the satellite prediction 
problem. The results are compared with those from previous studies which 
used finite differencing techniques to solve the transient satellite heating 
problem. 



Chapter 4 


Conclusions and Recommendations 

The purpose of this study was to determine the effect of the thermal 
imbalance force on the orbital behavior of a GPS satellite. This 
investigation used the techniques known as the Finite Difference Method 
and the Finite Element Method to determine the satellite surface temperatures 
which are used to compute the thermal imbalance force. It has been shown 
that the thermal imbalance force is caused by the exchange of thermal energy 
between the satellite and its surroundings and is dependent upon the 
satellite's temperature distribution. Heat transfer theory combined with 
satellite dynamics and statistical estimation theory has been used to develop 
a simulation test case of the problem of satellite heating and cooling by 
radiation. Although the main focus of this thesis is the GPS satellite, a 
TOPEX satellite model was used as a basis in order to compare results with 
several previous studies where the Finite Difference Method was used to 
determine the satellite temperature distribution on a TOPEX solar panel 
[Cook, 1989, Rosborough, et. al., 1989, Dennehy, et. al., 1990]. 

Although several assumptions have been made and the satellite thermal 
model has been simplified, the underlying physics of the problem remains 
unchanged. 

-LOG- 
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Section 4.1 Conclusion 

The result of this investigation compares well with that of the work 
conducted by Cook [1989] wherein the finite difference technique was 
used. In this study, two independent techniques, finite element and finite 
difference, have been used to determine the effects of the thermal imbalance 
force on a GPS satellite. It is clear that the thermal imbalance force can have 
a measurable effect on orbits of satellites like GPS and TOPEX for which 
orbit accuracy requirements are very stringent in order that they may be used 
for geophysical applications. The results indicate that centimeter level radial 
errors and meter level transverse errors could result from this perturbation in 
a one week arc for a GPS satellite. The estimation of parameters such as the 
solar radiation reflectivity can absorb some of this enor. 

Clearly the thermal analysis technique used to determine the 
temperature distribution of the satellite solar panel is important. Under 
certain conditions, the results may vary by as much as an order of 
magnitude when such variables as time step is varied. Both the finite 
difference approximation and the finite element method are adequate heat 
transfer solution techniques. Either of these could be used for solving 
satellite solar panel surface temperatures. Both methods have advantages 
and disadvantages. Speed and stability are two advantages of the finite 
difference approximation method. However, a finite element method is 
probably more appropriate for the following reason. This method is capable 
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of implementing varying material properties as a function of nodal spacing. 
This is highly favorable since satellite solar panels are composed of different 
layers of materials with different material properties such as thermal 
conductivity. When comparing the models used in this investigation, the 
author has a preference of the PDE/Protran finite element model over the 
finite difference method. The software developed by Sewell [1985] has 
undergone extensive testing and verification. It also has a system of checks 
which notifies the user when such variables as time steps are not acceptable 
to the algorithm. For small enough time steps, this method should be very 
accurate. A major difficulty in including the thermal imbalance perturbation 
in orbit computation is the need to solve the coupled partial differential and 
ordinary differential equations simultaneously. 

Section 4.2 Recommendations for Further Studies 

An effort should be made to determine the correct material properties 
for a GPS satellite solar array, or the specific satellite of interest Cook 
[1989] determined that the material properties such as surface emissivity and 
thermal conductivity will have a significant influence on the temperature 
distribution. The thermal imbalance force is directly related to e, the surface 
emissivity imbalance, and the heat equation is largely dependent on the 
thermal conductivity, K. Also, as Cook [1989] mentioned, instrumentation 
on board the spacecraft may provide direct measurement of some quantities 
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such as panel surface temperature. All these factors should be considered 
when implementing the thermal imbalance force into the general force model 
of a satellite. An accurate model of the thermal imbalance force for a 
particular satellite will be useful in reducing the orbit estimation errors 
where a precise orbit is critical for mission success. 
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THERMAL IMBALANCE FORCE MODELLING FOR A GPS SATELLITE 
USING THE FINITE ELEMENT METHOD 


Yvonne Vigue 
Jet Propulsion Laboratory 
California Institute of Technology 
Pasadena, California 

Bob E. Schutz 

University of Texas at Austin 
Center for Space Research 
Austin, Texas 

Geophysical applications of the NAVSTAR Global Positioning 
System (GPS) require the capability to estimate and propagate 
satellite orbits with high precision. An accurate model of all the 
forces acting on a satellite is an essential part of achieving high 
orbit accuracy. Methods of analyzing the perturbation due to 
thermal radiation and determining its effects on the orbits of 
GPS satellites are presented. The thermal imbalance force, a 
non-gravitational orbit perturbation previously considered 
negligible, is the focus of recent attention. The Earth’s 
shadowing of a satellite in orbit causes periodic changes in the 
satellite’s thermal environment. Simulations show that mis- 
modelling of the thermal imbalance force over several days 
gives meter level errors for those satellites which are eclipsing 
and can limit accuracy in orbit determination and in estimation 
of baselines used for geophysical applications. 

INTRODUCTION 

To study the thermal imbalance effect, it is important to understand the 
theoretical background and the laws of physics that govern the problem of radiative 
and conductive heat transfer. Radiative heat transfer between a satellite and its 
environment is the basis for the thermal force model. A satellite in Earth orbit is 
continuously showered by radiation, most of which comes from the sun. The thermal 
imbalance force is directly related to the temperature distribution of the satellite in its 
changing environment (Ref. 1). An uneven temperature distribution causes surfaces to 
re-radiate energy at different rates. Most of the thermal force on a satellite originates 
with the solar panels due to their large exposed area and low heat capacity. 

The satellite’s heated body re-radiates energy at a rate that is proportional to its 
temperature, losing the energy in the form of photons. By conservation of momentum, 
a net momentum flux out of the body creates a reaction force against the radiating 
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surface, and the net thermal force can be observed as a small perturbation that affects 
long term orbital behavior of the spacecraft (Ref. 2). The partial differential equations 
and boundary conditions describing the temperature distribution and the heat transfer 
between surfaces, along with the application of the finite element method are 
presented in this paper. A brief description of the statistical estimation technique used 
for studying the effect of the thermal imbalance force on GPS orbits is included. 

RADIATION AND HEAT CONDUCTION FORMULATION 

Two types of heat transfer that affect a spacecraft in orbit are radiation and heat 
conduction. The exchange of energy between the spacecraft and its surroundings is 
described by radiation heat transfer. Conduction is the transfer of heat by molecular 
motion within a solid medium. Figure 1 shows the type of heat transfer that affects an 
orbiting spacecraft. • 



Figure 1 General Heat Transfer Diagram for a Spacecraft 

The rate of radiant energy transfer for a real object is given by Stefan-Boltzmann 
Law (Ref. 3): x 


Er = eaT 4 0) 

where Er i s the energy radiated by a real body per unit area per unit time, 

e is the emissivity of the body, dependent on material properties, 
c is the Stefan-Boltzmann constant, 

T is the temperature. 

By conservation of momentum , the thermal force, or rate of change of momentum for 
a radiating surface element is expressed as (Ref. 1): 

^thermal * * \ dA t^ ... 

3 c (2) 

where t is time and c is the speed of light. 

The differential force must be integrated over the entire surface to determine the 
complete thermal force: 

Ftherma] “ ‘ 3 ^ I E^dA n a 

C (3) 


2 



Clearly, the thermal force cannot be calculated unless surface temperatures of the 
spacecraft are known. 

In general, the temperature at any point within a body satisfies the heat equation: 

KV ! T=pC p ^ 

dt (4) 

where K is the material thermal conductivity, 

p is the material density and specific heat, 

Cp is the specific heat, 

T is the temperature, 
t is the time, 

and ^ is the Laplace operator. 

The solution to this second order partial differential equation requires that the 
coefficients K, p and Cp be time independent. Also, the boundary conditions must be 
specified. The boundary conditions are described by thermal radiation and 
conduction. As given by the conservation of energy principle, the total amount of 
energy coming into a surface is equal to the total amount of energy leaving the surface. 
The boundary condition for the surface of the satellite can be obtained by using this 
condition as 

qin=qout (5) 

where qi n is the amount of incoming radiative energy due to external sources and 
internal conduction and qqut is the amount of radiative energy leaving the specified 
boundary due to re-radiation and conduction. Figure 2 shows the conservation of 
energy principle for the surface of a satellite solar panel. 



Figure 2 Conservation of Energy Diagram for a Satellite Surface 


Using this concept, the boundary conditions for each surface are constructed. The 
incident radiative solar energy received per unit area per unit time by side A and side 
B of the solar panel are represented by h a an( j h b (R e f. 2): 


KA^ = e b aA1 4 b -h b A 


Ka|£ = e a cATt - h a A 


(6 a) 


(6 b) 


The actual amount of incident radiative energy received by the each side of the soiar 
panel is a function of panel orientation and the orbit of the satellite. The subscript a 
represents the left boundary in local coordinates (cold side) and b is the right boundary 
(sunlight side). The term on the left side of the equal sign in equations 6a and 6b is the 
heat flux, energy per unit time per unit area, in the local x-direction, which is 
perpendicular to the solar panel face. The values used for some of the parameters 
described above are shown in Table 1. 


TABLE 1 

THERMAL AND ORBIT PARAMETERS 

Model Parameter 

Value 

Initial orbit radius: 

.26,550,000 meters 

surface emissivity £ a : 

0.78 

surface emissivity £^: 

0.83 

surface absorptivity CC a : 

0.77 - 14.1% (panel efficiency) 

solar panel surface area - A: 

10.832 meters square 

satellite mass - m: 

845 kg 

initial panel temperature (t = 0) : 

300 degrees Kelvin 

Stefan-Boltzmann constant - o: 

5.6699 E-08 Watts/m 2 °K 

speed of light - c: 

2.998 E+08 m/s 

solar constant - \|/ : 

1368.2 Watts/m 2 

total panel thickness (8 layers): 

0.01478 meters = 0.582 inches 
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The finite element method is used to solve the transient heat conduction and 
radiation problem. PDE-Protran, a finite element method program incorporates 
material properties to determine solar panel surface temperatures fRef. 4). Accurate 
and current knowledge of physical parameters such as surface emissivity, thermal 
conductivity, heat capacity and material density is required. These parameters, which 
may change due to degradation of materials in the environment of space (Ref. 5), 
influence the thermal force calculations and have an effect on the prediction and 
propagation of the spacecraft trajectory. 

ORBIT ANALYSIS TECHNIQUE 

In this investigation, the equations of motion for an Earth satellite are assumed to 
include the two body gravitational effect and the thermal imbalance forces only, and 
are given in vector form by: 


r = . + fthgnn 

r 3 m (7 ) 

where ? is the geocentric satellite position, 

I 1 is Earth’s mass multiplied by gravitational constant, 
m is the mass of the satellite. 


and thermal imbalance force perturbing the satellite, ftherm, is computed as: 

f,h«™ = 2 g !t (ebTt-e,H) (8) 

\ 

The effect of the thermal imbalance force on a satellite can be observed by 
comparing the perturbed orbit with the unperturbed two body orbit in time (Ref. 6). 
Since there is no closed form analytical solution for the perturbed equations of motion, 
a numerical integration technique is necessary to solve the ordinary differential 
equations of motion. The perturbed and unperturbed orbits can be set with the same 
initial conditions and then the displacement between them at a given time can be 
observed. A least squares estimation technique (Ref. 7) is used to determine the state 
of the satellite in its orbit at a specified epoch. The initial conditions of one orbit can 
be adjusted at a given time to eliminate the secular divergence between the perturbed 
and unperturbed orbits to observe the periodic behavior. 


DISCUSSION OF RESULTS 

In order to determine the direction and magnitude of the thermal force, the surface 
temperatures were calculated using a finite element method program. Several 
simulations were tested. The data input that was required for the satellite solar panel is 
shown in Table 2. This table lists the material properties for a Block II GPS satellite 
solar panel (Refs. 5, 8, 9, 10, 11, 12). 


TABLE 2 


GPS BLOCK II SOLAR PANEL PROPERTIES 


Panel Laver 

Composition 

Thickness 

(metersi 

Density 

(kg/m^ 

Specific Heit 
(J/kg°K) i 

i 

Conductivity 

(Watts/m°Ki 

fused silica 
coverglass 

0.00749 

2186.622 

753.624 

1.417 

adhesive 

0.00005 

1079.472 

1256.04 

0.116 

solar cell 

0.00025 

2684.84 

711.756 

147.994 

interconnect cell 
adhesive 

0.00018 

1051.793 

1256.04 

0.116 

Kapton cocured 

0.000076 

1 162.508 

1130.436 

0.1506 

graphite epoxy 

0.00019 

2186.622 

1373.27 

0.8706 

aluminum core 

0.00635 

2491.088 

1046.7 

250.966 

graphite epoxy 

0.00019 

2186.622 

1373.27 

0.8706 


The satellites of the Global Positioning System are distributed in six evenly 
spaced orbit planes. The proposed final constellation will consist of 24 satellites at an 
orbit altitude of 20,000 km and orbit period of 12 hours. These satellites are inclined 
in their orbits such that the spacecraft will experience an eclipsing season for only a 
few weeks every year. At all times, there is at least one orbit plane in eclipse mode. 
Eclipsing has a strong effect on the solar radiation environment those satellites. This 
is evident in the temperature of a GPS satellite solar panel over one orbit shown in 
Figure 3. The steady state temperature for the sun-facing side is approximately 317 
degrees Kelvin and the shaded side is 3 1 3 degrees Kelvin. These values compare well 
with the approximate value of 313 degrees Kelvin which was measured on the cold 
shaded side of the solar panel for a GPS satellite (Refs. 8,9). The face exposed to the 
sun was not directly measured and therefore the temperature difference between 
surfaces is not well known, but is believed to be approximately 5 degrees Kelvin (Ref. 
9, 10, 11). The eclipse period shows a slow decline to an panel equilibrium 
temperature of approximately 253 degrees Kelvin. 
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Figure 3. Temperature history simulation for a GPS solar panel 

Modelling the coverglass surface on the sun-facing side accurately has been 
difficult during this study since that information was not readily available. The 
thermal conductivity of the fused silica coverglass is much lower than any of the other 
7 solar panel stack layers (Ref. 10). This layer contributes most of the temperature 
imbalance because of its low thermal conductivity and high thickness as compared to 
the other solar panel layers. For example, two simulations were performed using 
identical solar panel parameters except for different thermal conductivities for the 
coverglass. The results shown in Table 3 describe the steady-state temperatures and 
thermal accelerations that were computed using the specified values for the coverglass 
thermal conductivity. 


TABLE 3. COVERGLASS THERMAL SIMULATION 

Test Case 1 

Test Cast 2 

K = 1.417 W/m°K 

K = 0.04327 W/m°K 

Hot side T a = 317.41 °K 

Hot side T a = 340.30 °K 

Cold side Tb = 313.66 °K 

Cold side Tb= 285.37 °K 

Thermal Acc. = 1.88 E-10 m/s^ 

Thermal Acc. = -8.01 E-9 m/s^ 


In this paper, the reference frame that is used is the spacecraft radial and along- 
track components. The along-track component is sometimes referred to as the 
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transverse or down-track direction and is defined in the direction of the satellite 
velocity vector. Figure 4 shows the radial and along-track, components of the 
acceleration due to thermal re-radiation over one orbit for an eclipsing satellite. These 
also compare well with studies which have shown un-modelled non-gravitational 
forces to cause errors of this magnitude (Ref. 1 1) 


CM 

< 


C/2 



Figure 4. Radial and transverse components for the thermal force over one orbit 

Figure 5 shows the differences between two orbits, one computed using two-body 
effects only and another trajectory was computed with two-body and thermal 
imbalance force for a satellite in an eclipsing orbit during one week. The radial rms is 
0.5 meters and the along-track rms is 5.2 meters. These results were computed using a 
technique similar to the method used to predict satellite orbits based on a set of initial 
conditions and a complete force model of the spacecraft, which could include the solar 
radiation pressure and thermal imbalance force. In this case for an eclipsing satellite 
after seven days, the along-track components differ by approximately 13 meters. 

Figure 6 also shows the differences between two orbits, one with two-body effects 
only and another computed using two-body and thermal imbalance force for a satellite 
which is not in an eclipsing plane. The radial rms is 0.5 meters and the along-track 
rms is 1.6 meters. It can be seen from these results that an eclipsing satellite 
experiences a larger perturbation in the along-track direction over the span of one 
week, than a satellite which is not in an eclipsing' orbit plane. For the non-eclipsing 
satellite case, after seven days, the along-track difference is approximately 5 meters. 
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The next two figures represent the results computed using a least squares 
estimation algorithm in which the simulated observation data contained only the two- 
body gravitational and thermal imbalance re-radiation forces. The force model used 
in the estimation algorithm contained the two-body gravitational force model with a 
solar radiation pressure model to observe the ability of the force model to account for 
thermal imbalance forces which have been difficult to model but exist in the 
observations. The best estimate of the satellite epoch state, in the least squares sense, 
is calculated which includes the satellite position, velocity and a solar radiation 
pressure scale factor. 

Figure 7 shows the orbit fit residuals for a satellite in an orbit plane that is 
regularly eclipsing. The radial rms is 5 centimeters and the along-track rms is 80 
centimeters. After seven days, the along-track orbit error is almost 2 meters. These 
results show that the solar radiation pressure scale factor in the estimation scheme is 
capable of absorbing most of the orbit error due to thermal re-radiation, but not all of 
the orbit error, especially in the along-track direction. 

Figure 8 also shows orbit fit residuals for a GPS satellite, the same estimation 
technique, but the satellite is in a non-eclipsing orbit plane. The radial rms is 9 mm 
and the along-track rms is 17 centimeters. After seven days, the along-track orbit error 
is approximately 40 centimeters. Clearly, the eclipsing of the satellites has an 
influence on the orbit errors, when a thermal re-radiation force is not included in the 
estimation force model. Larger orbit errors are calculated when the satellite is in an 
eclipsing orbit plane. 

Studies have shown that, for eclipsing satellites, the quadratic-like growth in the 
along-track direction can give errors as large as 50 meters after a one week prediction 
(Ref. 13). A one week prediction can be made using the satellite state computed for 
the best least squares estimate in Figure 7 and compared to the best least squares 
estimate for that predicted week. Although that specific test was not computed for this 
paper, errors of this magnitude are possible. 



Figure 7 Orbit fit residuals, with a solar radiation pressure scale factor, eclipsing satellite 


o 



Figure 8 Orbit fit residuals, with a solar radiation pressure scale factor, non-eclipsing satellite 
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CONCLUDING REMARKS 


The current analysis has shown that orbit errors larger than 10 meters occur when 
mismodelling non-gravitational forces such as the thermal imbalance force presented 
here. A finite element method technique has been used to calculate satellite solar 
panel temperatures which are used to determine the magnitude and direction of the 
thermal imbalance force. Although this force may not be responsible for all of the 
force mismodelling, conditions may work in combination with the thermal imbalance 
force to produce such accelerations on the order of l.e-9 m/s^. One possible 
contribution which is currently being studied is the solar panel misalignment, acting 
together with the thermal imbalance force which may account for much of the 
unmodelled perturbations. If sub-meter accurate orbits and centimeter-level accuracy 
for geophysical applications are desired, a time-dependent model of the thermal 
imbalance force should be used especially when satellites are eclipsing, where the 
observed errors are larger than for satellites in non-eclipsing orbits. This force has 
been observed for years, even on such satellites as LANDSAT and is still not 
completely understood. Forces such as the thermal imbalance force are dependent on 
the space environment and especially on such parameters as the satellite mass, cross- 
sectional area and material composition. Unfortunately, these parameters can change 
or degrade with long-term exposure in space. The results obtained using the finite 
element model used in this study agree with the work of others who have conducted 
studies for various spacecraft such as TOPEX and LANDSAT using the finite 
difference technique to determine spacecraft temperatures. 
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ABSTRACT 

CIGNETh 5 ^!? 8 Sy H em (GPS) 0r ^‘ de,crmination “sing data collected by the civilian tracking network 
. used laa Vanety 0f ex P eriments * dong with data collected from special campaigns to augment 
^ ^ ng data ; These experiments have been used to investigate the fidelity of the satellite force and 
kinematic models with the goal of improving the GPS ephemeris accuracy. Software has been developed for the 

Sti^ateirt! ° f ° PS P scud0 ' ran 8 e “d Phase measurements, including procedures for cycle slip corection. A 
low orbit dettmunation program, MSODP, allows for satellite-dependent force models, including high- 

w^ ne^fv alZl^rv , Asse ^ mcnts / orbit aeemacies derived from measurements collected in campaigns 
knLi frim JlR o^m S tJr been , perfor ™ d with a variety of techniques, including determination of baselines 
u C °JT pUtatl0n of longl multi -week arcs and the comparison of predicted ephemerides 

rl^T^ CphCmendCS : 71,656 resu,ts have P rwi dcd indications of inadequacies in the force mcScls 
especially with the nongravitational forces, particularly those aspects associated with eclipsing. 

INTRODUCTION 

SuSdS C a C rr° f aCCUra !f GP ! 0rbitS bas been widely acknowledged by an abundance of papers /l, 2/. The 
required orbit accuracy is dependent on the particular application and is further complicated by the fact that orbit 

SnSt 'Z y hlle *9 orblt accuracy can be described in terms of, for example, the positional accuracy in an 

dopted reference frame, the more commonly used measure of accuracy is the ability of an orbit, or ephemeris to 

leS sfrht CU i C P ^ am ^n* r CS f mati ° n ' 71,6 Iatter a PP lication describes the accuracy in terms of the baseline 
length, such as 1 part in 10* or 1 cm on a 1000 km baseline. 

accuracy of an orbit is dependent on a variety of factors, including the quality of the force and 
kinematic models as well as the precision and geographical distribution of the measurements used to determine the 

totMVal of LSS £ W1) ' *“ “* •"> *« 

° rbit accuracy usin 8 sevcral evaluation procedures. The data used for the study is 

n fi ? m the CASA Un0 experiment fit. Furthermore, this paper provides additional detail to the results 
published in /4/. 


DATA 

bLed on S the U orM^ ri rt nt WaS COnduC,ed from January 18 throu 8 h F « b ruary 5. 1988, using a tracking network 
based on the CIGNET, but augmented with several additional sites in the Pacific. The complete network of stations 

is given in Table 1, which consisted of TI-4100 receivers. Most of the sites were located near VLBI or SLR sites 
and the position vectors of the GPS receiver reference point with respect to the VLBI/SLR reference points has 
b “ S men (e ‘ 8 -’ 1 F ° r a11 of ** s,udies described in this paper, ionospherically corrected measurements 

2 lored bC satisfactorily s P e f ificd a P riori with respect to SLR and VLBI reference points 

a5 "° 6d ™‘ h *« astensk - ™ ese s,tes were treated as unknown and estimated during the solution process. In 

fhoulh ihl ■!? S,ttS f were selected t0 be estimated for the putpose of aiding in orbit accuracy assessment even 
though their coordinates are accurately known with respect to the SLR/VLBI systems 

m SS iT' " C0U “" <, “ 
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TABLE 1 Sites Used in Solution 


CIGNET 

Augmentation 

Mojave, CA 

American Samoa * 

Onsala, Sweden 

Black Birch, NZ* 

Richmond, FL 

Ft Davis, TX * 

Tromso, Norway * 

Kokee Park, HI 

Westford, MA 

Owens Valley, CA * 

Wettzell, FRG 

Tidbinbilla, Australia 

Yellowknife, NWT 


* denotes sites estimated in simultaneous solution with orbit 


FORCE AND KINEMATIC MODELS 

The force models used differed from the International Earth Rotation Service (IERS) Standards /5/ as follows. The 
adopted gravity field was GEM-L2, truncated at degree and order 8, instead of the IERS field, GEM-T1. 
Experiments have shown that the difference between these two fields has little effect in the orbit determination 
I^ocess. The solar radiation pressure model was ROCK4 with an estimated scale parameter and a y-bias parameter 
/6/. A zenith delay was estimated to account for tropospheric variations. The zenith delay was further divided into 
3-hr parameters for each station. 

The adopted coordinates of the tracking sites were based on the orientation and scale of a University of Texas SLR 
solution known as CSR 8901. Since several of the tracking sites were located in close proximity to VLBI sites, the 
Goddard Space Flight Center GLB484 was adopted. Using a set of common sites between SLR and VLBI the 
VLBI system was translated, rotated and scaled into the SLR system HL 

The orbit determination software, MSODP, is a multi- satellite derivative of the single satellite software used at 
CSR for processing of SLR data, UTOPIA. Data screening, editing and cycle slip corrections are accomplished in 
a separate suite of programs designed for efficient and interactive operation with graphical displays to facilitate 
examination of the data. 

ORBIT RESULTS 

The fidelity of the models used were examined by determining the GPS orbits in one week arcs in two consecutive 
weeks, namely, GPS Weeks 419 and 420 during the CASA Uno experiment In these experiments, the estimated 
parameters in Week 419 were used to predict the ephemeris into Week 420. This predicted ephemeris was directly 
compared with the ephemeris determined from the tracking data. Two classes of results were observed as 
illustrated in Figs. 1 and 2. The characteristics illustrated by Fig. 1 were shared by satellites with PRN identity of 
3, 8, 11 and 13 whereas the characteristics in Fig. 2 were common to PRN 6, 9 and 12. The element of 
commonality shared by these two groups of satellites is that PRN 3, 8, 11 and 13 orbits were experiencing 
eclipsing, whereas PRN 6, 9 and 12 were not. As is suggested by Figs. 1 and 2, the eclipsing satellites exhibit a 
quadratic-like error growth in the prediction, whereas the non-eclipsing satellite differences are bounded to about 
± 5 m. 
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Fig. 1. Comparison of Orbits: Week 420 estimated solution vs. Week 419 predicted solution for SV13 
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CONCLUSIONS 

For the CASA Uno data set examined in this paper, the existence of unmodeled effects, presumably in the force 
model, have been shown to exist that correlate with eclipsing of the sun. These effects are pronounced in 
ephemeris predictions and can be regarded as a limiting effect on the prediction accuracy. In spite of the existence 
of these unmodeled effects, the ability to obtain geodetic results that show an accuracy of 1-2 parts in 10 8 or better 
has not been diminished. The possibility that such unmodeled effects are constraining the ability to reach a few 
parts in 10 9 is under investigation. 
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Fig. 2. Comparison of Orbits: Week 420 estimated solution vs. Week 419 predicted solution for SV06 


The effects of eclipsing have been recently examined in two areas. First, the proper accounting of eclipsing in the 

" h T°L° f ? e diffCrCntial CqUati0nS ° f “ r ^ uircs care' wfih a fixed £,£%£ wht 
ossing the shadow boundary. A detailed study of this effect by Feulner /8/ showed that, while a significant effect 
was observed, it did not account for the observed effects in Fig. 1. significant effect 

m^els C of d a S GPS f ^ ** T** differen,ial ^ uations {ot conduction and radiation for simplified 
models of a GPS-type solar panel simultaneously with the orbit When the modeled satellite passed into shadow 

“ d *! r ! >e8inS 10 COOl> ^ SOlar paneIs Were found ^ « distinctly different tempentoSS? rise to a 
eimal imbalance or a net force due to the radiation properties. The study did show that this effect could have a 
significant (meter or tens or meter) effect on the orbit, however, additional work is £ pro“ 

^sf SSTSTm?? ° f d0UWC Ph3Se ** f ° r d * C3SeS — typically about 3 cm 

, .uT? ^ “ meaSUrC ° f orb,t accuracy, it cannot be implied that the orbit accuracy is 

formation of die double dif/* 1 * 1 ** RMS ' In f3Ct ’ bccause of cancellations of parameters other than clocks in the 
GEODETIC RESULTS 

^ t d ^c^v d S«WtfrST^ Uan ^ CS * SUCh a V eCt ° r baselines ’ « an measure of the sensitivity to 

.. , y ’ Kts d mlafiye position vectors for the sites estimated in Table 1 are given in Table 2 in terms of 

SLR or VLBlTho ws^th^ mdication °£. thc P rccision of estimates. Comparison of the average values with 
Suit W^K Tv! fences given in Table 3. It is particularly important to note that the (Lala/Tromso 
een JT th ° b r ed 545016 VLBI measured baseline was available, thus constituting a "blind test" of GPS In 
comparison with independently determined results from SLR and VLBI provide a measure of the 

parts iSlO 8 ^ bluer. ^ ^ ***“ tCchniques ' From Acsc comparisons, the GPS results have an accuracy of 1-2 


TABLE 2 Week-to-Week Estimate Repeatability, RMS about Mean (m) Over Three Weeks 


Baseline 

dx 

dy 

dz 

dl 

Tidbinbilla/American Samoa 
Tidbinbilla/Black Birch 
Onsala/Tromso 
Kokee Park/Ft. Davis 
Mojave/Owens Valley 

0.044 

0.005 

0.023 

0.024 

0.012 

0.144 

0.072 

0.029 

0.047 

0.033 

0.101 

0.020 

0.026 

0.022 

0.026 

0.177 

0.066 

0.011 

0.048 

0.007 

— 3 . Comparison of Estimates with SLR or VLBI (m) 

Baseline 

Longitude 

Latitude 

Height 

Length 

Mojave/Owens Valley 
Kokee Park/Ft. Davis 
Onsala/Tromso 

-0.014 

-0.013 

-0.019 

-0.007 

-0.008 

-0.025 

0.004 

0.017 

0.030 

0.000 

0.012 

0.026 
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Chapter 1 
Introduction 

The modeling of forces for the numerical integration of near Earth satellites 
is an important aspect of the orbit determination procedure. In applications using 
satellites such as Lageos, Topex, or the Global Positioning System (GPS), orbits 
must be known very precisely to obtain the most benefit from each satellite. To 
recover the orbits, many small perturbing forces must be modeled correctly 
within the numerical integration portion of the orbit determination procedure. 
This study is an investigation of one of the important perturbing forces, the solar 
radiation pressure (SRP) and its modeling when a fixed-mesh/fixed-order class II 
general formulation multistep numerical integrator is used to propagate the orbit. 

There are two important aspects of the SRP force and its modeling: (1) the 
shadow model and (2) the nature of the numerical integration algorithm in which 
the SRP force is to be modeled. For most applications, the shadow of the Earth 
can be modeled as a cylinder extending behind the Earth, with respect to the Sun, 
enclosing the region of total shadowing, the umbra. In this model, there is no 
region of partial shadowing (penumbra), thus a true step function in the force 
model is caused by the SRP when satellites pass into and out of the Earth’s sha- 
dow. For more precise applications, the shadow can be modeled as two cones, 
one extending inward behind the Earth enclosing the umbra and the other 
extending outwards enclosing the penumbra. The partial shadow region, the 
penumbra, provides a smooth transition from total sunlight to total shadow. 
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Multistep numerical integrators are frequently used in the orbit determina- 
tion procedure because of their efficiency relative to single step methods [Lund- 
berg, 1981]. A general formulation multistep integrator is derived assuming the 
forcing function / of the problem 

-^=y (2) =/(r,y,y (1) ) (1-1) 

dt l 

is smooth and continuous so that the integrals of / can be approximated by 
integrating an m-1 order polynomial formed using m nodes, where the time step 
between each node is referred to as the step size . If the second integral of the 
polynomial is formed explicitly, as opposed to reducing the second order dif- 
ferential equations into twice as many first order equations, the integrator is 
referred to as a class II integrator. When the the conical model is used in the 
numerical integration, the SRP force appears smooth and continuous as long as 
the step size is small enough to sample the penumbra more than a few times. 
However, the force model may have an apparent discontinuity which violates the 
assumption of a smooth, continuous function of the numerical integrator. This 
discontinuity may occur if (1) the step size is longer than the time satellites 
spend in the penumbra, meaning the numerical integrator can step at most once 
in the penumbra, or (2) a cylindrical model is used. Thus, the polynomial does 
not accurately represent the forcing function and the integrator does not step 
correctly across the shadow boundary. GPS satellites which .pass through the 
shadow take on the order of 60-100 seconds to cross the penumbra. However, to 
reduce the computational load, it is desired to use a 500 second numerical 
integration step size in the orbit determination process. Consequently, GPS 
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satellite orbits are good candidates for studying the effect that this problem has 
on the numerical solution for the orbits and techniques that can be used to allevi- 
ate this problem. 

One possible solution to this problem is to reduce the step size used by the 
numerical integrator, but as was mentioned previously, it is desired to have a 
larger step size to reduce the computational load. Another possible solution is to 
integrate up to the boundary of the cylindrical model and restart the integrator at 
the boundary using the new level of the SRP force. This solution increases the 
computational load to a lesser degree than decreasing the step size, but it also 
introduces errors in the restarting process and is not an appropriate alternative 
when a conical shadow model is used. 

The modified back differences (MBD) method, proposed by Hubbard, and 
described by Anderle (1973), corrects for most of the error by introducing 
modifications to the state at the leading node and modifying the back differences 
of the forcing function so the polynomial approximation is again valid. Lund- 
berg (1985) also presents some results of this method. Modifications to the 
MBD method will also be examined including a new approach to the leading 
node state modifications. 

Two types of tests were performed to analyze the MBD methods: (1) 
Numerical integration tests were performed using different step sizes, shadow 
models and numerical precision to analyze the direct effect of the MBD method 
on the propagation of the state given true initial conditions; (2) Batch filter esti- 
mations were performed to analyze the effects of the MBD method on the pred- 
iction of GPS orbits. 
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Chapter 2 discusses the numerical integrator, the orbit determination pro- 
cedure, and the two applications to be used in the evaluation of the MBD 
method. Chapter 3 describes the Global Positioning System, including the orbits 
of the satellites and some of the data types available for orbit determination. 
Chapter 4 discusses the MBD method proposed by Hubbard as well as proposed 
modifications to this approach. Chapter 5 presents results from integration tests 
and batch filter estimations using both simulated and real data. Finally, Chapter 
6 contains a discussion of the major conclusions and some suggestions for future 
work. 



Chapter 6 
Conclusions 


The purpose of this work was to investigate solutions to the problem of 
integrating across the step function-like solar radiation pressure boundaries when 
large integration step sizes are used. Two shadow models were examined: the 
cylindrical model, which has true step function boundaries between full sunlight 
and full shadow, and the conical model, which has a smooth transition from full 
sunlight to full shadow only when the penumbra is sampled at least a few times 
per crossing. In the orbit determination of GPS satellites, it is desired to use an 
integration step size of at least 5 times as long as it takes the satellite to cross the 
penumbra, so the conical shadow model will appear in some cases to be a step 
function. Only when an integration node happens to lie near the middle of the 
penumbra will the step function nature be lessened, although the state at the next 
node will still not be correct. 

Hubbard’s original MBD method was the basis for most of the algorithms 
used to alleviate this problem at large step sizes. The following were proposed 
modifications to the original MBD method: the extension of the method to the 
conical shadow model, including two different approximations to the fractional 
acceleration in the penumbra; the interpolation method for leading node state 
corrections; and corrections due to the point mass term. The other method exam- 
ined to alleviate this problem was restarting the integrator at each cylindrical 
shadow boundary crossing. 

In Sections 5.1.1 and 5.1.2, it was found that Hubbard’s MBD method with 
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integral leading node state corrections is effective in reducing the error. Further- 
more, the proposed point mass term corrections at the leading node are necessary 
to obtain an accurate solution. Also from Sections 5.1.1 and 5.1.2, it was found 
that the MBD method with the proposed interpolated leading node state correc- 
tions is the most effective algorithm in reducing the error. Again, the point mass 
term corrections were found necessary for an accurate solution but the choice of 
penumbra approximations does not affect the results. Section 5.1.3 showed that 
restarting the integrator at each cylindrical shadow model boundary was very 
effective in reducing the error, but from Table 5-1, significant additional com- 
puter time was required to restart the integrator. Also from Section 5.1.3, it was 
found that the cylindrical shadow model, when the MBD method with interpo- 
lated leading node state corrections was applied, ' was an accurate approximation 
to the conical model. 

Section 5.2 showed how the magnitude of the errors incurred by estimating 
trajectories with no special action taken at the shadow boundaries could be 
reduced by using the MBD method with interpolated leading node state correc- 
tions applied to either shadow model or to a lesser degree, the MBD method with 
integral leading node state corrections applied to the cylindrical model. Table 
5-5 showed that the character of crossing the shadow boundaries with no special 
action has RMS values of less than 1 m in the radial and normal directions and 
up to 20 m in the transverse direction. 

Section 5.3 applied the MBD method to actual observations to reduce the 
prediction errors. Table 5-7 showed that the RMS values in the radial and nor- 
mal directions were as large as 4 m indicating that there is still mismodeled or 
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unmodeled force(s) in the GPS force model causing errors at least as large as that 
of crossing the shadow boundaries with no special action. Possible error sources 
include: (1) modeling the Y-bias acceleration as a constant for satellites which 
passed through the shadow, including the satellites which were shut down during 
periods of eclipse, (2) Earth radiation or albedo pressure, (3) the use of a constant 
area model instead of a variable area model or (4) inclusion of the Moon’s sha- 
dowing effect. 

Additional unanswered questions or possible improvements have arisen in 
the course of this study. The first one would be modeling the Earth as a nons- 
pherical body. The effective radius of the Earth used in this study, 6371.0 km , is 
within about 7 km of the true radius of the Earth everywhere, but this can result 
in crossing times which are off by as much as 2-3 seconds. The implications of 
this was not examined. Second, the shadowing effect of the Moon’s shadow 
should be included. Although MSODP has an option to include the Moon’s sha- 
dow in the variable area model, the constant area model does not. At the Earth’s 
distance, the Moon’s umbra has almost vanished, so the MBD method should not 
have to be used. Finally, it might be worth while to actively check for 
sunlight/penumbra/sunlight crossings which might be skipped in the numerical 
integration. If the value of 0+a 5 +a e is saved at the past three integration nodes, 
a second order polynomial can be formed using these values. For example, if it 
is assumed the satellite is in sunlight and a boundary crossing has not been 
detected, the penumbra may have been crossed if the minimum value of the 
polynomial is both less than 0 and between the nodes of the polynomial. This 
can also be applied to the cylindrical model by monitoring the values of rsina 
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and cosa. 

A further application of this method could be to model impulse maneuvers. 
The trajectory would be integrated past the impulse by one step and interpolation 
would be used to find the state immediately before the maneuver. The AV would 
be added directly to the velocity at the time of the impulse and the back differ- 
ences should be modified to reflect the maneuver. Finally, interpolation would 
be used to find the new state at the leading node. 
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IMPROVING THE NUMERICAL INTEGRATION SOLUTION OF 
SATELLITE ORBITS IN THE PRESENCE OF SOLAR RADIATION PRESSURE 
USING MODIFIED BACK DIFFERENCES ^ 

J. B. Lundberg,* M. R. Feulner,** P. A. M. Abusali, + and C. S. Ho^ 


Numerical integration techniques are frequently used to form the 
solution of the equations that describe the orbital motion of a satellite. 
The basic errors associated with numerical integration can be described 
in terms of round-off and truncation errors resulting from a finite number 
of significant digits used in the computations and the functional 
approximation used in deriving the numerical integration algorithm. 
Additional numerical integration errors can result when the time scale of 
some of the perturbing forces is smaller than or on the same order as the 
integration step size. This type of error can occur when solar radiation 
pressure is included in the force model and the trajectory crosses through 
the shadow of a planet. The method of modified back differences is a 
technique that significantly reduces the numerical integration errors 
associated with crossing shadow boundaries using a fixed— mesh 
multistep integrator without a significant increase in computer run time. 


INTRODUCTION 

Many near Earth satellite missions require the most precise trajectory solutions 
possible from the available tracking information and solution techniques with some orbit 
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requirements ranging from a few meters to a few decimeters or less in position. Precise 

“T 5' P ‘ C!d ! 5 ! Involve forming estimates of the orbit state parameters as well as 

nS* 1 ? emP, ™ aI P f™ eters used ,0 model 1* forces acting on the satellite. 
Any numencal errors introduced into the orbit determination process can affect the 

toends „l Se ,h eSUma,ed P r ameterS ' Si " Ce ** predsion orbit determination precess 
depends upon the numencal integration of ihe equations of motion, it is important to 

SS effiCien ' imegra,i0n "" is —e with 

the n,,men>,T m ? riCal iMegrati f erTors m unav ° ida bln, the issue of how to carry out 

compurer’ S.’ST" •“ r d "“ d 10 Se,eC ' in8 “ ***** *" which *> -Sting 

EmrsateZ L?! T Ca OTOTS 3,6 ,0lerable - Thc C ‘) ua,i0ns of "-Orion for a near 
cann satellite can be written as 


d 2 ^ 


dt 


2 ^ ^ r, r ) — f(t, r, r ) 


( 1 ) 


?" d ” vera ! jmpottant characteristics about the numerical solution of Eq. (1) have 

moTeffimtuhan ' ?n ' I™' 3 CIa “ “ imegratOT which s0lves ,his ***** directly is 
more efficient than a Class ! integrator which would require writing Eq (1) as a system 

single' sten r Second ’ multistep integrators are more efficient than 

ccltncU „n ”!r*u Ut,a) a g T' hmS in solvi "« ritis system. Finally, when th“ 
i7 “2 ° f ° rb “ ,S sma11 ' ,he K > uano " s of "“ion are sufficiently smooth that step 

'° maintain ‘ he SOlU,i ° n Consequently, a Class £ 

fixed mesh, multistep integrator of either the summed form (Second Sum 

orbi^detemdnation ^mf *” (AdamS ' ^ « wide * ia ” a ^ P**™' 

Eo (l!h Z!li, S iI eP 1 , mesralor is based -PO" -h« assumpdon that the right hand side of 
SviomL T a COmmU ° US a " d Can ** a PP rox '—ated m sufficient accuracy by a 

SongTe m l Zt US Z a “ ,° f diS a CTete values ° f f P°"-P" Kd a ‘ *• --^adonUes 

force this a i T' 20 S ,° lar radiatlon P ressure < SRP ) is included as a perturbing 

satellite nasS fT T 7 When 3 cylindrical shadow model is used and thf 

a con ca shaHo T ?• I W ° the Umbra * ™ s assumption is also violated when 

step size is larg^ 7?^ if JJ Sed ^ n combination with a fixed-mesh algorithm in which the 

revolution Thf -°^ ■ * ^ tbe P enumbra is sampled only once or twice at most per 

len^ ofonte f atl0I f C3n le3d t0 3ppreciable numerical i^tion errors for arc 
lengths of only a few tens of revolutions. 

reduced neL?A°^ fiCati ° nS t0 the fixed ~ mesh > multistep integrator can be employed to 
method of When crossin S shadow boundaries. However, the 

shadow l?ef dl nt rf CrenCCS ?“ ^ ^ ^ ^ thC Cylbdrical 0r Co;ical 
size and can he r a m‘ h a ° y alteratl0n of the multistep integrator’s mesh or step 

thL dTscussion ofT “ T ,y a SHght inCreaSC in COm P uter To begin 

propo tS Z ult* mC ?°? ° f m0dified back d ^rences, the original method, as 
proposed by Hubbard and described by Anderle, 1 is discussed. Improvements to 
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Hubbard s approach for the case involving a cylindrical shadow model will be offered. 
The extension of the modified back difference approach to cases using a conical shadow 
model will also be demonstrated. 

The integrator to be used in this investigation is a general formulation (Adams’) 
Class II integrator of the form 

yn+s = y n + (sh) y„ + (sh) 2 £ a jtS V-j -1 f n (2) 

j=l 

yim = y n + (sh)£Pj.sV j - 1 f n 
j=» 


where h is the step size, are the back differences of the force, f(t, y, y), formed at 
|n, a j.s and Pj.s are the coefficients used in the prediction (s=l), correction (s=0), or 
interpolation (s<0) stages of the integrator. Given a value for s, the coefficients are 
computed using a recursive algorithm. 2 ’^ The selection of the general formulation 
integrator instead of the summed formulation was arbitrary; the method of modified back 
differences is applicable to'either formulation. 


HUBBARD’S METHOD 

The method of modified back differences as proposed by Hubbard is an attempt to 
reduce the numerical errors associated with the crossing of the cylindrical shadow model 
boundary between sunlight and shadow which introduces a step function change in the 
forcing function. In the cylindrical model, a satellite is assumed to be in shadow when 
cos e < 0 and (r sin e) < R e are satisfied where e is the angle between the Earth-satellite 
and the Earth-sun lines, r is the radial distance of the satellite from the origin, and R e is 
the effective radius of the Earth. The method of modified back differences consists of 
three steps. 

Step 1 

To begin the application of the method of modified back differences, it is necessary 
for the numerical integrator to detect the crossing of a shadow boundary and to solve for 
the time and state at the boundary crossing. Although a shadow boundary is crossed, the 
integration should be propagated across the boundary as if there were no boundary. This 
will assure that the integrated solution up to the time of the shadow boundary crossing is 
correct. After the crossing, the trajectory will have errors that are proportional to the 
difference between the actual SRP model and the value used to form the integrated 
solution. If the last three integration nodes are taken to be tj, to and t_i with tc, the time 
of crossing, in the interval [to.tj ], the shadow boundary can be found by using a quadratic 
approximation of the function T(t) = r(t) sina(t) where T(t c ) = R c . By evaluating the 



4 


function T(t) at the nodes, the shadow boundary crossing is found by evaluating the 
expression ' & 


8t = 


2V 2 T! 


2Vt i + V 2 T, - p £ (2VT, + V 2 T,) 2 - 8V 2 T! (Tj - RjJ 


1/2 


(3) 


where V Tl = Tj - T 0 , V 2 T, = T, - 2T 0 + T_! , 0 = 1 for VT, > 0, (3 = -1 for VT, < 0, 
and from this tc = tj - St. The accuracy of ^ is dependent on the spacing between the 
nodes, i.e. the integration step size, because of the the second order polynomial 
approximation. If a more accurate solution is desired, t<, is taken to be the last solution of 
the boundary crossing time, the range [t_! ] is reduced, the trajectory is interpolated at 

t_i , ^ and tj to evaluate T(t), and St is recomputed using Eq. (3). This algorithm can be 
iterated until the value of ^ converges. 

Step 2 

At this point, the back differences, V ! f n , represent the forcing function evaluated at 
thelevel of SRP before the boundary crossing. To avoid the step function nature of the 
cylindrical shadow model boundary, it is necessary to modify the back differences to be 
consistent with the new level of SRP at the leading node, t n ._ To accomplish this, the 
back differences of the SRP are formed at the leading node, V j f SRPn , where the SRP has 
been evaluated at each node as if the satellite were in full suniight. The SRP back 
differences are added to (subtracted from) the back differences of the total forcing 
function when the satellite is passing out of (into) the shadow. 

Step 3 

Since the state at the leading node (r n , r n ) is in error, the integrated effect of the 
SRP from ^ to t n must be used to modify the state at the leading node. If the SRP is 
assumed to be constant over one integration step, corrections of 


2^ tn ~ t c) 2 fsRP.n 
Sr = y (t n -tc) fsRp, n 


must be applied to the position and velocity at the leading node, respectively, where 

fsRP.n is the SRP evaluated at the leading node and y is +1 (-1) if the satellite is moving 
out of (into) shadow. 


r e m 1 e< ^ differences algorithm is represented in Fig. 1 where the solid 
me represents t e actual trajectory, the dashed line represents the trajectory if there were 
no oun ary, and the dotted line represents the new "past" trajectory after the method of 
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modified back differences is applied. The technique of using an approximate integral of 
the SRP to correct the state at the leading node will be referred to as the integral method. 


True Trajectory 

Predicted Trajectory Before Correction 
Past Trajectory After Correction 



a SRP 


0-L 


*«— 1 


Fig. 1 Cylindrical SRP Force Model with Trajectory 
Representation of the Method of Modified Back Differences 


MODIFICATIONS TO HUBBARD’S METHOD 

Two important modifications to Hubbard’s method have been identified. The first 
comes as a result of the leading node state correction. Once the leading node state 
correction is made, the force at the leading node, f n , is no longer consistent with the state. 
Since the leading node state correction is small, the only important force to re-evaluate is 
the point mass term (PMT) of the geopotential. Using the corrected state, the PMT is 
re-evaluated and is used to update f n . However, since f„ has been altered, the state 
should be corrected for this small change in the force function. If the difference in the 
PMT is assumed linear with a value of zero at the shadow boundary crossing to the value 
of the computed difference at the leading node, the force can be integrated to yield the 
leading node state corrections 
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^ ~ "g" (^“tc) 2 SfpMT,n 
5F = -J ( l n~tc) SfpMjn 


Where Sf P ^ n is the difference in the PMT evaluated at the leading node. Since another 
leading node state correction has been made, the PMT should be re-evaluated once more. 

econ re-eva uation can then lead to another leading node state correction, so this 
process can be iterated more than once. 

The second important modification to Hubbard’s method is how the leading node 

inte^TH^T 5 arC made ' SmCe the integTal corr ections are only valid over short 
ntervals, an . f tematIve method is needed to increase the accuracy at larger integration 

h 2L ! S ?ff ' AftCr “f t tmg t0 find thC b0undary crossing time and state and correctSg the 

crassinfr^ f r. fU " C,i0n fOT ,he Cha "«' SRP ' thc * <■» bound^ 

? T d ln * 1S ,nter POlation equations of numerical integrator to solve for 
the state at the leadmg node. Selecting Eq. (2) can be written as 


m 


yn=yc-(t n -t c )£p J , s V>- 1 f n (4) 

j=l 

y n = y c ~ (t n -tc) y n - (t n -tc) 2 £ Oj.sV 1 - 1 ^ 

j=l 

where y c and y c are the position and velocity at t<., respectively, and s = (tc - t n )/h. This 

C TT ng - the State at 1116 leadin S node wil1 be called the interpolation 

™ ° . aPP i 7 the P° int mass term corrections to the interpolation method, it is only 

]p , ary !° updat ^ the back differences with the new PMT, then solve for the new 
leading node using Eq. (4). 


APPLICATION TO THE CONICAL SHADOW MODEL 


cylindrical' S ^ dOW T 1 ' 1 iS m0rc rcpresentative of the Earth's shadow than the 

exurntL udn^th In “ 8 ' “ parame,ers used 10 d =«tribe the conical model can be 

expressed using the following equations 


cx s = tan 1 



ctg — tan 



0 = cos 1 


r e * r s 

|r e | |r s | 


where R s is the effective 


radius of the sun, r s is the vector from the satellite to the sun, R e 
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is the effective radius of the Earth, and 7 e is the vector from the satellite to the Earth. 
The satellite is assumed to be in the region of partial shadowing, the penumbra, when 
(a c +a s ) > 8 > (cc,.— ctj) and in the full shadow, the umbra, when 8 < (o^— a s ). 

To implement the method of modified back differences with the conical model, the 
time and state at both the sunlight/penumbra and umbra/penumbra boundaries must be 
found. By defining the function T(t) as T(t) = 8-Oe-ya s , where y is +1 for the 
sunlight/penumbra crossing and -1 for the umbra/penumbra crossing, then the boundary 
crossings satisfy the condition that T(t c ) = 0. The times of the boundary crossings can 
then be found by writing Eq. (3) in the form 



where P - 1 for VTi > 0, (3 - — 1 for VTj <0, and tc = tj — St as before. As in the case of 
the cylindrical shadow model, the numerical integration is carried out assuming that the 
satellite is in full sunlight (shadow) until both of the shadow boundaries are crossed and 
the satellite is in full shadow (sunlight). 

At this point both shadow boundary crossing times are found and the corrections to 
the state at the leading node are applied. To extend Hubbard’s integral approach to the 
conical model, the fractional SRP encountered in the penumbra as well as the appropriate 
SRP from the second boundary crossing to the leading node must be modeled. If the 
SRP in the penumbra is modeled as a linear function with a value of 0 at the 
umbra/penumbra boundary to a value of fsRp, n at the penumbra/sunlight boundary, the 
leading node corrections can be written as 


8r = y 
5r = y 


\ (tn-t c2 ) 2 + J (tc2-tcl) 2 


(tn~tc2) + J (tc2~t c i) 


f SRP,n 
fsRP.n 


where t^ and t^ are the times at which the boundaries, umbra/penumbra or 
sunlight/penumbra, were crossed and y is +1 (-1) if the satellite is passing out of (into) 
shadow. 

If the interpolation method is used, an integral correction is applied to the state at 
the second boundary crossing, and this corrected state is used in the interpolation 
equations to solve for the leading node state. The integral corrections applied to the state 
at the second boundary crossing are 
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^ “ 5 y ) 2 fsRP.n 
S* = yY(t c 2-tcl) ^SRP.n 


of Vn? 0f A ? e pe L numbra 1S used where y is +1 (-1) if the satellite is passing out 
no ii a/pST’ After the mte f al COITection is used to update the state at the leading 
applied correctlons t0 the *aie and force at the leading node are iteratively 

70 invest >gate the accuracy of using a linear model for the shape of the solar 

thfn^° n h K Pre f U 7 u penumbra ’ a non hnear model was derived by assuming that in 

the neighborhood of the solar disk the Earth’s limb is linear and the satellite’s velocity is 
nearly constant. Under these conditions, the fractional SRP can be written as 


f SRP = 


_ cos-^l-gs) 


n 


- — (1 — 2s) (s-s 2 ) 3/2 


MdT- f l n ° ndime " si °" al Pieter such that s = 0 at the umbra/penumbra boundary 
and s-1 at the sunltght/peuumbra boundary. For this nonlinear model, the integral 
correction for the penumbra effect becomes ^ 


Sr - 32 Y ( l c 2 - tel ) 2 fsRP.n 

Sr = ~2 Y 0c2 - tel ) fsRP.n 


rZ™ er ’ di< ! n0 ‘ provide any si S nificam ™pr»vemeut over the linear model 

Consequently, only results using the linear model are given. 


TEST RESULTS 

aleorithm° Fir,? ;? f ^ Were performed to evaluate the modified back differences 
the he<;t • . as necessar y to use numerical integration comparisons to identify 

were used^^^chfilter ** d “ fa ~ Then ’ the best variations 

o rb U tCStS t0 dCtermine thC CffeCt ° f the al S orith ™ on the 

of those in the GPS 00^1^“ * * * ** investi g ation « representative 
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Algorithm Comparison 

The numerical integration comparisons were done in two groups, one using a 
cylindrical shadow model and one using a conical shadow model. The initial conditions 
were the same for all tests: 

semimajor axis = 26,000 km ascending node = 0° 

eccentricity = 0.01 argument of perigee = 0° 

inclination = 55.0° mean anomaly = 0° 

The time span of the tests was taken to be two weeks or 28 orbital revolutions. The 
dynamical force model included only a constant area SRP model and a geopotential 
model including terms through degree and order eight. For the SRP model reflectivity 
coefficient (rj) was 0.22, the effective radius of the Earth (Re) was taken to be 6371.0 km 
and the effective radius of the Sun (R s ) was taken to be 696000.0 km. The Sun was 
assumed to be in the equatorial plane with the Earth— Sun direction initially in line with 
the line of nodes of the satellite. The number of function evaluations represented in the 
back difference table (m in Eq.s 2 and 4) was 14. These test were performed on a dual 
Cyber 170/750 computer utilizing a 60 bit word with a 48 bit mantissa. 

To evaluate the various permutations of the modified back difference algorithm, 
reference trajectories were computed for each shadow model. The reference trajectory 
for the cylindrical shadow model comparisons was computed using a 500 second 
integration step size which was restarted at each cylindrical model shadow boundary. 
This was chosen over a double precision integration using a 5 second step size because 
even at small step sizes, the solution is flawed by the step function introduced in the force 
model at the cylindrical shadow boundary. The reference for the conical shadow model 
comparisons was computed using a 5 second step size and double precision (28 
significant digits) arithmetic. This integration does not suffer from a step function in the 
force model because it samples the penumbra more than 12 times per pass. 

The transverse (alongtrack) errors for the selected algorithms involving both the 
cylindrical and conical shadow models are shown in Figs. 2 and 3. Only the transverse 
error is displayed because it is known to be much larger than the errors in the normal 
(crosstrack) or radial directions. For each shadow model, the interpolation method with 
two iterations of the PMT correction results in the best ephemeris with a maximum error 
of less than 3 cm when applied to the conical model and a maximum error of less than 1 
cm when applied to the cylindrical model. It is also shown that the PMT corrections are 
important in obtaining high accuracy solutions, although it was found that no more than 
two PMT corrections are necessary. In comparison, the cylindrical model with no 
special action taken at the boundaries had a maximum transverse error of 7 m, and the 

conical model with no special action at the boundaries had a maximum transverse error 
of 20 m. 
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□Integral MBD 
with 2 PMT corrections 

+ Interpolation MBD 
with 2 PMT corrections 


* Integral MBD 
a Interpolation MBD 



Time (days) 


Fig. 3 Results of the MBD Algorithm 
Using the Conical Shadow Model 


□Integral MBD 
with 2 PMT corrections 


+ Interpolation MBD 
with 2 PMT corrections 


* Integral MBD 


A Interpolation MBD 
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The method of modified back differences will be an effective tool only if it can be 
implemented at a reasonable cost in terms of computer time. To evaluate the additional 
cost of using the method of modified back differences, trajectories spanning one week 
(14 orbital revolutions) were computed in which no special action was taken at the 
boundaries for both the conical and cylindrical models, and the execution time of these 
trajectories were used as the baseline for evaluating the method of modified back 
differences. As in the previous tests, the force model included only SRP and a 
geopotential model complete to degree and order eight An integrator step size of 500 
seconds was used for all cases. It was found that the method of modified back 
differences increased the execution time by about 10% for the cylindrical case and about 
20% for the conical case over the baseline cases; whereas restarting the integrator at the 
boundary of the cylindrical model represented an increase of over 100%. The relative 
cost of the method of modified back differences will decrease as the complexity of the 
force model increases since its cost is proportional to the number of boundary crossings. 
The relative cost of restarting the integrator at the boundaries will remain high since it is 
proportional to the cost of evaluating the forcing function as well as the number of 
shadow crossings. 

Batch Filter Tests 

To measure the effect that this numerical integration error has on batch filter 
solutions and the extrapolation of these solutions, tests were conducted in which the only 
source of error is the result of using a relative large step size to cross the shadow 
boundary. The scenario for the batch filter tests involved the satellites PRN 3, 6, 8, 9, 1 1, 
12, and 13 of the GPS constellation during Weeks 419 and 420 which correpsonds to the 
CASE UNO experiment. 4 PRN 3, 8, 11 and 13 are in one orbital plane which eclipses 
during these weeks, and PRN 6, 9 and 12 are in another orbital plane that is in full 
sunlight. The locations of the twelve ground stations involved in the experiment are 
shown in Table 1. The observations used in this experiment were double differenced 
carrier phase measurements computed along simulated "true" trajectories for the seven 
satellites using the time tags of the actual observations made during this two week 
period. Simulated "true" trajectories and noiseless observations were used 1) to remove 
the effect of other error sources on the batch filter solutions and 2) to provide a true 
solution to evaluate the technique. The time tags of the actual observations were used to 
provide a spatial and temporal data distribution representative of actual tracking. There 
were over 80,000 observations which were simulated along the true trajectory that was 
generated using a 20 second integration step size with a conical shadow model. The step 
size was chosen to have at least three steps in the penumbra per pass to avoid introducing 
the step function change in the SRP at the shadow boundary which occurs with larger 
step sizes. 

The force model in all cases included a constant area SRP model and a 
geopotential model complete to degree and order eight. The estimated parameters in 
these tests included the initial conditions for each satellite in the constellation and the 
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Table I 


GPS GROUND STATIONS 


Station 

Longitude (°E) 

Latitude (°N) 

American Samoa 
Black Birch, New Zeland 
Onsala, Sweden 
Tidbinbilla, Austrailia 
Tromso, Norway 
Wettzell, FRG 
Yellowknife, Canada 

Ft. Davis, TX 
Kokee Park, HI 
Mojave, CA 
Owens Valley, CA 
Westford, MA 

-170.7147 

173.8056 

11.9255 

148.9786 

18.9383 

12.8786 

-114.4689 

-103.9470 

-159.6649 

-116.8882 

-118.2836 

-71.4933 

-14.3322 

-41.7452 

57.3953 

-35.3976 

69.6627 

49.1446 

62.4769 

30.6360 

22.1263 

35.3315 

37.2320 

42.6133 


Ihe back difference Taui (mT wm 14 '’S'co'^ nUmber of nodes "presented in 
X-MP/24 using singie arithmetic (64 bit w^dST/g ££££3“ 0 " * °* 

applied^ ^ 0> "if, “ 

cylindrical shadow model »nH ™ mt<;r Polation method as applied to the 

shadow model. Each of these algorithm lnt ^ r PP ad ° n m ethod as applied to a conical 
correction. For comDarison the f r a • vanatl0ns incI uded two iterations of the PMT 
action at the boundaries wem ^ and conicaI shadow models with no special 
estimating the trajectories using 1 on 1v th L the tests * The batch ft 1161 " tests involved 
estimated trajectories were P vtl , y *? C data fr0m the ^ week of the study. The 
reference trajectory. The results fromToth^h * if Week ^ com P ared with the 

.0 evaluate the effectiveness of the modl^Sj^J^ST^ ” 

of PRN 3 which'fe Slipslnfi^d PRN^wh'h' n™™ for ‘ hC es,imatcd >rajectories 
results were obtained Esing the inte ,S “l/ 1 ? S ? nllght during week 419. These 
conical shadow model to d»T T ™ ethod w,th PMT corrections and a 
orbit errors for PRN 3 and PRN 6 under th^’ SinuIa [! ) !’ Flg * 5 gives m example of the 
is taken at the boundaries These result*; H condlt ‘ ons exce Pt that no special action 
differences can be used to reduce^ 7h Iff iemonstrate dial the method of modified back 

batch fiber solutions by two orders of mag^ fmm " ^Z‘°" °" 

estimated and trae°toj rcton^'fo ° f ‘ he difftrenc “ between the 

trajectories for the week 419 estimates. All of the MBD algorithms 
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Fig. 4 Transverse Orbit Errors for the Batch Filter 
Solutions of Week 419 Using the Inteipolation MBD Method 
with a Conical Shadow Model 



PRN 6 
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Fig. 5 Transverse Orbit Errors for the Batch Filter Solutions 
of Week 419 When No Special Action is Taken at the Shadow 
Boundaries of a Conical Model 
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Ih« double difference measurement retdduTtor toe i"i. I?o 3 . shows the RMS of 
MBD algorithms result in RMS vl« of tss L I Zl n A " of *« 

special action was taken at the hnunHanV u mm e l ^ e cases in which no 

suggest that the MBD method should r»m ’d^ \ UCS greater 7 cm - These tables 
week 420 than the model s in which no special Irionk ° f ““ trajeCt0r> ' int0 

Table 2 

rc ^ the differences between the 

ESTIMATED AND TRUE TRAJECTORIES OF WEEK 419 

PRN 8 PRN 9 PRN 1 1 PRN12 PRN 13 


Model 
Cylindrical 


Component PRN 3 


No Special 
Action 
Cylindrical 
Integral 
MBD 
Cylindrical 
Interpolation 
MBD 
Conical 
No Special 
Action 
Conical 
Interpolation 
MBD 


Radial 

Trans 

Normal 

Radial 

Trans 

Normal 

Radial 

Trans 

Normal 

Radial 

Trans 

Normal 

Radial 

Trans 

Normal 


0.176 

1.201 

0.100 

0.001 

0.003 

0.003 

0.010 

0.004 

0.003 

0.154 

0.652 

0.112 

0.010 

0.004 

0.003 


PRN 6 
0.081 


0.460 

0.027 

0.002 

0.003 

0.002 

0.002 

0.002 

0.002 

0.090 

0.217 

0.090 

0.002 

0.002 

0.002 


0.305 

2.626 

0.020 

0.003 

0.007 

0.004 

0.002 

0.003 

0.002 

0.183 

1.109 

0.099 

0.002 

0.003 

0.002 


0.065 

0.433 

0.080 

0.002 

0.003 

0.002 

0.001 

0.002 

0.002 

0.099 

0.267 

0.111 

0.002 

0.002 

0.002 


0.174 

0.811 

0.085 

0.002 

0.004 

0.003 

0.003 

0.004 

0.003 

0.151 

0.671 

0.129 

0.003 

0.005 

0.003 


0.067 

0.300 

0.111 

" 0.002 

0.003 

0.002 

0.002 

0.002 

0.002 

0.008 

0.191 

0.121 

0.002 

0.002 

0.002 


0.236 

0.929 

0.104 

0.003 

0.005 

0.003 

0.002 

0.003 

0.002 

0.171 

0.623 

0.130 

0.002 

0.003 

0.002 
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without the MBD approach. Also, the interpolation method applied to either shadow 
model has RMS trajectory differences of less than 1 cm while the cylindrical shadow 
model with integral corrections has an RMS less than 9 cm. However, the orbit 
predictions in which no special action is taken at the boundaries have RMS values as 
high as 30 m for the cylindrical model and as high as 9 m for the conical model. 





Table 4 





RMS (m) OF THE DIFFERENCES BETWEEN THE 
EXTRAPOLATED AND TRUE TRAJECTORIES OF WEEK 420 


Model 

Component 

PRN 3 

PRN 6 

PRN 8 

PRN 9 

PRN 11 

PRN 12 

PRN 13 

Cylindrical 
No Special 
Action 

Radial 

Trans 

Normal 

0.337 

17.316 

0.098 

0.185 

1.073 

0.027 

0.939 

31.361 

0.018 

0.070 

0.904 

0.079 

0.757 

7.759 

0.095 

0.182 

0.838 

0.111 

0.517 

3.244 

0.107 

Cylindrical 

Integral 

MBD 

Radial 

Trans 

Normal 

0.006 

0.020 

0.004 

0.003 

0.005 

0.003 

0.005 

0.085 

0.004 

0.003 

0.006 

0.003 

0.004 

0.047 

0.004 

0.003 

0.005 

0.003 

0.004 

0.051 

0.004 

Cylindrical 

Interpolation 

MBD 

Radial 

Trans 

Normal 

0.004 
' 0.016 
0.004 

0.002 

0.003 

0.002 

0.004 

0.007 

0.004 

0.002 

0.003 

0.002 

0.002 

0.003 

0.002 

nu 

6.003 

0.009 

0.004 

Conical 
No Special 
Action 

Radial 

Trans 

Normal 

0.285 

1.700 

0.113 

0.088 

0.420 

0.090 

0.459 

9.471 

0.099 

0.078 

0.290 

0.111 

m 

H 


Conical 

Interpolation 

MBD 

Radial 

Trans 

Normal 

0.005 

0.009 

0.004 

0.003 

0.004 

0.003 

0.003 

0.006 

0.003 

0.002 

0.003 

0.002 

0.004 

0.008 

0.004 

0.003 

0.005 

0.003 

6.004 

0.013 

0.004 


CONCLUSIONS 

It has been demonstrated that the method of modified back differences applied to 
either the cylindrical or conical shadow model significantly reduces the numerical 
integration errors associated with the crossing of shadow boundaries. While Hubbard’s 
integral approach can produce significant improvements to the trajectory solution, the 
interpolation method provides the best overall results. Also, it has been shown’ that 
iterating on the PMT correction is also important for achieving the best overall results. 
In addition, it has been demonstrated that the method of modified back differences can be 
implemented with only a small increase in execution time. 
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GPS OR aJTm A 1 ?? BAseLESTE: EXPERIMENTS: 
MINI-MAC/TI COMPARISONS 


B.E. Schutz, P.A.M. Abusali, C.S. Ho, B.D. Tapley 
, Center for Space Research 
University of Texas at Austin, Austin, Texas 78712 USA 


ABSTRACT 


Major changes in the Cooperauve International GPS Network (CIGNET) tracking hardware took place in 

«?• JJ* 1° 1989 ’ the CIGNET so,e, y consisted of TI-4100 receivers except for the Mini-Mac receiver 
at Tsukuba. In January, 1989, the National Geodetic Survey (NGS) installed Mini-Mac 2816AT receivers 
at Mojave, Westford and Richmond which, with the Tsukuba receiver, brought the network to an almost 
equal number of 11 and Mini-Mac receivers. Simultaneous operations of TI and Mini-Mac receivers were 
conducted by the NGS at the three U.S. sites for a period of several weeks. Using this simultaneous data 
set m addition to the other CIGNET data, an investigation has been conducted into the compatibility of the 
- and the mixed receiver configuration. The preliminary results have verified the existence of a 1 ms 
time tag discrepancy in the Mini-Mac data reported by others, however, a 2.5 cm unresolved discrepancy 
was found on the short TI to Mini-Mac baseline at Mojave. When the discrepancy in GPS time was 
corrected, several meter RMS orbit differences were found,-however, the 1400 km test baseline between 
nsala and Tromso demonstrated less than two parts in 10 8 agreement between the two cases. The results 
re S0 rt ^ nonstrate d die existence of a change in Tromso antenna position from the location previously 


INTRODUCTION 


eginning in 1987, the National Geodetic Survey has been active in the development of the CIGNET. The 
site locations as of November 1988 are illustrated in Fig. 1 [CSTG Bulletin, Vol. 1, No. 3, 19881. All site 

excep ‘ T s “ kuba - consisted of TI-4100 receivers using the CORE software developed by the 
Applied Research Laboratory of the University of Texas (UT). 

In January, 1989, Mini-Mac 2816AT receivers were installed at Mojave, Westford and Richmond. For 
vera 1 weeks after this installation, the TI receivers were operated simultaneously with the Mini-Mac. In 
ttion, NGS performed local surveys to determine the position vector of the Mini-Mac antenna with 
respect to other geodetic markers in the vicinity of vety long baseline interferometry ( VLBI) systems. 

al 0b T^! d fr ° m mi ™ d receiver have tee" repo" 6 * 1 by (Gurtner et al., 1989), (Rocken et 

ai. 1989) and (Chin, private communication, 1989). In the former case, a 1 ms time tag discrepancy with 
GPS time was observed when TI and Mini-Mac data were mixed. The source of the discrepancy was 
us^d l ° thC Mim MaC software: how ever, the discrepancy is not apparent when only Mini-Mac data were 


Because of the observed problems with mixing receiver types, the purpose of this study was to 
investigate the compatibility between TI and Mini-Mac in the use of CIGNET data for GPS orbit and 
baseline determinations, especially for longer baselines than those used in the previous studies. It was 
believed that this study was essential for the understanding of results obtained using CIGNET data after the 
TI receivers were removed and to aid in establishing historical links to prior campaigns using only TI 
receivers. The results presented in this paper are preliminary, and further detailed studies are planned. 


DATA 


Jamiarv'^O^lQRO Kjif “!* 7! rcceivcrs M Mo i ave > Westford and Richmond began on 

January 30, 1989, m GPS Week 473 using Mim-Mac receiver software Version 1.49. The TI operations 

continued for approximately two weeks at Westford and Richmond after the Mini-Mac installation. At 
Mojave, the UT receiver operated until June. At all three sites, the displacement between the TI and Mini- 
Mac antennas was about 5 m or less. 

_ P ° r ‘ h f P relmmar y analysis reported in this paper, only data from the first three days of simultaneous TI 

Ae Ss? T use f (Janu fy 30 10 February 1, 1989). Since prior wotk had suggested that 
data from Richmond would not be usable, these data were initially excluded from the analysis. 

R.VhmonH fU ™ er 1 " vesUgaUon sl J owed existence of a significant amount of usable TI data from 

n 11 f SC WCTe USCd in SOme ex P eriments . along with TI data from Kauai, Yellowknife, 
Wettzell, Onsala and Tromso. 


PROCEDURE 

S7L an JT Min i‘^ aCd xf 31 ^ tfliee US - sites were used *“ two separate, but complementary analyses. 
First, the TI and Mini-Mac data were used for short baseline tests (<5 m) in which the TI to Mini-Mac 

vector baseline was determined at each of the three sites. The purpose of this test was to assess the 
reported time tag discrepancy and phase center of the respective antennas. Because of the very short 
ase ines, L\ only solutions were obtained, as well as ionospherically corrected solutions. For these 
analyses the orbits were fixed to a nominal ephemeris obtained from either pseudo-range or double- 
ditrerence phase data. All baseline solutions were obtained with phase data 
s econ d , the GPS orbits were determined from double-difference phase data using a) only TI data and b) 
Mim-Mac data from Mojave, Westford and Richmond plus-TI data from the other sites. In these cases, the 
typical estimated parameter set consisted of the epoch position and velocity for each satellite in the three- 
ff 3 (j arc * 3 scale P arameter for the ROCK4 radiation pressure (Fliegel et al., 1985), a parameter for y-bias 
torce, three-hour troposphere zenith delay, appropriate biases and selected antenna coordinates. No a 
priori constraints were used on any estimated parameters. Hie modeling of the GPS force and kinematic 
charactensucs and the site coordinates are given by (Schutz etal, 1989a). In all cases, Tromso was treated 
as a sue with unknown coordinates. 


SHORT BASELINE TESTS 

** Sh ° rt b3Seline between ** 11 ante,ma 3,1(1 ^ Mini-Mac antenna at each 
5L Mi?S a U f S, T' L hese 4031x865 verified ^ exislence of the I ms time tag correction required for 
at Mojave^ " addlUOn ’ a more detailed ao^ysis was performed with the short baseline estimates 

The orbit used for this analysis was derived from double difference phase data using only TI instruments 
(Case B m the next section). This orbit was fixed and the coordinates of the Mojave Mini-Mac were 
estimated with respect to the TI. The TI antenna at Mojave was the FRPA-2 and further information on the 

Z u -\ 8 , 1Ven ,n . the CSTG GPS Bulletin for May-June, 1989. The results are shown in Table 1 

h Mmi-Mac position vector is given with respect to the FRPA antenna used on the TM100. 
The Mini-Mac coordinates were obtained from the average of three individual solutions, each comprising a 
l 6 3X 0 0(1 e differences between the TI and the Mini-Mac. This result was compared with the 
vector difference of positions for the respective antennas obtained on previous surveys. The resulting 2.5 
cm vertical discrepancy is under investigation. No significant time tag adjustment to the Mini-Mac was 
obtained after application of the 1 ms discrepancy. Furthermore, the estimate was obtained with and 
without bias-fixing; however, both cases gave the same result to within a few millimeters. 


ORBIT AND LONG BASEUNE TESTS 

Several cases were used for the orbit and long baseline tests. These cases were: 

A. TI only data from Mojave, Westford, Richmond, Kauai, Wettzell, Onsala and Tromso, 

B. Same as Case A except Richmond TI data were excluded, 

C. Mini-Mac data from Mojave and Westford, no Richmond data and TI data from other sites. 

The orbits that resulted from Case B and Case C are compared in Fig. 2 for PRN-9 and a tabular summary 
for all satellites is given in Table 2. For comparison, the tabular summary for the differences between Case 
A and Case B are shown in Table 3, showing the influence of Richmond TI data. Comparison of Table 2 
and 3 indicates agreement between the three cases at the few meter level 
Since Tromso cannot yet be linked directly into a VLBI coordinate system (the mobile VLBI occupied a 
sue at Tromso in July, 1989), Tromso has been treated as a site with unknown coordinates. The results 
obtained from the three cases considered are shown in Table 4, where the coordinates are given with 
respect to Onsala. The table also contains the number of double-difference measurements and the RMS of 
the global set of measurement residuals. All three cases give Tromso coordinates that agree to within 1.6 

cm in ( x , y, r) and 1.5 cm in the 1400 km baseline length, thus suggesting internal consistency at the level 
of apart in 10 8 . 7 

Comparisons of the Onsala to Tromso baseline with previously determined results are shown in Table 5. 
The previous results are given by (Schutz et al., 1989a; 1989b). While the previous results agreed with 
each other within about 8.5 cm in all components and in the baselines, the Week 473 results disagreed with 
the others by about 81 cm in height. This discrepancy prompted further investigation to identify the source. 
(Engen, private communication, 1989) indicated that a new pillar far the GPS antenna had been 
constructed at Tromso in November, 1988. The results from earlier weeks given in Table 5 were obtained 
with data collected prior to November, 1988, whereas die Week 473 results of this paper used data 
recorded after the pillar was installed. Detailed measurements of the pillar by the Norwegian Mapping 
Agency showed a change in height of about 81 cm, thereby justifying the conclusion that the Tromso 
iscrepancy was caused by a change in the antenna height. This result is of further interest since the initial 
GPS result was obtained without knowledge of the antenna height change, thereby providing a "blind test" 
An additional test of the GPS determination of Tromso will be possible when the VLBI result is available 
at the end of 1989. 


SUMMARY AND CONCLUSION 

Data collected for a period of three days using the TM100 and the Mini-Mac 2816AT receivers collocated 
at Mojave and Westford were processed along with data from other CIGNET sites. Three-day orbit 
comparisons indicate that the impact of replacing TI-4100 with Mini-Mac on the orbits was less than four 
meters RMS. However, the effect on the Onsala to Tromso baseline determination was insignificant A 
previously reported Mini-Mac one millisecond time tag discrepancy from GPS time was verified in the 
process of determining short (-5 m) baselines between the TI-41(X) and the Mini-Mac antennas collocated 
at Mojave and Westford, however, an unresolved 2.5 cm discrepancy exists in the Mojave result The 
baseline between Tromso and Onsala (-1400 km) was determined with a repeatability between three cases 
of 1 part in 10 . An 81 cm change in antenna height at Tromso (due to new construction) was detected 
without any prior knowledge of the change. While this test can be regarded as a "blind test," an essential 
component of the test will be provided by the VLBI determination. Based on the analysis of three days of 
data, it can be tentatively concluded that TI-4100 and Mini-Mac 2816AT CIGNET receivers are 
sufficiently compatible for orbit determination to support a part in 10 8 baseline determinations. However, 
orbit differences at the few meter level require further investigation. 


Acknowledgement. The data were preprocessed by D. Carroll, L. Burton and H. Rim at the Center for 
bpace Research. The information provided by G. Mader and M. Chin at NGS for this study is gratefully 
acknowledged. This research has been supported by NASA. 


REFERENCES 


Gunner, W.G., Beutler, G., and Rothacher, M., (1989), Combining GPS observations made with different 
receiver types. Proceedings of the Fifth International Geodetic Symposium on Satellite Positioning, Las 
Cruces, New Mexico, March 13-17. 

Fliegel, H F., Feess, W.A., Layton, W.C., and Rhodus, N.W., (1985). The GPS radiation force model. 
Proceedings of the First International Symposium on Precise Positioning with Global Positioning 
System, NOAA, Rockville, Maryland, April 15-19. 

Schutz, B.E., Tapley, B.D., Ho, C.S., Rim, H.J., and Abusali, P.AM., (1989a), GPS orbit determination: 
experunents and results, Proceedings of the Fifth International Geodetic Symposium on Satellite 
Positioning, Las Cruces, New Mexico, March 13-17. 

Schutz, B.E., Ho, C.S., and Bevis, M„ (1989b), Analysis of southwest pacific campaign data: July 1988, 
Proceedings of the Fifth International Geodetic Symposium on Satellite Positioning, Las Cruces, New 
Mexico, March 13-17. 

Rocken, C. and Meertens, C.M., (1989), GPS antenna and receiver tests: multipath reduction and mixed 
receiver baselines, Proceedings of the Fifth International Geodetic Symposium on Satellite Positioning, 
Las Cruces, New Mexico, March 13-17. 


Table 1. Mojave short baseline test (TI to Mini-Mac). 



L i Average of three days* -1.197 -2.478 

RMS of three days 0.001 0.001 

Independent solution** -1.181 -2.466 | -4.196 [ 5 001 

*The three days in Week 473 were days 30, 31 and 32 
**The independent solution was obtained by TI to TI and Mini-Mac 
to Mini-Mac local surveys at different times 
(GPS Bulletin, Vol. 2, No. 3, May-June 1989) 

Difference (m) between the three-day and the independent solu tion 
North I East P 

0-0013 -0.0054 0.0239 































Table 5. Onsala to Tromso baseline comparison with other results. 


Components and length (m) 



A z 

Ay 

A z 

l 

Week 473 
Case B 

-1267718.149 

9692.552 

608406.661 

1406187.014 

Week 419 
1 week soln 

-1267718.396 

9692.560 

608405.885 

1406186.900 

Week 420 
1 week soln 

-1267718.417 

9692.533 

608405.862 

1406186.909 

Week 444 
1 week soln 

-1267718.414 

9692.512 

608405.882 

1406186.915 


Week 444 
3 day soln 
w/o S. Pacific 


-1267718.480 


9692.566 


608405.899 


1406186.983 












Figure 1. CIGNET Tracking Sites [CSTG GPS Bulletin] 
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Figure 2. Comparison of solution trajectories in J2000 system for PRN-9 
(Case C - Case B) Epoch - Jan.29,1989 00:00:00 



APPENDIX X 


"Observations of selective availability and the effect on precision geodesy," presented at the 
Fall Meeting of the American Geophysical Union, San Francisco, Calif., December 1990 
(Abstract; Figures) 
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wouM Kara cauaad auddan krwarta* of land in ooaatal 
Waahinfton, Tha foaaU traaa and graaa add to 
pro victual? roporUd avldanoa that auch aarlhquahaa did 
occur, av«n though non# hava happanad In tha 200 
raara alnca Europaano arrived in tha Pacific NorthweaU 


G12B CA: 303 Mon 1440h 

The Impact of Selective Availability (SA) 
on GPS Geodesy 

Presiding: W Melbourne, Jet Propulsion 
Lab; K Feigl, MIT 

G12B-1 1440) INVITED 

An Overview of CPS Selective Availability 

1-P.Yunck (Jet Propulsion Laboratory, California Institute of 
Technology. Pasadena. CA 9! 109) 

Selective Availability, or S A, is now active nearly full time on 
the GPS Block 11 satellites, ten of which are due to be in service 
by December of 1990. S A is designed to degrade the real time 
positioning accuracy itivnable by those not authorized to use 
the GPS Precise Positioning Service. Notably, it is within the 
constraints of unauthorized use that a number of scientific 
groups are now achieving geodetic accuracies ofl : HP or better 
over distances ranging from hundreds to many thousands of 
kilometers — undoubtedly the most precise GPS mea surem ents 
being made in any conditions. 

There is nonetheless a good deal of concern over SA among the 
unauthorized, a situation that has been encouraged by die rather 
elusive manner in which S A has been officially ponrayed to the 
outside world. This presentation will attempt io address com- 
mon SA questions by reviewing GPS signal structure, defining 
the two methods (data message corruption and “dithering”) by 
which SA errors are imposed, giving quantitative examples of 
the effect of SA on point positioning and baseline determina- 
tion, describing strategies by which SA effects can be mini- 
mized. and illustrating these topics with data and results taken 
from recent experiments. 


G12B-2 150R) INVITED 

Global Positioning Syates (GPS) Policies 

Lt Col Jules G. McNeff, USAF (Military 
Assistant to the Assistant Secretary of 
Defense for Command, Control, 
Communications and Intelligence, 

3Dl7 *' Washington, DC 20301- 
3040, 703-695-6123) (Sponsor; Dr. 

William Melbourne) 

This presentation will cover the overall 
process for radionavigation system plan 
and policy development which culminates 
in the biennial Federal Radionavigation 
Plan. It will highlight current plans 
and policies affecting the Global 
Positioning System (GPS). These include 
availability of and access to GPS 
services, implementation of selective 
availability and anti-spoofing features 
and requirements for access to precise 
GPS capabilities. 


G12B-3 15151) INVITED 

Ethical Considerations for the Custodians of GPS 

Roger Bilham, Department of Geological Sciences and ORES. 
University of Colorado, Boulder CO 80309. 

Since its inception, the military architects and custodians of the 
Global Positioning System (GPS) have maintained a consistent 
policy concerning the availability of the precise positioning code 
to the civilian community. "Access to the code can and will be 
denied". Where alternative methods exist this policy has 
moderated the wholesale adoption of CPS geodesy, because of the 
financial risk of investing in the acquisition of commercial GPS 
receivers that may become unusable in the near future. Where 
alternative methods are not available because of the unique 
qualities of GPS positioning, the possible denial of future 
availability continues to alarm the scientific community. Several 
ongoing crustal deformation programs have as a fundamental 
goal the mitigation of seismic and volcanic hazards. The denial 
of the precise positioning code to such programs raises ethical 
and moral issues that were unforeseen at the time that selective 
availability was being designed into the CPS. 


A policy to permit continued access to the GPS precise 
positioning code In its unencrypted original form would provide 
two benefits to the custodians of the GPS: it would be perceived 
as a generous contribution by the US to hazard mitigation In 
many nations of the world, and its free availability would inhibit 
the proliferation of alternative navigation systems by other 
nations who for commercial reasons would not wish to include 
limited access. It is clear that the GPS precise positioning code has 
enormous benefits to civilization beyond any military 
application. If the AS and SA limited access signals are 
Implemented and cannot be successfully overcome in geodetic 
applications, it may be desirable for the world scientific 
community to launch its own satellite system for geodetic control 
In sea level, earthquake and volcano studies. 


G12B-4 153ft 

Monitoring Si A and Hi Effect on GPS Geodesy 

P QCkBl C. tlM * mi CM. Meeneni"** 

> UNAyCQ/ClRES. CB 216 Univ. of Colorado, Boulder CO 80309 

80309 W ° e P t of G* 0 * * Geophys., Univ. of Utah, Salt Lake City. UT. 841 12 . 

S3 ? H^r 10 7°' 1 , H |5 r doubIe differences and with an uncertainty 
ohreroJ^fJ? Smg C dlffcr ® nces lf a rubidium oscillator is used. We 
observed frequency offsets due to S/A of +/- 3 6 Hz from the nominal 

25 (See figure bclow) - Monitoring S/A may aid 

1 h i‘ gh * cune y P° ,m positioning applications Our 
limited data set -Iso show, that S/A monitoring is generate, re' 
muucu on any terrestrial baseline for h>gh accuracy relative GPS k- 
JSJJJ: A*/ 1 * “ obu,ncd from ttrobaseline tests Indicate that ci<5c 
dithering causes maximum double-difference phase residuals (dd ) 
as a function of receiver clock synchronization error of; dd - 004 

nnok^ t"^,L fuOCt T of bucline lengTh of; dd ^ - 0014 cm/ 
(100 km). Since these ire linear relations for maximum observed val- 

MMUefft aiiO Ou*™! V* L° & r’ *** cffcci on ,h « Inline will be 

ihS *£ ,ne “P^don of 100 km. We con- 

C ° Ck dl,henn * cffec « can be monitored and corrected 

wwStiiSl'if m^rp?^ nven,en,ly nt * lected *" double-difference 
j*ocess»ng if the GPS receiver, remain synchronized to better than 10 



G12B-5 15451) 

A Schama for Reducing tha Effact of Salactiva Availability 
on PraclM GPS Carrtar Phaaa Maasuramanta 
KUBILfBSL ROBERT W. KWG, THOMAS A HERRING. AND MARKUS 
ROTHACHEfT; MIT 54-610. Cambridge MA 02139; 

(617) 253-8872 (* visiting from UnNarstty of Bama) 

In early 1990, the Department of Defense began dithering the carrier 
frequency transmitted by Block II uteflftes of the Global Petitioning System 
(GPS). We can charactertza the dithered oscillator frequency by 
differentiating tha phase received at station, equipped with atomic 
frequency standards, it appears that the deviations from tha nominal 
frequency of 1.5742 GHz are typically 1-2 Hz or several parts in 10®. 
Moreover, these deviations vary rapidly with time, often reversing sign 
wltNn a span of 5 to 10 minutes. H left unmodeled, the deviations produce 
quati-periode arches in one-way residuals with amplitude of order 100 
C 7 C *®*- Fof data sampled at the same nominal receiver time, the w"p Mud e 
ol tha arches in (Sortie differences Is small, less lhan 0.02 cycles (or 
transcontinental baselines. For asynchronously sampled data, however, 
these arches can map into the double deferences with an arroiituda ot up 
to several cycles. The magnitude of the problem is roughly proportional to 
the difference in sampling times between the receivers. For example, the 
extreme, but typical case ol a Minimac 2B6AT (sampling at tha even GPST 
minute) and a TI4100 (sampling 0 92 seconds earlier), produces doubly- 
ditferencad residuals with an RMS of about 0.6 cyclas ("broadcast" in 
Fifrjre). To estimate the frequency deviation from the nomteat value, we 
differentiate the phase sampled at 30 second intervals by the CIGNET 
receiver at Mojave At Ms station, tha local ©senator is a hydrogen maser, 
but we expect that any atomic standard stable at the level of ter 1 1 would 
be adequate ter this purpose For convenience, we parameterize the 
deviation of the sateile oeeftator using the tame three coeflttents (ato, 
■It and aa)ai broadcast by fttMteWte In *, navtgafton message. Using 
one set of these vateet every X seconds, we can reduce the RMS in 
asynchronously sampled dortty-dffsrenced residuals by a (actor ol 8. to 
01 cycles rmodetsd* in Figure j. This scatter m UN worse than ter 
Btotft t saMiies wehou Seieetive Availabmty. but nVgrt be improved if ona 
ot the CtQNET stations sampled more frequently Although ou approach 
can modal most oithe effect ol Selective AvtiisMty on asynchronously 
sampled data, a simpler approach tor future surveys it to program the 

receivers to sampts lha phase as dosa to tirro it aneouafy as postible For 

consistency with the CIGNET stations in North America, title means 
sanplng af QA t after the even minute GPST. 
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Figure. Dortty-difte fenced phase residuals usinQ vatees tor the catenae 
frequency as broadcast in the navigation message and as modeled at 
Moiave. The data were transmitted by PRN 16 (Block H) end PRN 3 (Block !) 
■nd received at Mo)ava and Owens Valley. California. The receivers 
sampled at times separated by 0 92 seconds 
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The Impact of Recent Selective Availability on Precise 
Relative Positioning with Static GPS 

EjtocaLUlL Y. Bock 1 , J-B. Minster and D.C. Agnew 
0‘jrxi?’t i 7 >PPS l bls ^ ut30n Oceanography. La Jolla, CA 

Avn^^Ky^A^^Block D dw*y 

P° sltion estimates of the Permanent 
HTf Tfo' (PG OA).imi»Qiifon^ W,cenchZh”1!,^; 

“cStSTs* 3h,?f, b IT °™“* s •* miinuintd between 
S* 8 alterations of the GPS signals do not degrade 
^«wmng accuracy whe n using several hours of doubly-differenced 
1x1 faC t addin * options Of the SA 
wlSi Bto£l lability (even in CaJtforma 

-n.? ^ Wmctncal covcrt « c « quite favorable). 

opcrit,oni1 11 lbou « the same time that the 
Department of Defense initiated SA. Data are collected dailv at 

uKrai^' l Jiww° f ° c . eino « n P h y (SIO) and the Jet Pr^ultion 
JRL SIO^.2 Vf ™ SNR * 8 recciv « deployed^™ 
^ P,flon F,tl °^rvKory, and from three Mini^Mac 
2816AT receivers at the CIGNET tracking station, incTform. 
Massachusetts, and Florida. In the first few week sof ™ ,£ 

SS bwd) "* Mini - Mac «SS?S»5 

, Je - b / 5 S4C0nds - Rogue receiver time tags 
^*te nUy r m ° < ^ , ^ d 10 66 GPS ,imc b «ed. We examin^ 
the repeatability of position determinations of JPL Pirton and SIO 
SS* i?, ,he rtxed CIGNET sites, using five days duri^t^nhS 
pen«J of ume ug discrepancy and four days after the time uu were 
l^n^r^^^tes^efrecu of SA on relative positioning within 

the dlfferei ? t improvement schemes, and arrive at 

However bribed in the firs, paragraph! 

C4m ? d tb,t cumnl SA ha * v «7 adverse effects 
observations, and makes data 
editing more problematic, even with P-code based receivers We also 

^!- Cffe ^ U SA iterations in the navigation message by 
comparison with orbit improved solutions. message by 
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OBSERVATIONS OF SELECTIVE AVAILABILITY AND THE 
EFFECT ON PRECISION GEODESY 

B. E. Schutz (Center for Space Research, The University of Texas at 
Austin, Austin. Texas 78712) 

Y. C. Chao. R. Gutierrez, and C. Wilson 

From July 13-JuIy 22. several experiments with GPS receivers 
were conducted on zero length baselines at McDonald Observatory 
(TX), Planeville (CO), and Quincy (CA) as part of a GPS campaign 
m West Texas. New Mexico, and Nevada. TI-4100 and Tnmble 
4 000 ST receivers were used on the zero length baselines with 
recording epochs that differed by 0.92 s. The McDonald 
Observatory TI was connected to a cesium oscillator, while all other 
receivers used internal oscillators. Evidence of selective availability 
is apparent in the double-differenced residuals; for example, the 
differences between TRN 3 and 20 exhibit periodicities within rms 
of about 10 cm, a characteristic not observed when only Bloc) . 
satellites are used. Similar characteristics were observed at ti 
different sites. The nature of the periodicities and techniques f 
their removal are presented, and the influence of these characteristics 
on long baselines is examined. 


G12B-8 163ft 

The Effect of SA on Baselines from the California GPS 
Array in Southern California 

Frank H . Webb . Ulf J Lindqwister, Geoffrey Blewitt. and James 
Zumbcrge (All at; Jet Propulsion Laboratory. California Institute of 
Technology, Pasadena, CA 91 109; 818-354-4670) 

GPS observables are routinely collected from an array of 
continuously operating Rogue receivers in southern California. The 
receivers are located at Pinyon Flat in the Santa Row Mountains, 
Scnpps Institution of Oceanography in La Jolla, and at JPL in 
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RESIDUALS, PRN 20, Week 549, 167-168 Hours 





Double Difference Residuals Showing SA-Induced Effects 
Caused by Nonsynchronized Receivers 




